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Preface 


An  understanding  of  cathode  sheath  effects  in  gas  discharges  is 
necessary  for  the  design  of  efficient  high-energy  plasma  switches.  Such 
switches  are  essential  for  the  rapid  switching  of  high-energy  inductive 
energy  stores.  Experimental  investigations  of  this  type  of  switch  are 
being  conducted  by  Dr.  Peter  Bletzinger  of  the  Aero  Propulsion  Laboratory, 
Energy  Conversion  Branch.  The  theoretical  investigations  of  the  present 
study  were  done  in  support  of  his  work  with  such  plasma  switches. 

A  number  of  people  provided  guidance  and  encouragement  to  me  during 
the  course  of  this  study.  In  particular,  thar <s  are  given  to:  Dr.  Peter 
Bletzinger,  who  helped  me  realize  that  numbers  and  diagrams  without  physi¬ 
cal  insight  are  meaningless;  Lt.  Col.  William  Bailey,  who  insisted  on 
consistent  boundary  conditions  and  also  helped  me  find  reality;  Capt. 

*  Gary  Duke,  who  provided  empathy  and  consolation  through  his  analytic 

*  * 

investigation  of  a  similar  problem;  Capt.  Greg  Schneider,  who  listened 
critically  and  spoke  helpfully;  and  especially  my  wife  Cheryl  and  children 
Christopher  and  Joshua,  who  helped  me  establish  priorities  and  keep  my 
perspective  throughout  the  ordeal. 

Marc  R.  Hallada 
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Abstract 

Lowke' s  theoretical  investigation  of  gas  discharges  sustained  by 
an  external  ionization  source  was  extended.  In  addition  to  the  processes 
of  electron  diffusion,  electron-ion  recombination,  and  ionization  of  the 
ground  state  considered  by  Lowke,  electron- impact  ionization  of  metastable 
states  and  secondary  emission  from  the  cathode  were  included  in  the  pre¬ 
sent  analysis.  Current-voltage  (I-V)  characteristics  were  obtained  for 
discharges  in  argon  and  methane.  At  current  densities  of  lA/cm,,  Town¬ 
send  ionization  in  the  cathode  sheath  resulted  in  positive-ion  number 
densities  which  significantly  exceeded  those  in  the  positive  column.  I-V 
characteristics  obtained  revealed  a  sharp  ris*  in  the  current  at  an 
"ignition"  voltage.  The  slope  of  the  high- current  region  of  the  char¬ 
acteristic  increased  for  an  increase  in  secondary  emission.  With  the  addi¬ 
tion  of  metastable  ionization,  the  ignition  voltage  and  the  entire  I-V 
characteristic  shifted  to  lower  discharge  voltages.  Below  the  ignition 


voltage  and  at  low  external  ionization  source  strengths,  an  unstable 
region  of  negative  differential  conductivity  was  observed.  _ 


*J.  J.  Lowke,  D.  K.  Davies,  J.  Appl.  Phys.,  ^8(12),  499  (1977). 
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EXTERNALLY  IONIZED  GAS  DISCHARGES 


I.  Introduction 


Background 

High-energy  switches  with  short  switching  times  and  long  operation¬ 
al  lifetimes  are  essential  in  many  areas  of  current  scientific  interest. 
In  particular,  the  realization  of  both  controlled  fusion  and  directed 
energy  weapons  will  require  such  switches.  Although  spark-gap  type 
switches  have  been  used  for  rapid  switching,  -.hey  cannot  function  at 
high  repetition  rates  over  an  extended  period  of  time.  In  addition, 
spark-gap  switches  cannot  switch  "off"  high  energies  (Ref.  1:1).  Such 
"off"  switches  are  required  for  inductive  energy  storage  schemes,  which 
are  very  attractive  in  many  high-energy  applications.  Inductive  energy 
storage  is  much  more  compact  (and  thus  more  mobile)  than  capacitive 
energy  storage.  As  a  consequence,  the  Air  Force  is  particularly  Inter¬ 
ested  in  the  development  of  high  energy  "off"  switches  for  use  in  the 
power  supplies  of  directed-energy  weapon  systems. 

One  of  the  most  promising  approaches  to  producing  such  "off"  switch¬ 
es  involves  the  use  of  an  external  electron  beam  to  produce  ionization 
within  a  gas  discharge.  The  gas  discharge's  operating  parameters  are 
selected  such  that  the  electron-beam-produced  ionization  is  essential 
to  maintain  the  discharge.  As  a  consequence,  the  external  electron  beam 
can  be  used  to  control  the  zmes'  flowing  through  a  gas  discharge.  The 
discharge  acts  as  a  volumetric  conductor  only  while  the  electron  beam 


irradiates  the  active  region  of  the  discharge  (Ref.  1:2).  A  discharge 
Which  exhibits  this  type  of  behavior  Is  known  as  an  externally-sustained 
gas  discharge. 

Literature  Review 

Although  external  sources  of  ionizations  have  been  widely  used  to 
control  gas  discharges,  their  application  to  the  area  of  high-energy 
switching  is  relatively  recent.  A  large  portion  of  the  work  concerning 
external  ionization  has  dealt  with  applications  to  high  pressure  dis¬ 
charges  in  gas  lasers  (Refs.  2:3:4).  Consequently,  early  switching  ex¬ 
periments  made  u;e  of  the  electron  guns  and  gas  mixtures  associated  with 
gas  discharge  lasers.  These  high  current  density  electron  guns  allowed 
gas  mixtures  with  significant  quantities  of  atomic  or  molecular  species 
with  large  electron  affinities  to  be  used.  Taese  "attaching"  species, 
in  turn,  produced  very  short  switch  rise  and  fall  times.  However,  in 
these  experiments,  the  energy  expended  by  the  electron  beam  was  usually 
of  the  same  magnitude  as  the  energy  switched  (Ref.  5:1). 

More  recently,  electron-beam-sustained  switches  with  much  higher 
gains  (ratio  of  energy  switched  to  switching  energy)  have  been  investi¬ 
gated,  primarily  experimentally,  by  several  authors.  Koval’chuk  et.al. 
reported  (1976)  achieving  gains  as  high  as  100  by  decreasing  the  electron 
beam  current  density  (i.e.  the  external  ionization  source  strength)  in  a 
nitrogen  gas  discharge  (Ref.  6).  In  1980,  Bletzinger  experimentally 
Investigated  several  gas  mixtures  (Ar,  Ar  /Oj,  and  CK^)  at  these  lower 
electron  beam  current  densities.  He  observed  a  discontinuous  transition 
to  a  high  conductance  region  as  the  discharge  voltage  was  increased 
beyond  a  certain  value.  Bletzinger  found  methane  to  be  one  of  the  most 
promising  switching  gases  (achieving  gains  as  high  as  1000) ,  because  of 
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its  low  "on-voltage"  and  high  dielectric  strength.  He  also  discovered 
that  methane's  conductivity  and  current  gain  were  increased  by  the  addi¬ 
tion  of  argon  (Ref.  5) .  The  discontinuous  transitions  in  the  current- 
voltage  (I-V)  characteristics  obtained  by  Bletzinger  were  also  observed 
by  Averin  et.al.  in  a  concurrent  investigation.  Averin  maintained  that 
their  experimental  parameters  allowed  them  to  "trace  out"  this  discon¬ 
tinuity,  which  they  termed  a  region  of  "negative  differential  conductivity" 
(Ref.  7). 

Several  theoretical  investigations  of  a  general  nature,  dealing  with 
electric  gas  discharges  sustained  by  various  cources  of  ionization,  have 
also  been  conducted.  A  very  thorough  study  of  discharges  with  external 
ionization  produced  by  photoionization  of  the  gas  was  conducted  by  Ward- 
law  and  Cohen  in  1972.  They  analytically  investigated  the  properties  of 
such  discharges  as  a  function  of  applied  electric  field  and  ionization 
source  strength.  Although  they  included  losses  due  to  diffusion,  they 
neglected  recombination  losses  because  their  ionization  source  strengths 
were  quite  small  and  thus  their  charged  particle  densities  were  also  very 
small  (Ref.  8).  Previously,  in  1958,  Ward  calculated  the  properties  of 
discharges  sustained  by  photoelectric  emission  from  the  cathode.  He 
found  that  at  high  current  densities  a  maximum  was  produced  in  the  posi¬ 
tive-ion  nymber  density  in  the  cathode  sheath  (Ref.  9).  Although  such 
positive  space  charges  were  previously  observed  in  self-sustained  dis¬ 
charges,  Ward  was  the  first  to  calculate  such  an  effect  in  externally- 
sustained  discharges. 

This  peak  in  the  positive-ion  number  density  was  also  observed  in 
later  theoretical  investigations  of  discharges  sustained  by  uniform 
sources  of  ionization.  In  addition  to  finding  such  a  "peak,"  Zakharov 


et.al.  (Ref.  10)  also  calculated  Che  "negative  differential  conductivity" 
region  later  seen  in  experiments  by  Averin,  Bletzinger,  and  Koval' chuk 

(Refs.  7,  5,  6).  Another  computational  investigation,  by  Lowke  and 

\ 

Davies  (Ref.  11),  was  reported  in  1977.  They  also  found  the  positive 
ion  density  maximum  in  the  cathode  sheath,  although  it  only  occurred  for 
the  highest  current  densities  they  considered.  (Their  currents  were,  in 

if 

general,  much  lower  than  those  considered  by  Zakharov.)  The  rise  in  cur¬ 
rent  theoretically  predicted  by  Lowke  and  Davies  was  not  as  rapid  as  that 
observed  experimentally  by  Averin,  Bletzinger,  and  Koval 'chuk.  However, 
Lowke  and  Davies'  analysis  did  not  account  for  two  potentially  important 
processes:  secondary  emission  of  electrons  from  ion  bombardment  of  the 
cathode  and  ionization  of  metastables  in  the  gas  (Ref.  11).  A  very 
thorough  theoretical  analysis  of  externally-ionized  gas  discharges 
(Including  secondary  emission  from  the  cathode)  was  reported  by  Alek¬ 
sandrov  et.al.  in  1978.  They  analytically  modeled  an  externally-sustained 
discharge  using  asymptotic  expansions  in  each  of  several  longitudinal 
regions.  Although  their  analysis  did  not  consider  diffusion,  it  did 
include  electron-ion  recombination,  Townsend  (impact)  ionization,  and 
secondary  emission.  Aleksandrov  attempted  to  describe  the  causes  of  the 
observed  rapid  current  rise  by  identifying  the  dominant  physical  pro¬ 
cesses  in  several  regions  near  the  cathode  (Ref.  12:143).  Both  Alek¬ 
sandrov  and  Lowke  dealt  with  the  effects  of  external  ionization,  recom¬ 
bination,  and  electron  impact  ionization.  The  present  investigation 
extended  the  calculations  done  by  Lowke  and  Davies  to  include  the  addi¬ 
tional  processes  of  secondary  emission  and  metastable  ionization. 

Problem 

The  objective  of  this  study  was  to  investigate  the  properties  of 
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an  electric  gas  discharge  which  is  maintained  by  a  spatially  uniform, 
external  source  of  ionization.  A  one-dimensional  model  of  a  discharge 
between  two  parallel  conducting  electrodes,  the  anode  and  cathode,  was 
adopted.  Particular  emphasis  was  placed  upon  the  discharge  properties 
in  the  vicinity  of  the  cathode.  In  order  to  carry  out  this  investiga¬ 
tion,  the  system  of  differential  equations  describing  the  physical  pro¬ 
cesses  occurring  in  such  a  discharge  was  solved  numerically  and  compared 
with  experimental  results  and  simplified  analytic  solutions. 

Assumptions 

A  number  of  assumptions  were  made  in  ordc  r  to  make  the  numerical 
solution  of  this  set  of  equations  less  complicated  without  significantly 
sacrificing  reality.  Additionally,  it  was  hoped  that  these  assumptions 
would  ease  subsequent  attempts  to  gain  physical  insight  into  the  meaning 
of  these  solutions.  Two  of  the  assumption  dealt  with  the  physical  environ- • 
ment  of  the  discharge,  but  the  majority  were  concerned  with  limitations 
on  the  types  and  nature  of  the  processes  considered. 

The  externally-sustained  gas  discharge's  transverse  dimension  (per¬ 
pendicular  to  the  anode-cathode  axis)  was  taken  to  be  large  enough  for  a 
one-dimensional  analysis  to  be  appropriate.  Planar  electrodes  perpendi¬ 
cular  to  this  axis  were  the  boundaries  between  which  the  discharge 
occurred.  The  region  between  the  electrodes  contained  a  source  of  posi¬ 
tive  ion  and  electron  pairs  which  was  assumed  to  be  both  spatially  uni¬ 
form  and  constant  for  all  values  of  the  electric  field. 

The  majority  of  the  assumptions  made  are  concerned  with  the  types 
of  physical  processes  allowed  to  occur  in  the  discharge.  Initially,  the 
interactions  of  only  electrons,  neutral  atoms  and  molecules,  and  singly- 
charged  positive  ions  were  considered.  Later,  this  assumption  was 
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relaxed  and  the  Influence  of  metastable  neutrals  was  Included.  Thus,  the 
processes  considered  at  each  point  in  the  Investigation  were  limited  by 
the  types  of  particles  available.  Initially,  the  only  kinetic  processes 
allowed  were  Townsend  (impact)  ionization  and  electron-ion  recombination. 
With  the  inclusion  of  metastable  neutrals,  the  process  of  metastable 
ionization  was  added.  Throughout  the  investigation,  the  diffusion  of 
positive  ions  was  neglected:  because  of  the  relatively  low  fields 
encountered,  the  ion  mobilities  were  much  less  than  those  of  electrons. 

The  final  assumption,  which  affected  nearly  all  of  the  processes 
mentioned  and  greatly  simplified  the  calculations,  was  that  the  discharge 
was  in  equilibrium  with  the  local  sustaining  electric  field.  This  as¬ 
sumption  permitted  the  use  of  tabulated  values  for  equilibrium  trans¬ 
port  coefficients  (local  drift  velocities  and  diffusion  and  ionization 
coefficients)  obtained  from  drift  tube  measurements  or  calculated  by 
solving  the  collisional  Boltzmann  equation.  The  validity  of  this  as¬ 
sumption  is  particularly  questionable  in  the  electrode  regions  since 
the  equilibrium  coefficients  are  characterized  by  uniform  electric  fields 
and  uniform  neutral  number  density  distributions,  requirements  which  are 
certainly  not  satisfied  near  the  electrodes .  Additional  assumptions  were 
also  made  concerning  the  boundary  conditions  at  the  electrodes  and  at 
various  points  within  the  discharge.  These  boundary  condition  assump¬ 
tions  and  the  previously  listed  assumptions  are  explained  in  more  detail 
as  they  are  encountered  in  the  development  of  the  theory  in  the  following 
chapter. 

General  Approach 


Based  upon  the  assumed  quasi-equlllbrlum,  the  gas  discharge  was 
modeled  by  a  set  of  four  coupled  first-order  differential  equations. 


These  equations  are:  the  continuity  equations  for  electrons  and  posi¬ 
tive  ions,  the  electron  current  density  equation,  and  Poisson's  equa¬ 
tion.  In  certain  regions  of  the  discharge,  this  particular  set  of 
equations  becomes  "stiff,"  i.e.  some  terms  are  changing  very  rapidly 
while  others  are  changing  very  slowly.  To  solve  this  stiff  system  of 
equations,  a  computer  program  called  Gear,  written  by  C.  V.  Gear  (Ref. 

25),  was  used.  This  program  uses  an  implicit  linear  multistep  method, 
termed  a  "backward  differentiation  formula"  method.  Using  Gear,  the 
charge  and  current  densities,  the  electric  field,  and  the  voltage  were 
calculated  as  functions  of  position  within  the  inter-electrode  gap. 
Compiling  the  remits  from  a  number  of  these  calculations,  current- 
voltage  (I-V)  discharge  characteristics  were  sroduced  for  several  source 
strengths,  pressures,  and  secondary  emission  coefficients  for  argon  and 
methane  discharges.  The  computer  code  written  to  model  these  discharges 
was  validated  by  comparing  the  calculated  I-V  curves  and  the  spatial 
variation  of  the  field  and  number  densities  within  particular  discharges 
with  those  calculated  by  Lowke  and  Davies  (Ref.  11).  In  addition,  the 
I-V  curves  calculated  for  argon  were  also  compared  with  experimental  I-V 
characteristics  obtained  by  Leffert  et.al.  (Ref.  13).  The  computer  code 
was  then  modified  to  Include  the  effects  of  secondary  emission  of  electrons 
from  the  cathode  and  ionization  of  gas  atoms  in  metastable  states.  This 
modified  code  was  used  to  assess  the  relative  importance  of  these  pro¬ 
cesses  in  various  gas  discharges.  The  results  of  these  calculations  were 
compared  with  experiments  performed  by  Bletzinger  (Ref.  5)  and  theoreti¬ 
cal  calculations  by  Zakharov  et.al.  (Ref.  10). 

Sequence  of  Presentation 


The  physical  processes  considered  in  the  computer  calculations 


performed  are  presented  In  the  following  chapter,  entitled  "Theory." 

The  expressions  used  for  the  various  reaction  rates  and  coefficients  for 
these  processes  are  described  and  any  approximations  made  are  noted. 

More  specific  information  is  given  in  Appendix  A.  In  addition,  the 
classification  and  characteristics  of  gas  discharges,  in  general,  are 
also  presented  in  this  chapter. 

The  physics  described  in  the  theoretical  chapter  is  transformed  into 
a  model  of  a  gas  discharge  in  chapter  III,  "Numerical  Method."  In  this 
chapter,  a  system  of  differential  equations  which  describes  the  inter¬ 
action  of  the  physical  processes  in  a  gas  discharge  is  developed.  The 
boundary  conditions  imposed  on  this  set  of  eqrations  are  also  described. 

The  solutions  to  these  equations  for  a  range  of  discharge  parameters 
are  presented  in  chapter  IV,  "Results."  These  results  are  presented  in 
two  forms,  in  "profiles"  of  the  spatial  variation  within  a  discharge  of 
the  independent  variables  and  in  current-voltage  (I-V)  "characteristics" 
obtained  by  plotting  the  voltage  drop  across  the  discharge  for  a  range 
of  discharge  currents.  These  I-V  characteristics  are  compared  to  similar 
characteristics  which  were  measured  experimentally. 

However,  because  the  calculations  and  experiments  do  not  completely 
agree,  indicating  that  other  physical  processes  may  need  to  be  considered, 
suggestions  for  further  research  are  presented  in  the  final  chapter, 
"Conclusions  and  Recommendations."  Also,  a  brief  summary  of  the  physical 
insights  revealed  by  the  calculations  and  their  implications  for  gas  dis¬ 
charge  switch  design  are  included  in  this  final  chapter. 
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II.  Theo: 


A  gas  discharge  may  be  defined  as  a  gas  with  free  charged  particles 
which  allow  the  gas  to  carry  a  current.  Gas  discharges  can  occur  over  a 
wide  range  of  gas  pressures  and  can  carry  currents  ranging  over  more 
than  twelve  orders  of  magnitude:  from  one  millionth  of  an  amp  to  one 
million  amps.  Gas  discharges  may  exist  only  transiently,  as  in  sparks 
or  lightning,  or  they  may  operate  in  the  steady  state,  as  in  neon  lights. 
In  addition,  the  discharge's  behavior  is  often  greatly  influenced  by  the 
nature  of  the  electrical  circuit  of  which  it  Is  a  part.  In  all  cases, 
however,  since  gases  are  normally  insulators,  charged  particles  must  be 
produced  in  some  way  to  carry  the  current  in  a  discharge.  In  addition, 
an  electric  field  must  be  provided  to  drive  these  charged  particles  in 
order  for  them  to  form  a  current. 

Charge  Carrier  Kinetics 

Cosmic  and  ultraviolet  radiation,  natural  radioactivity,  and  thermal 
collisions  among  gas  particles  are  responsible  for  producing  a  background 
level  of  charged  particles  (~1000  positive  and  negative  ions  per  cm  )  in 
all  gas  samples  under  normal  conditions  (Ref.  15:47).  The  higher  levels 
of  ionization  required  for  efficient  current  conduction  may  be  produced 
by  an  external  source  of  charged  particles  and/or  by  the  imposition  of 
an  electric  field  on  the  gas.  This  electric  field  is  usually  produced 
by  applying  a  voltage  across  two  electrodes  placed  in  the  gas.  The  field 
exerts  a  force  on  the  charged  particles  which  is  parallel  to  the  field. 
The  fluid  equation  of  motion  for  one  species  of  such  particles  is: 

n^ft  +  “  ±enE  -  VP  -  mnvW  (1) 
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Where  p  is  the  momentum,  P  is  the  pressure,  m  is  the  mass,  n  is  the 
number  density,  W  is  the  velocity,  v  is  the  average  collision  frequency, 

and  E  is  the  electric  field.  In  one  dimension,  in  a  steady  state,  and 

\ 

for  negligible  pressure  gradients,  this  equation  reduces  to: 


HP 

0  ■  ±  enE  -  - - nmVW 

dx 


(2) 


Using  the  ideal  gas  law,  the  drift  velocity  may  be  represented  as 


/X. 


A  eNT,  .KT.  dn/dx 
(mV  )E  “  V  n 


(3) 


Where  T  is  the  temperature  and  K  is  the  Boltzmann  constant.  In  this 
equation,  the  coefficient  of  the  electric  field  is  called  the  mobility 
and  the  coefficient  of  the  number  density  gradient  term  is  called  the 
diffusion  coefficient. 


(4) 


(5) 


These  two  coefficients  are  connected  by  what  is  known  as  the  Einstein 
relation: 


K 


±  eD 
KT 


(6) 


Also,  in  both  coefficients,  the  collision  frequency,  v,  for  a  particle 
in  a  gas  of  neutral  number  density,  N,  is  given  by 

v  -  N<va>  (7) 


where  <Va>  is  the  average  of  the  product  of  the  particles'  relative 
velocity  and  the  collision  cross-section.  Thus,  in  the  case  of  negli¬ 
gible  diffusion,  equation  3  may  be  expressed  as 
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(8) 


W  - 


■  -■A.  (I) 
m«jv>  VN/ 


Because  of  the  direct  relationship  between  number  density  and  collision 

frequency,  this  type  of  functional  dependence  on  the  ratio  (E/N)  is  very 

2 

common  in  the  study  of  gas  discharges.  (The  units  for  E/N  are  volt-cm 

-17  2 

or  Townsend,  where  1  Townsend  (Td)  ■  1  x  10  volt-cm  .) 

The  functional  dependence  of  the  mobility  and  the  diffusion  coeffi¬ 
cient  on  the  collision  cross-section,  which  is  generally  a  non-linear 
function  of  E/N,  requires  complicated  analytic  fits  or  tabulated  data  in 
order  to  be  modeled  correctly.  Tabulated  values  for  the  drift  velocity 
and  diffusion  coefficient  for  electrons  in  argon  and  methane  were  obtained 
from  solutions  of  the  collisional  Boltzmann  equation  and  experimental  mea¬ 
surements.  However,  it  was  experimentally  foind  that  the  mobility  of 
positive  ions  could  be  described  independent  of  the  field  strength  (Ref. 
11:4994).  More  detailed  information  on  the  actual  values  used  for  these 
''material  functions"  is  contained  in  Appendix  A. 

The  electrons  and  positive  ions,  whose  behavior  in  a  gas  discharge 
is  described  using  these  material  functions,  must  be  produced  in  some 
manner  for  a  current  to  flow  in  a  gas,  which  is  normally  an  insulator. 

How  these  charged  particles  are  generated  and  lost  in  the  gas  discharges 
investigated  in  this  study  is  the  subject  of  the  next  section. 

Fundamental  Processes 

To  accurately  model  the  behavior  of  a  gas  discharge,  it  is  necessary 
to  consider  the  detailed  interactions  of  individual  gas  particles.  In 
this  analysis,  "average"  properties  are  used  to  describe  particles  of 
the  same  type.  These  average  particles  are  assumed  to  interact  with 
each  other  and  with  the  electrodes  in  ways  which  lead  to  the  loss  or 
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production  of  charged  particles.  Those  processes  which  are  considered 
in  the  generation  of  charged  particles  include:  impact  (Townsend)  ion¬ 
ization  of  neutral  particles  in  the  ground  and  metastable  states  and 

\ 

secondary  emission  of  electrons  from  the  cathode.  Those  processes  con¬ 
sidered  in  the  loss  of  charged  particles  include:  electron-ion  recom¬ 
bination  and  electron  diffusion. 

Impact  Ionization.  A  small  voltage  applied  to  two  parallel  elec¬ 
trodes  immersed  in  a  gas  will  cause  the  electrons  and  ions  naturally  pre¬ 
sent  in  the  gas  to  flow  to  the  oppositely-charged  electrodes  and  thus 
produce  a  current.  At  a  given  applied  voltage,  an  electric  field  is 
produced  between  the  electrodes  which  accelerates  the  charged  particles 
to  the  electrodes.  At  higher  applied  voltages,  the  strength  of  this 
electric  field  is  increased,  accelerating  the  charged  particles,  parti¬ 
cularly  the  electrons  (because  of  their  greater  mobility),  even  more.  At 
voltages  larger  than  a  certain  value,  the  electrons  gain  enough  energy 
between  collisions  to  ionize  neutral  particles  with  which  they  collide 
(see  Table  I) .  The  electrons  produced  may  also  be  accelerated  by  the 
electric  field  and  collide  with  other  particles,  producing  more  ionizing 
collisions.  At  this  point,  the  gas  is  said  to  have  "broken  down."  In 
an  attempt  to  theoretically  describe  the  increase  in  current  which 
results,  Townsend  introduced  the  quantity  a,  known  as  Townsend’s  first 
ionization  coefficient.  This  coefficient  is  defined  as  the  average 
number  of  ionizing  collisions  made  by  an  electron  in  traveling  one  cm 
in  the  direction  of  the  electric  field.  Thus,  the  change  in  the  number 
density  of  electrons,  dng,  occurring  in  a  distance  dx  is  given  by: 

dn  *  n  adx  (9) 

e  e 
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Table  I 


Kinetic  Reactions 


Excitation 

-  p  ^ 

A+e  +  A*  +  e  +  A+  e  +  hv 

Ionization 

A  +  e  A+  +  Townsend  (impact)  ionization 

A*  +  e  -*■  A+  +2  metastable  ionization 

e 


Recon  ?ination 

A"  +  B  +  e-*,A+B 
AB  +  e  -*•  A  +  B 
A+  4-  e  A  +  hv 


three-body  recombination 
dissociative  recombination 
radiative  recombination 


Attachment 

A  +  e  ■*  A  +  hv 
AB  +  e  ■*  A  +  B 
2A  +  e  -*■  A”  +  A 


radiative  attachment 
dissociative  attachment 
three-body  attachment 
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(if  recombination  and  diffusion  are  neglected) .  It  may  be  assumed  that 
an  electron  will  ionize  a  gas  molecule  if  the  energy  it  acquires  between 

collisions  equals  or  exceeds  the  ionization  energy,  eV^,  (V^  is  the  ion- 

\ 

lzation  potential)  of  the  gas  molecules : 


eEx  >  eV±  (10) 

where  x  is  the  path  length  between  collisions  for  an  electron  and  L  is 
the  electric-field  (Ref.  33:514).  Also,  the  average  number  of  ionizing 
collisions  which  will  occur  in  a  given  path  length  is  inversely  propor¬ 
tional  to  the  electron  mean  free  path.  Thus,  Townsend's  ionization 
coefficient  may  be  expressed  as  the  inverse  of  the  mean  free  path,  L, 
times  some  function,  F,  of  the  energy  gained  by  an  electron  (Ref.  33:515). 


a 


F(Ex) 

L 


(ID 


However,  as  mentioned  earlier,  only  those  electrons  whose  path  lengths 
between  collisions  exceed  x  will  have  energies  large  enough  to  ionize  gas 
molecules.  The  fraction  of  these  electrons  with  free  path  lengths 
exceeding  x  is  obtained  from  kinetic  theory  and  is  given  by 


f(x)  -  x/L 


(Ref.  33:515)  (12) 


Thus,  the  number  of  electron  mean  free  paths  which  exceed  x  in  a  given 
unit  length  is  just 

,  -  x/L 

-  £  e  (Ref.  33:515)  (13) 

This  expression,  then,  may  be  used  to  represent  Townsend's  ionization 
coefficient.  Recalling  that  the  mean  free  path  is  inversely  proportional 
to  the  gas  pressure,  and  using  the  expression  for  x  in  equation  10, 


df  (x) 
dx 
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equation  13  may  be  rewritten: 

-  BVi 

a/P  -  A  e  (Ref.  33:516)  (14) 

where  A  and  B  are  constants  determined  by  the  gas  in  the  discharge  and 
P  is  the  gas  pressure.  In  the  present  analysis,  an  analytic  expression 
such  as  this,  given  by  Lowke  and  Davies  (Ref.  11:4994),  was  used  to  cal¬ 
culate  the  Townsend  ionization  coefficient  in  argon.  However,  for 
methane,  Cookson's  experimental  values  of  the  "net  ionization  coeffi¬ 
cient,"  which  included  electron  losses  by  attachment,  were  used  (Ref. 
14:18).  From  tha  appropriate  Townsend  ionization  coefficient,  the  rate 
of  such  ionization,  i.e.  the  number  density  of  electron-ion  pairs 
produced  in  one  second  by  a  single  electron,  is  obtained  from  the  pro¬ 
duct  of  a  and  the  average  drift  velocity  of  an  electron  betweem  colli¬ 
sions,  W  : 

e 

v^l/sec)  *  a(l/cm)  •  We(cm/sec)  (15) 

Metastable  Ionization. 

Townsend’s  first  ionization  coefficient  accounts  for  ionization  from 
the  ground  state  of  an  atom  or  molecule.  Another  means  must  be  used  to 
include  the  effect  of  ionization  from  metastable  excited  states  of  atoms 
or  molecules.  This  metastable  ionization  is  particularly  important  for 
rare  gases,  such  as  argon  and  krypton,  and  molecular  gases,  such  as  (>2 
and  N2.  To  include  this  effect,  the  number  density  of  metastable  neutrals 
and  the  rate  of  ionization  from  these  states  had  to  be  determined.  The 
rate  of  ionization  of  metastables  was  approximated  by  assuming  that  the 
electrons  had  a  Maxwellian  ssizrzy  distribution  and  that  the  magnitude  of 
this  cross-section  was  similar  to  that  of  the  ground  state  ionization 
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cross-section,  although  this  cross-section  was  shifted  to  much  lower 
energies.  The  number  density  of  neutral  argon  metastables  was  approxi¬ 
mated  by  using  a  Boltzmann  distribution,  at  the  electron  temperature, 
between  the  ground  and  metastable  states.  The  effect  of  metastable 
ionization  was  included  in  the  investigation  of  argon,  but  not  in  the 
Investigation  of  methane. 

Secondary  Emission.  One  additional  source  of  charged  particles 
considered  in  this  investigation  was  secondary  electron  emission  from 
the  cathode.  (Although  electrodes  can  also  emit  positive  ions  in  certain 
circumstances,  the  rate  of  positive  ion  emission  is  usually  very  low  and 
was  therefore  neglected  in  this  analysis.)  These  electrons  are  produced 
by  the  bombardment  of  the  cathode  by  positive  ions.  In  a  metallic  elec¬ 
trode,  which  has  "free"  electrons,  an  electron  is  emitted  if  it  is  given 
enough  energy  by  the  incident  ion  to  overcome  the  work  function  (potential 
barrier  at  the  surface)  of  the  material.  This  effect  was  also  described 
by  Townsend,  who  defined  a  secondary  emission  coefficient,  y,  which  is 
the  cathode  yield  in  electrons  per  incident  icn.  The  value  of  y  Is  deter¬ 
mined  by  the  cathode  material,  the  nature  of  the  gas  in  the  discharge, 
and  the  positive  ion  energy,  which  is  a  function  of  the  electric  field 
strength  in  the  discharge.  The  values  for  y  given  in  the  literature  for 
positive  argon  ions  in  the  ground  and  metastable  states  range  from  less 
than  one  percent  to  almost  forty  percent  for  various  values  of  electric 
field  strength  and  cathode  material.  Thus,  a  range  of  secondary  emission 
coefficients,  each  constant  with  respect  to  field  strength,  was  used  in 
this  analysis  (Ref.  19:118). 

Recombination.  Recombination  may  occur  when  a  positive  ion  collides 
with  an  electron  or  negative  ion.  The  free  electron  or  the  extra  electron 
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from  the  negative  ion  recombines  with  the  positive  ion,  producing  one 
or  two  neutral  particles.  The  initial  kinetic  energy  of  the  particles 
and  the  potential  energy  of  ionization  may  be  dissipated  in  one  of  three 
major  ways  (see  Table  I).  The  energy  may  be  radiated  away  (radiative 
recombination) ;  it  may  be  transferred  to  a  third  particle  (an  electron 
or  a  heavy  particle)  as  kinetic  energy  (three-body  recombination) ;  or  it 
may  be  used  to  dissociate  the  neutral  molecule  formed  from  a  molecular 
ion  and  an  electron  (dissociative  recombination) . 

The  probability  of  recombination  decreases  as  the  relative  velocity 
of  the  oppositely-charged  particles  increases.  Consequently,  the  prob¬ 
ability  of  electron-ion  recombination  is  usually  much  less  than  that  of 
ion-ion  recombination  because  electrons  typically  have  much  higher  velo¬ 
cities  than  ions.  However,  in  methane  and  argon,  their  very  small  attach¬ 
ment  cross-sections  result  in  correspondingly  low  number  densities  of 
negative  ions  and  thus  low  rates  of  ion-ion  recombination.  Therefore, 
electron-ion  recombination  was  considered  to  be  the  dominant  recombination 
mechanism  in  both  methane  and  argon. 

By  comparing  typical  rates  for  the  various  types  of  electron-ion 
recombination,  it  was  found  that  dissociative  recombination  was  the 
dominant  recombination  mode  for  the  discharge  conditions  considered  in 

both  argon  and  methane.  The  radiative  recombination  rate  for  electrons 

—13  3 

and  ions  is  estimated  by  von  Engel  to  be  of  order  10  cm  /sec  (Ref.  17: 

141).  In  contrast,  because  of  the  multiplicity  of  energy  absorption 

possibilities  in  molecules,  dissociative  recombination  coefficients  are 
—7  3 

usually  of  order  10  cm  /sec  (Ref.  17:142-143).  Also,  at  the  inter¬ 
mediate  gas  pressures  (240  and  760  Torr)  and  low  electron  number  densi¬ 
ties  investigated,  the  rate  of  three-body  recombination,  with  an 
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electron  as  the  third  body,  is  approximately  10  cm  /sec.  When  the 
third  body  is  a  heavy  particle,  the  rate  is  even  lower,  at  10  car/ 
sec  (Ref.  18:36). 

As  noted  previously,  the  electron-ion  recombination  rate  is  greater 
for  lower  relative  velocities  between  electrons  and  ions.  As  a  result, 
recombination  is,  in  turn,  greater  at  lower  field  strengths.  In  the 
Investigation  of  argon,  the  results  using  various  recombination  coeffi¬ 
cients  were  compared.  Both  the  constant  coefficient  used  by  Lowke  and 
Davies  and  a  coefficient  with  a  functional  dependence  on  electron  tem¬ 
perature,  offered  by  them  as  an  alternative,  were  used  (Ref.  11:4994). 

In  addition,  a  recombination  coefficient  with  a  field  strength  depen¬ 
dence,  given  by  Kline,  was  used  for  comparison  (Ref.  14:13).  The  inves¬ 
tigation  of  methane  made  use  of  a  recombination  coefficient,  also  with  a 
functional  field  dependence,  which  was  calculated  by  Kline  (Ref.  14:19) 
using  measured  recombination  cross-sections.  Detailed  information  on 
the  values  and  analytic  expressions  used  for  these  rates  are  provided 
in  Appendix  A. 

Attachment .  An  electron  which  collides  with  a  neutral  particle  may 
become  "attached"  to  the  neutral,  forming  a  negative  ion,  in  the  process 
called  attachment.  In  order  for  this  process  to  occur,  the  neutral 
particle  must  have  an  electron  affinity  resulting  from  a  nearly  filled 
atomic  electron  shell  or  molecular  orbital.  Strong  attachment  rates  are 
exhibited  by  only  a  few  species  of  particles,  including  atomic  and  mole¬ 
cular  oxygen  and  the  halogen  gases  (e.g.  fluorine  and  chlorine).  Similar 
to  recombination,  the  kinetic  energy  of  the  incident  particles  and  the 
binding  energy  (affinity)  of  the  resulting  neutral  particle  must  be 
dissipated  in  some  way.  Consequently,  attachment  may  take  one  of  three 
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forms:  radiative  attachment,  dissociative  attachment,  or  three-body 
attachment  (see  Table  I)  (Ref.  15:29).  However,  In  this  analysis,  these 
forms  of  attachment  were  not  distinguished  from  each  other  and  attach¬ 
ment  itself  was  only  implicitly  accounted  for  in  methane.  For  methane 
calculations,  an  experimentally  measured  "net  ionization  coefficient" 
was  used  which  considered  attachment  as  a  reduction  in  the  Townsend  ion¬ 
ization  coefficient.  For  argon  calculations,  attachment  effects  were 
completely  neglected  because  of  argon's  inert  nature. 

Diffusion.  Charged  particles  are  also  lost  through  the  process  of 
diffusion.  The  rate  of  this  loss  is  directly  proportional  to  the  dif¬ 
fusion  coefficient 

D  -  —  (16) 

mv 

as  shown  previously  in  the  discussion  of  charge  carrier  kinetics.  In 
this  expression,  T  is  the  temperature  of  the  particle  species,  m  is  the 
particle  mass,  and  v  is  the  average  collision  frequency  of  the  particle. 
Thus,  diffusion  is  greater  for  higher  temperatures  and  smaller  masses,  and 
collision  frequencies  (and  thus  lower  gas  pressures) .  Because  the  mass 
of  an  ion  is  usually  over  2000  times  that  of  an  electron,  the  diffusion 
of  positive  ions  was  neglected  in  this  analysis.  Values  for  the  longi¬ 
tudinal  diffusion  coefficient  of  electrons,  as  a  function  of  E/N,  were 
obtained  from  calculations  (for  argon)  and  from  experimental  measurements 
(for  methane) .  The  sources  for  these  values  are  listed  in  Appendix  A. 

Classification  and  Character! zatlcs  of  Gas  Discharges 

In  a  steady-state  gas  discharge,  the  voltage  applied  to  the  elec¬ 
trodes  determines  how  the  charged  particle  generation  and  loss  processes 
interact  in  the  production  of  a  steady  state  distribution  of  charged 
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particles  in  the  discharge.  This  distribution  of  charged  particles  then 
determines  the  magnitude  of  the  discharge  current.  Discharges  may  be 
divided  into  two  basic  categories  using  the  breakdown  voltage,  mentioned 
previously  in  the  discussion  of  Townsend  ionization,  as  the  dividing 
point.  At  voltages  below  breakdown,  a  discharge  is  said  to  be  non-self- 
sustained,  while,  above  breakdown,  the  discharge  is  self-sustained.  Non¬ 
self-sustained  discharges  require  an  external  source  of  ionization  to 
provide  sufficient  numbers  of  current-carrying  charged  particles.  A 
self-sustained  discharge  generates  enough  current-carrying  charged  par¬ 
ticles,  from  internal  discharge  processes  such  as  Townsend  ionization, 
secondary  emission,  or  thermionic  emission,  to  compensate  for  any  loss 
processes  in  the  discharge. 

Non-Self -Sustained  Discharge.  If  no  ionization  is  produced  by  the 
discharge  itself,  and  the  only  source  of  charged  particles  is  natural 
radioactivity  and  cosmic  rays,  at  very  low  applied  voltages  the  discharge 
current  is  measured  in  random  bursts.  This  intermittent  current  reaches 
a  maximum  when  the  applied  voltage  is  increased  to  the  point  that  the 
field  thereby  created  sweeps  all  of  the  ions  and  electrons  to  the  elec¬ 
trodes  before  they  can  recombine.  (See  Figure  1,  which  shows  the  varia¬ 
tion  of  current  with  voltage  (I-V  characteristic)  in  a  low  pressure  (~1 
Torr)  gas  discharge.  Thus,  if  the  external  ionization  source  produces 

electron-ion  pairs  at  a  rate  S  (1/cm  /sec) ,  between  electrodes  separated 

2 

a  distance  d(cm),  the  saturation  current  density  J  (A/cm  )  is  given  by 

s 

Jg  -  edS  (17) 

where  e  is  the  electronic  charge  in  Coulombs  (Ref.  15:48). 


Current,  amps 


Wich  a  further  increase  in  the  applied  voltage,  Townsend  ionization 
begins  to  occur,  additional  electrons  and  ions  are  produced,  and  the  dis¬ 
charge  current  Increases  from  its  saturated  value.  At  a  certain  voltage, 
the  number  of  electrons  and  ions  produced  by  Townsend  ionization  and 
secondary  emission  are  sufficient  to  replenish  the  number  lost  by  recom¬ 
bination,  diffusion,  attachment,  and  the  external  circuit  current.  At 
this  voltage,  the  discharge  is  said  to  have  "broken  down,"  and  is  thus 
self-sustained . 

Self-Sustained  Discharge.  Breakdown  occurs  at  the  beginning  of  what 
is  known  as  the  Townsend  discharge  region.  This  discharge  is  also  known 
as  a  dark  discharge  since  it  emits  very  little  visible  radiation.  The 
Townsend  discharge  and  current  saturation  regions  are  indicated  in  Fig.  1. 

As  the  discharge  current  increases,  the  voltage  across  the  electrodes 
drops  considerably  as  the  discharge  transitions  to  a  more  efficient  oper¬ 
ating  regime.  This  region  of  the  I-V  characteristic,  where  the  voltage 
falls  as  the  current  increases,  is  known  as  a  subnormal  or  corona  dis¬ 
charge.  (See  Fig.  1.)  In  this  region,  the  effective  electrical  resist¬ 
ance  of  the  discharge  is  negative  (Ref.  15:82). 

As  the  voltage  drops,  the  discharge  begins  to  glow.  At  a  certain 
current,  the  discharge  voltage  stabilizes:  and  remains  relatively  constant 
as  the  current  is  further  increased.  In  this  region  of  the  I-V  char¬ 
acteristic,  known  as  a  glow  discharge,  the  luminous  column  covers  only 
a  portion  of  the  cathode  surface.  As  the  current  in  the  discharge  is 
increased,  this  "cathode  spot"  increases  in  size,  keeping  the  current 
density  constant  (Ref.  15:82).  When  the  cathode  spot  covers  the  entire 
cathode  surface,  a  further  increase  in  current  requires  the  voltage 
across  the  electrodes  to  increase  rapidly  to  a  maximum.  This  region  of 
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the  I-V  characteristic  is  termed  an  abnormal  glow  discharge  (Ref.  15: 

82). 

If  Che  discharge  is  driven  to  higher  currents,  the  discharge  volt- 

\ 

age  drops  by  over  an  order  of  magnitude  as  a  new  mechanism  of  electron 

production  becomes  operative.  At  these  higher  currents,  electrons  are 

-  ■* 

also  produced  by  thermionic  emission  at  the  cathode.  This  type  of  gas 

:p 

discharge,  which  is  known  as  an  arc  discharge,  is  much  more  luminous 
than  a  glow  discharge.  In  addition,  the  current  carried  by  an  arc  dis¬ 
charge  is  determined  primarily  by  the  external  electric  circuit  of  which 
the  discharge  is  a  part,  while  the  current  carried  by  a  glow  discharge 
is  also  very  dependent  upon  the  properties  of  the  discharge  itself  (Ref. 
15:88).  This  difference  occurs  because  the  physical  processes  respon¬ 
sible  for  producing  a  glow  discharge  are  very  sensitive  to  the  discharge 
current  and  voltage,  while  the  thermionic  emission  which  allows  an  arc 
to  form  is  very  stable  once  it  is  initiated. 

Externally-sustained  Discharge.  However,  none  of  the  previous  dis¬ 
charge  descriptions  allow  for  strong  external  ionization  sources.  If 
such  sources  are  considered,  it  is  found  that  although  a  discharge  may 
appear  to  resemble  a  glow  discharge,  it  may  not  necessarily  be  self -sus¬ 
tained,  as  required  by  the  previous  description.  In  fact,  the  discharges 
produced  under  strong  external  ionization  may  not  exactly  resemble  any 
of  the  various  self -sustained  discharges  mentioned  previously.  However, 
since  the  currents  carried  by  glow  and  arc  discharges  are  similar  to 
those  desired  for  electron-beam-sustained  switches,  some  of  the  features 
of  these  discharges  will  be  presented.  These  characteristics  will  later 
be  compared  to  those  predicted  by  this  investigation  of  externally-sus¬ 
tained  gas  discharges. 
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Low  Pressure  Glow  Discharge.  A  low  pressure  glow  discharge  may  be 

divided  longitudinally  into  several  distinct  regions,  as  shown  in  Fig.  2. 

Not  all  of  these  regions  produce  visible  radiation.  Light  is  emitted 

only  in  those  regions  of  the  discharge  in  which  excited  ions  or  neutrals 

transition  to  lower  energy  levels  or  in  which  electrons  and  ions  radia- 

tively  recombine.  The  major  features  in  each  of  four  longitudinal  regions 

if  ' 

of  a  glow  discharge  are  described  in  the  following  sections. 

Cathode  sheath.  The  electrons  required  to  maintain  the  discharge 
are  produced  primarily  by  positive-ion  bombardment  of  the  cathode.  These 
electrons  are  quickly  accelerated  by  the  strong  negative  electric  field 
through  the  region  known  as  Aston’s  dark  space  (see  Fig.  2).  Incoming 
positive  ions  colliding  with  the  slowest  of  the  exiting  electrons  emit 
light  as  they  r 3 combine  in  the  region  referred  to  as  the  cathode  glow. 
These  electrons  may  also  produce  excited  atoms,  which  emit  light  of  dis¬ 
crete  frequencies  as  they  return  to  their  ground  states.  The  remaining 
higher-energy  electrons  are  quickly  accelerated  through  the  cathode  dark 
space,  leaving  a  net  positive  space  charge  in  this  region.  (No  light  is 
emitted  from  this  region  because,  as  mentioned  previously,  the  recombina¬ 
tion  coefficient  is  very  small  at  high  relative  velocities.  The  net 
positive  space  charge  produces  a  strong  electric  field  in  these  three 
regions,  which  are  collectively  known  as  the  cathode  sheath  (Ref.  16:83). 
Typically,  most  of  the  voltage  drop  across  the  discharge  occurs  in  the 
cathode  sheath  (of  thickness  dc)  and  is  referred  to  as  the  cathode  fall, 

V 

Negative  Glow.  The  large  negative  electric  field  in  the  cathode 
sheath  causes  large  numbers  of  electrons  to  be  created  through  Townsend 
ionization.  However,  these  electrons  are  rapidly  accelerated  out  of  the 


Figure  2.  Glow  Discharge  Parameter  Profile 
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sheath  region  and  into  the  negative  glow  region,  where  they  produce  a 
slightly  negative  net  space  charge.  With  a  much  smaller  field  in  this 
region,  the  electrons  quickly  lose  their  energy  through  ionization  and 
recombination.  The  resulting  recombination  and  de-excitation  is  respon¬ 
sible  for  producing  the  "negative"  glow.  As  the  swarm  of  electrons  pass 
through  this  region,  recombination  gradually  reduces  the  net  space 
charge  to  zero  as  the  electrons  enter  the  next  region,  the  positive 
column  (Ref.  16:85). 

Positive  column.  In  spite  of  its  name,  the  positive  column  has 
little  or  no  net  space  charge.  Although  the  region  is  neutral  overall, 
the  more  rapid  diffusion  of  electrons  to  the  ’/alls  of  the  discharge  tube 
causes  a  "column"  of  positive  charge  to  develop  along  the  axis  of  the 
tube.  To  maintain  this  overall  neutrality,  the  small  negative  field 
within  the  positive  column  is  of  sufficient  strength  for  Townsend  ion¬ 
ization  to  balance  electron  loss  by  recombination  and  diffusion  (Ref. 
16:86).  These  electrons  also  cause  electron  impact  excitation,  so  that 
de-excitation  and  recombination  combine  to  produce  a  luminous  positive 
column.  However,  this  luminous  region  does  not  perform  any  essential 
function  in  the  production  of  a  glow  discharge.  It  merely  serves  to 
connect  the  cathode  and  anode  sheaths  ends  of  the  discharge.  Thus,  as 
the  discharge  parameters  (gas  pressure,  tube  length,  etc.)  are  varied, 
the  positive  column  adjusts  its  length  to  compensate  for  any  required 
'Charges  in  the  sheath  thicknesses. 

Anode  sheath.  A  short  distance  from  the  anode,  the  positive  column 
gives  way  to  the  anode  sheath,  which  somewhat  resembles  a  compressed 
cathode  sheath.  (The  anode  sheath  thickness  is  often  less  than  IX  of 
the  cathode  sheath  thickness.)  (Ref.  11:4999)  Also,  in  the  anode  sheath. 
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the  space  charge  produced  is  negative  since  electrons  are  attracted  to 
this  electrode.  This  negative  space  charge  causes  the  negative  electric 
field  in  the  discharge  to  increase  in  magnitude  in  the  sheath.  Some 
electrons  accelerated  by  this  enhanced  field  gain  enough  energy  to 
excite  or  ionize  neutral  gas  particles.  The  resulting  de-excitation  and 
recombination  produce  another  luminous  region,  known  as  the  anode  glow, 
some  small  distance  from  the  anode  surface  (Ref.  16:87).  The  increase 
in  magnitude  of  the  anode  field  produces  a  voltage  drop,  VA,  across  the 
anode  sheath,  of  thickness  d^.  Since  the  field  is  negative  throughout  a 
glow  discharge,  the  voltage  drops  across  the  cathode  and  anode  sheaths 
and  the  positive  column  are  all  positive.  Thus,  the  total  discharge 
voltage,  VT>  if  the  sum  of  these  voltage  drops  (See  Fig.  2). 

High  Pressure  Discharge.  If  the  pressure  is  increased  in  a  normal 

glow  discharge  (in  the  same  external  electrical  circuit) ,  the  thickness 

(d  )  of  the  cathode  sheath  will  decrease  and  the  current  density  in  the 
c 

sheath  will  increase.  In  addition,  the  field  in  the  positive  column 
will  be  greater:  the  increase  in  recombination  losses  at  higher  pres¬ 
sures  must  be  compensated  for  by  an  Increase  in  electron  impact  ioniza¬ 
tion  and  thus  an  increase  in  the  electric  field.  However,  at  higher 
pressures  (atmospheric  or  greater) ,  a  normal  glow  discharge  may  not  exist 
as  a  stable  discharge  (Ref.  16:88):  the  drop  in  voltage  upon  entering 
the  normal  glow  region  (see  Fig.  1)  may  no  longer  exist.  As  a  result, 
a  dark  discharge  may  transition  directly  to  an  arc  discharge  at  high  gas 
pressures . 

Arc  Discharge.  Although  an  arc  discharge  may  be  divided  into  regions 
similar  to  those  in  a  glow  discharge,  an  arc  discharge  is  able  to  carry 
a  much  higher  current,  at  a  lower  voltage,  than  a  glow  discharge.  (The 
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total  voltage  drop  across  an  arc  discharge  is  roughly  an  order  of  magni¬ 
tude  less  than  the  drop  across  a  glow  discharge.)  The  higher  current 
that  an  arc  carries  is  also  carried  at  a  higher  current  density.  This 
high  current  density  causes  the  discharge  to  be  intensely  luminous,  which 
prevents  the  observation  of  variations  in  the  sheath  regions.  The  sheath 
regions  are  compressed  and  the  voltage  drops  across  them  are  decreased. 
However,  the  electric  field  in  what  may  still  be  called  a  positive 
column  is  greater  than  for  a  glow  discharge. 

Experimentally-measured  I-V  Characteristics.  The  longitudinal  varia¬ 
tion  of  the  electric  field,  voltage,  and  number  and  current  densities  in 
a  gas  discharge  can  usually  be  obtained  experimentally,  although  with  some 
difficulty.  Hrwever,  it  is  usually  the  total  behavior  of  a  Discharge 
which  is  of  interest  in  practical  applications.  I-V  characteristics  of 
the  variation  c  f  the  discharge  (circuit)  current  with  total  discharge 
voltage  drop  are  experimentally  obtained  using  a  circuit  as  shown  in 
Figure  3.,  The  I-V  characteristic  is  "traced  out"  by  varying  the  supply 
voltage,  VQ,  for  a  particular  resistor  value,  R.  Thus,  given  Vq  and  R, 
and  measuring  the  external  circuit  current,  I,  the  discharge  voltage, 

Vj,(I)f  is  obtained.  The  relationship  between  these  circuit  parameters  is 
represented  by  what  is  termed  a  "load  line."  Several  load  lines  are 
superimposed  on  a  general  discharge  characteristic  in  Fig.  3.  As  VQ  is 
varied,  the  load  line  "scans"  the  characteristic.  If  the  characteristic 
has  a  complicated  behavior,  some  ambiguity  may  occur  if  the  resistance 
in  the  external  circuit  is  too  small,  as  shown  by  load  line  3  in  Fig.  3. 


Thus,  particularly  large-valued  resistors  must  be  used  to  trace  out 
reversals  in  I-V  characteristics,  as  in  the  Townsend  discharge  region. 


Experimental  I-V  Measurements  in  Gas  Discharges 


III.  Numerical  Method 


Externally-sustained  gas  discharges  may  contain  some  features  of 
each  of  the  types  of  discharges  just  described.  This  chapter  explains 
how  an  externally  sustained  gas  discharge  between  parallel  electrodes 
was  modeled  using  the  continuity  equations  for  electrons  and  ions  and 
Poisson's  equation.  With  the  boundary  conditions  imposed  by  the  elec¬ 
trodes,  this  set  of  equations  was  solved  numerically  to  obtain  the 
steady-state  spatial  distributions  of  the  electric  field  and  electron 
and  ion  number  densities  and  current  densities.  Three  methods  of  solu¬ 
tion,  subsequently  referred  to  as  Cases  Z,  II,  and  III,  were  used.  These 
methods,  each  of  which  corresponds  to  a  particular  type  of  discharge, 
were  used  by  Lowke  and  Davies  in  their  analysis  (Ref.  11:4999-5000). 
However,  before  these  individual  methods  are  described,  the  physical 
setup  of  the  components  of  the  discharge  tube,  the  sign  conventions 
used,  the  differential  equations  describing  the  discharge’s  behavior,  and 
the  boundary  conditions  at  the  electrodes  are  presented.  Additional 
boundary  conditions  are  discussed  separately  for  each  case. 

Setup  and  Sign  Conventions 

The  components  of  the  model  discharge  tube  are  shown  in  Fig.  4. 

The  problem  was  restricted  to  one  dimension  between  plane,  parallel 
electrodes  placed  a  distance  d  apart:  the  cathode  was  at  x  «  0  and  the 
anode  was  at  x  *  d.  It  was  also  assumed  that  an  external  ionization 
source  created  electron-ion  pairs  uniformly  throughout  the  electrode 
gap  at  a  constant  rate  S(l/cm  /sec).  (In  some  electron  beam-sustained 
discharges,  this  source  is  produced  by  directing  a  high  energy  electron 
beam  through  a  foil  at  the  cathode  end  of  the  discharge  tube.  To  allow 
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the  electron  beam  to  penetrate  the  entire  interelectrode  gap,  the  cathode 
le  In  the  form  of  a  conducting  grid.) 

The  positive  directions  of  the  electric  field,  current  densities, 
end  spatial  coordinates  are  defined  to  be  from  cathode  to  anode.  Thus, 
when  subjected  to  a  negative  electric  field,  an  electron,  with  a  negative 
charge,  will  have  a  positive  drift  velocity  and  move  toward  the  anode. 
Conversely,  a  positive  ion  will  have  a  negative  drift  velocity  and  move 
toward  the  cathode.  Thus,  when  an  electron-ion  pair  is  created  by  Town¬ 
send  ionization  or  an  external  ionization  source,  the  charged  particles 
migrate  to  opposite  electrodes.  Although  the  electron  drifts  much  faster 
than  the  ion,  and  would  thus  seem  to  carry  more  current,  it  exists  as  a 
charge  carrier  for  a  time  inversely  proportional  to  its  velocity.  Thus, 
the  total  contributions  to  the  current  of  the  electron  and  ion  are  exactly 
the  same.  Consequently,  the  current  carried  by  an  electron-ion  pair 
created  in  the  center  of  the  interelectrode  gap  is  equivalent  to  that 
carried  by  an  electron  created  at  the  cathode  or  a  positive  ion  created 
at  the  anode  (Ref.  17:154). 

However,  the  current  carried  at  a  particular  longitudinal  point  in 
the  Interelectrode  gap  may  be  shared  unequally  by  electrons  and  ions  if 
their  number  densities  do  not  exactly  compensate  for  differences  in  their 
drift  velocities.  At  all  longitudinal  points,  however,  the  total  current 
density,  which  is  the  sum  of  the  magnitudes  of  the  electron  and  ion  cur¬ 
rent  densities,  must  be  constant.  Thus,  in  terms  of  particle  fluxes, 
and  recalling  the  sign  convention  adopted,  the  total  particle  flux  at  a 
point  x  is: 

J/e  -  je(x)  -  j+(x)  (18) 
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Figure  4.  Setup  and  Sign  Conventions 


Equations 

The  continuity  equations  for  electrons  and  positive  ions  and  Pois¬ 
son's  equation  were  solved  simultaneously  to  obtain  the  spatial  variation, 
in  the  steady-state,  of  the  electric  field,  the  electron  and  positive  ion 
number  densities,  and  the  electron  and  positive  ion  fluxes.  Not  all  of 
the  physical  processes  described  in  the  previous  chapter  were  included  in 
every  calculation.  Initially,  only  those  processes  considered  by  Lowke 
and  Davies  (Ref.  11),  namely  electron  ionization  and  recombination,  were 
used.  Both  argon  and  methane  discharges  were  modeled  in  this  way.  Sub¬ 
sequent  calculations  included  the  effect  of  metastable  ionization  in 
argon  discharges  and  the  effect  of  secondary  mission  in  both  argon  and 
methane  discharges.  Including  terms  describing  metastable  ionization, 
the  continuity  equations  for  electrons  and  positive  ions  are  (secondary 
emission  affects  the  boundary  conditions  but  not  the  form  of  the  dif¬ 
ferential  equations): 


3ne  3je 

IT  “  S  +  “|wJne  -  Ven+  +  Rimne%  "  *3 T 


(19) 


3n+  3 j+ 

IT  "  8  +  a!wei ne  -  Rrnen+  +  Rimnenm  "  HT 


(20) 
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where  n^,  n+  -  number  density  of  electrons  and  positive  ions  (1/cm  ) 

3 

je,  J+  ■  flux  of  electrons  and  positive  ions  (1/cm  /sec) 
no  “  number  density  of  metastable  neutrals  (1/cm  ) 
x  ■  distance  from  cathode  (cm) 
t  ”  time  (sec) 

3 

S  ■  electron-ion  pair  production  rate  (1/cm  /sec) 
a  ■  Townsend's  first  ionization  coefficient  (1/cm) 
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3 

■  metastable  ionization  rate  coefficient  (cm  /sec) 

Rr  ■  electron-ion  recombination  rate  coefficient  (cm  /sec) 

We  «  electron  drift  velocity  (cm/sec) 

\ 

Equation  19  states  that  the  temporal  rate  of  charge  of  the  number  den¬ 
sity  of  electrons  at  a  particular  point  in  the  discharge  is  equal  to  the 
rate  at  which  electrons  are  produced  by  the  source  function,  Townsend 
ionization,  and  metastable  ionization,  minus  the  rate  at  which  electrons 
are  lost  by  recombination  and  the  gradient  in  electron  flux.  Similarly, 
equation  20  states  that  the  temporal  rate  of  charge  of  the  number  density 
of  positive  ions  is  equal  to  the  rate  at  which  positive  ions  are  produced 
by  the  source  function,  Townsend  ionization,  and  metastable  ionization, 
minus  the  rate  at  which  they  are  lost  by  recombination  and  the  gradient 
in  positive  ion  flux.  Thus,  the  continuity  equations  for  electrons  and 
positive  ions  are  identical.  In  the  steady-state,  the  temporal  deriva¬ 
tives  vanish,  and  the  result  is: 


3x  3x 


S  +  a|We[ne  -  Rrnen+  +  RimV 


m 


(21) 


In  equation  21,  an  expression  for  the  flux  of  electrons  at  a  parti¬ 
cular  point  in  the  discharge,  j£,  may  be  obtained  from  the  fluid  equation 
of  motion  (equation  1)  for  electrons.  Assuming,  as  before,  one  dimension, 
steady-state,  and  negligible  pressure  gradients,  and  using  the  ideal  gas 

law: 

-en  E  -  KT(dn  /dx)  +  m  n  VW  (22) 

6  6  6  6  6 


Dividing  by  m^v: 


“^e  KT  **ne 

mgV  11  vmev'  v  dx  ' 


n  W 

e  e 


(23) 
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Rewriting  this  equation  in  terms  of  the  conductivity,  0, 


2 

-ej  -e  n 
Je  _ e 

E  a  v 
e 


and  the  diffusion  coefficient  (equation  5), 
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0  e 

-  —E  ■  -D  ■—  +  n  W 
e  e3x  e  e 


3n 

+je  -  4-  neWe  -  DeTr 


(24) 

\ 


>.(25) 


(Ref.  11:4992)  (26) 


Equation  26  states  that  the  flux  of  electrons  at  a  particular  longitudinal 
point  is  due  to  the  electron  drift  velocity  (which  is  positive  for  nega¬ 
tive  electric  fields)  and  diffusion,  which  opposes  the  gradient  in  the 
electron  numbei  density.  Because  the  mobility  of  ions  is  usually  several 
orders  of  magnitude  less  than  the  mobility  of  electrons,  and  the  diffusion 
coefficient  is  directly  proportional  to  mobility,  the  diffusion  of  posi¬ 
tive  ions  was  neglected  in  this  analysis  (Ref.  16:34).  Thus,  the  flux 
of  positive  ions  was  assumed  to  be  given  by 

j+  “  n+W+  -  n+K+E  (27) 

(The  experimentally  determined  mobility,  k+,  of  A^+  ions  was  also  used 
for  the  mobility  of  positive  ions  in  methane,  although  their  relative 
masses  Indicate  that  the  mobility  of  positive  methane  ions  should  be 
approximately  five  times  greater  than  that  of  argon  ions.) 

The  potential,  V,  in  the  region  between  the  electrodes  is  given  by 
Poisson's  equation: 


(28) 
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where  the  independent  variables  are  j  ,  ne,  n+,  and  E.  The  potential, 

V,  is  obtained  directly  from  E  by  using  equation  29.  Also,  j+  is 
obtained  directly  from  n+  and  E  by  using  equation  27. 

In  using  this  set  of  equations,  the  value  for  the  external  ionization 
source  strength  had  to  be  stipulated  and  the  values  for  the  reaction 
rates,  ionization  coefficients  and  the  electron  diffusion  coefficient 
and  drift  velocity  had  to  be  known  as  functions  of  the  electric  field. 

The  values  used  for  some  of  these  terms  were  obtained  by  experimental 
aeasurement  of  the  quantities  in  a  uniform  field.  The  values  used  for 
other  terms  were  obtained  from  solutions  of  the  collisional  Boltzmann 
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equation,  in  which  electron  velocity  distributions  for  constant  fields 
and  electron  number  densities  were  used.  The  validity  of  using  the 
values  obtained  by  these  methods  was  thoroughly  discussed  by  Lowke  and 
Davies  in  their  report'.  They  also  emphasized  that  the  alternative  to  this 
equilibrium  approach,  which  would  probably  require  a  Monte  Carlo  techni¬ 
que,  would  be  considerably  more  complex  without  significantly  affecting 
the  major  features  of  the  solutions  obtained  (Ref.  11:4993). 

Methods  of  Solution 

Three  methods  were  used  to  solve  the  system  of  equations  describing 
an  externally-s  istained  discharge.  These  methods  were  used  by  Lowke  and 
Davies  in  their  investigation,  which  covered  a  very  wide  range  of  ion¬ 
ization  source  strengths.  In  the  first  two  Methods,  which  were  appli¬ 
cable  in  low  current  and  source  strength  sitaations,  Lowke  and  Davies 
found  that  no  plasma  region  was  produced.  However,  a  plasma  region  was 
produced  at  higher  source  strengths  and  currents,  which  required  the 
third  method  of  solution.  Although  all  of  Lowke  and  Davies'  methods 
were  successfully  duplicated,  the  high  currents  and  source  strengths  of 
Interest  in  switching  applications  made  their  third  method  much  more 
useful  than  the  others  in  this  investigation. 

In  all  three  methods,  integrations  were  begun  at  a  particular  point 
(or  points)  in  the  discharge  at  which  all  of  the  independent  variables 
were  either  known  or  iteratively  assigned  a  value.  The  position  of 
these  starting  points  depended  upon  the  case  being  considered:  in  Case 
It  the  anode  was  the  starting  point;  in  Case  II,  the  point  at  which  the 
field  vanished  was  the  starting  point;  and  in  Case  III,  two  starting 
points,  at  opposite  ends  of  the  positive  column,  were  used.  With  these 
starting  values,  and  with  all  of  the  transport  coefficients  and  material 
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functions  known  as  a  function  of  the  field,  the  integration  procedure 
of  C.  W.  Gear  (Ref.  25)  was  used  to  obtain  values  for  the  independent 
variables  at  various  points  in  the  discharge.  Gear's  method  was  parti¬ 
cularly  useful  because  it  automatically  decreased  the  step  size  only  in 
those  regions  for  which  the  system  of  equations  was  stiff. 


Boundary  Conditions 

For  all  three  methods  of  solution,  the  boundary  conditions  at  the 
electrodes  were  deduced  from  physical  requirements  on  the  ratio  of  elec¬ 
tron  to  ion  fluxes  at  the  electrodes.  The  resulting  boundary  conditions 
are  the  same  as  those  used  by  Lowke  and  Davie  3,  but  their  rationale  for 
them  is  based  upon  the  electrodes'  perfect  absorptivity  of  electrons 
(Ref.  11:4993).  The  boundary  conditions  used  in  this  analysis  will 
first  be  described  at  the  cathode.  It  will  be  seen  that  a  similar  devel¬ 
opment  applies  for  positive  electric  fields  at  the  anode,  while  some 
modification  is  required  if  the  field  is  negative  at  the  anode. 

At  the  cathode,  the  applied  voltage  produces  a  negative  field  which 
drives  positive  ions  into  the  electrode.  As  a  consequence,  electrons 
are  liberated  through  the  process  of  secondary  emission.  Assuming  that 
there  are  no  other  sources  of  electron  or  ion  currents,  the  secondary 
emission  yield  of  the  cathode  determines  the  ratio  of  the  electron  and 
ion  current  densities  at  that  point. 


Je  "  YJ+  “  (^r)  i 


(The  last  equality  is  obtained  from  the  requirement  that  the  total 
current  density  be  a  constant:  J/e  ■  j  -  j+.)  This  requirement  may 
also  be  expressed  in  terms  of  the  value  of  n£  at  the  cathode.  Since 
J+  ■  n+K+E,  »nd  using  Piosson's  equation  (33)  to  rewrite  n+  in  terms  of 
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a  and  —s — ,  equation  34  becomes 
e  ox 
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We  -  YW+ 


(35) 


Then,  in  the  case  of  no  secondary  emission,  y  *  0,  and  the  boundary  con¬ 
ditions  simplify  to: 

ie  -  0  (36) 


or 


n 

e 


(37) 


and  the  entire  current  is  carried  by  positive  ions.  In  most  of  the  dis¬ 
charge  calculations,  this  value  for  n£  was  over  eight  orders  of  magnitude 
less  than  the  maximum  value  for  n  in  the  discharge,  thus,  this  boundary 
condition  was  approximated  as 

ne  *  0  (38) 


which  agrees  with  Lowke  and  Davies'  development  (Ref.  11:4999).  This 
approximation  resulted  in  exaggerated  values  for  the  cathode  fall,  V£, 
although  they  were  usually  in  error  by  less  than  one  per  cent. 

At  the  anode,  if  the  electric-field  there  is  positive,  positive  ions 
are  driven  toward  the  electrode  and  it  should  behave  as  a  cathode.  Thus, 
the  number  density  of  electrons  at  the  anode,  for  no  secondary  emission, 
is  the  same  as  at  the  cathode: 


n 


e 


-£  _£  z  o 
fix 


(39) 


Although  it  is  also  required  that  J  -  0,  this  alternate  boundary 
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condition  could  not  be  satisfied  because  of  an  error  in  the  treatment  of 
electron  diffusion.  This  error  occurs  whenever  there  is  appreciable 
electron  diffusion  into  a  retarding  electric  field.  Although  these  elec¬ 
trons  should  lose  energy,  if  the  diffusion  is  into  a  region  with  a 
stronger  electric  field,  the  material  functions  used  assign  the  electrons 
a  larger  diffusion  coefficient.  This  problem  is  particularly  evident  in 
Case  II  discharges,  in  which  diffusion  has  a  significant  effect  11: 

4998).  Thus,  as  seen  in  Figure  7,  the  electron  flux  continues  to  increase 
in  the  anode  "sheath"  although  the  electron  number  density  is  decreasing 
and  the  retarding  field  is  also  increasing.  A  similar  problem  does  not 
occur  in  the  cathode  region  because  the  electric  field  there  decreases 
with  increasing  x.  The  anode  sheath  in  Case  III  discharges  suffers  from 
the  same  problem,  although  its  overall  effect  is  not  as  important  because 
the  anode  sheath  is  usually  much  smaller  than  the  cathode  sheath  in  this 
case.  Thus,  if  the  anode  electric  field  is  positive,  the  most  that  can 
be  required  of  the  particle  fluxes  at  the  anode  is 

je  -  j+  -  j/e  (40) 

If  the  anode  electric  field  is  negative,  electrons  are  driven  toward 
the  cathode  and  the  ratio  of  positive  ion  to  electron  fluxes  is  again  a 
constant 

J+  -  13. 

In  this  case,  however,  n  is  always  zero  because  the  probability  for  an 
electron  to  induce  the  emission  of  a  positive  ion  is  vanishingly  small. 
Thus,  j+  «  0,  and  the  entire  current  is  carried  by  electrons.  In  terms 
of  number  densities,  then,  the  number  density  of  positive  ions  at  the 
anode  is: 
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0 


(41) 


while  the  number  density  of  electrons  Is  (using  Poisson’s  equation): 


-e 


_ o  6E 

*e  e  fix 


(42) 


If  fiE/fix  at  the  anode  is  very  small,  as  in  a  Case  I  discharge,  this 
boundary  condition  becomes 


ne  s  0 


(43) 


These  electrode  boundary  conditions  are  summarized  in  Table  II.  Addi¬ 
tional  boundary  conditions  at  the  ends  of  a  Case  III  positive  column 
and  at  the  center  of  a  Case  II  discharge  are  discussed  as  they  are 
encountered  in  the  following  descriptions  of  the  three  methods  of  solu¬ 
tion. 

Case  I.  This  method  was  used  for  low,  saturated  current  densities 
and  low  ionization  source  strength  discharges.  In  this  situation,  the 
sheath  thickness  is  greater  than  the  inter-electrode  gap  and  there  is 
no  plasma  region  with  ng  -  n+.  In  fact,  ng  Is  generally  much  less  than 
n+  throughout  the  inter-electrode  gap.  An  example  of  this  situation  is 
shown  in  Fig.  5,  for  conditions  identical  to  those  used  by  Lowke  and 
Davies  (Ref.  11:4995).  At  this  low  applied  voltage,  the  field  in  the 
discharge  is  correspondingly  small  and  Townsend  ionization  is  minimal. 
Thus,  the  maximum  current  density  that  can  be  carried  at  these  low 
voltages  is  just  the  saturation  current  density,  *  edS  (Equation  17) 
However,  the  Case  I  method  is  only  designed  to  handle  discharges  at  the 
saturation  current  density.  At  voltages  lower  than  that  required  to 
just  sweep  all  electrons  and  ions  to  the  electrodes  before  recombining, 
the  Case  II  method  should  be  used.  This  dividing  point  can  also  be 
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described  as  Che  discharge  voltage  for  which  the  cathode  sheath  thick¬ 
ness,  d£,  equals  the  electrode  separation,  d.  Thus,  at  voltages  for 

which  Townsend  ionization  is  negligible.  Case  I  applies  if  d  >.  d  and 

c 

Case  II  applies  if  dc  <.  d.  The  regions  of  applicability  of  these  methods 
are  shown  in  Fig.  6,  which  is  taken  from  Lowke  and  Davies'  report. 


Fig.  6.  I-V  Characteristic  for  Case  I  and  Case  II  (S  *  8  x  10^*) 
(Ref.  11:4995) 


In  order  to  solve  the  system  of  equations  (30-33)  using  the  Case  I 

o 

method,  a  "shooting"  technique  is  used  because  three  boundary  conditions 
are  known  at  the  anode  and  one  is  known  at  the  cathode.  At  the  anode, 
as  discussed  previously,  the  improbability  of  positive  ion  emission 
requires  j+  ”  0,  which  in  turn  requires  that  n+  *  0  and  j  *  J/e.  Then, 
since  n+  ■  0,  Poisson's  equation  requires,  at  the  anode,  that  (equation 


“S  6E 
e  6x 


42): 


(44) 


Boundary  Conditions 


In  a  Case  I  discharge,  however,  6E/6x  at  the  anode  Is  essentially  zero, 
so  that  another  anode  boundary  condition  Is  that  ng  -  0.  Also  at  the 
anode,  the  voltage  Is  arbitrarily  assigned  a  value  of  zero  and  the  field 
is  Iteratively  assigned  values  In  order  to  satisfy  the  cathode  boundary 
condition  that  ng  *  0,  or  j  ■  yj+  if  secondary  emission  is  included. 

Thus,  with  the  total  discharge  current  density  (J)  the  only  Input  para¬ 
meter,  the  spatial  variation  of  nfi,  n+,  j+,  E,  and  V  is  obtained, 

with  the  total  voltage  drop  across  the  discharge,  V^,  calculated  by 
taking  the  difference  between  the  anode  and  cathode  voltages. 

Case  II.  At  voltages  below  the  value  for  which  d£  -  d,  the  current 
in  the  discharge  falls  below  the  saturated  value  and  the  spatial  variation 
of  the  independent  variables  takes  on  a  form  very  different  from  that  of  a 
Case  I  discharge.  (See  Fig.  7,  which  duplicates  Lowke  and  Davies'  solu¬ 
tion  for  a  Case  II  discharge.)  The  existence  of  central  maxima  in  n  and 
n+  is  a  result  of  space  charge  effects,  which  are  much  more  important 
than  in  Case  I:  while  the  electric ^field  in  a  Case  I  discharge  Is  farily 
uniform,  positive  space  charges  at  both  the  electrodes  in  a  Case  II  dis¬ 
charge  cause  the  field  to  be  negative  near  the  cathode  and  positive  near 
the  anode.  This  field  reversal  forces  electrons  to  the  center  of  the 
discharge  and  positive  ions  to  the  electrodes,  creating  maxima  in  n£  and 
n+  at  the  point  where  E  passes  through  zero.  However,  these  maxima  in 
the  number  densities  are  prevented  from  becoming  a  plasma,  with  ne~n+, 
because  of  the  "influence  of  diffusion  together  with  the  boundary  condi¬ 
tions  that  n  -  0  at  the  electrodes  (Ref.  11:4995).  Because  positive 
e 

ion  diffusion  is  neglected,  the  number  density  of  positive  ions  may  be 
expressed  as 
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Figure  7.  Case  II  Discharge  (Argon,  240  Torr,  S  ■  8  x  101X  1/cnr/sec 
J  -  .022  yA/cnT) 


Thus,  at  the  point  where  E  ■  0,  j+  must  also  be  zero  to  avoid  a  singu¬ 
larity  in  n+.  The  value  for  n+  at  this  point  is  obtained  by  using 
L'Hospital's  rule: 


dj.  dE 

■v  ■  (-dr>/(k+f> 


Making  use  of  the  continuity  equation  for  positive  ions  and  Poisson 
equation,  the  expression  obtained  for  n+  is: 


(46) 

’s 


(k+e/c  -  y)ne  +  Ky  -  k+e/e)2ne2  +  Ak^eS/e]1^2 
__ 


(47) 


Thus,  with  a  trial  value  for  n  at  this  point,  n,  is  also  determined. 

In  addition,  with  j+  *  0,  j  must  equal  J/e.  Integrating  toward  the 
cathode  from  this  central  point,  the  position  of  the  cathode  is  defined 
to  be  at  that  point  for  which  ne  »  0  or  J  *  yj+,  depending  upon  whether 
secondary  emission  is  allowed  or  not.  Integrating  toward  the  anode,  the 
anode's  position  is  defined  to  be  where  ng  -  0.  The  value  of  ne  at  the 
central  point  is  Iteratively  varied  so  that  the  calculated  distance 
between  the  anode  and  cathode  is  exactly  the  physical  separation  of  the 
electrodes  (Ref.  11:4999). 

Case  III.  This  method,  which  is  applicable  at  higher  current  den¬ 
sities  and  much  higher  source  strengths  than  the  previous  two  methods, 
produces  two  types  of  solutions  determined  by  the  magnitude  of  the  dis¬ 
charge  current  (see  Figs.  8  and  9).  In  both  types,  a  plasma  region  is 
created  with  n  -  n,  and  an  electric  field  much  less  than  that  in  the 

€  T 

2 

cathode  sheath.  At  current  densities  below  lmA/cm  (for  S  -  3.6^),  a 
constant  number  density,  positive-ion  space  charge  is  formed  in  both  the 
cathode  and  anode  sheaths  as  a  result  of  the  differing  drift  velocities 
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Figure  8.  Case  III  Discharge  (Argon,  240  Torr,  S  -  3.6  x  1010  l/cnr/sec 
J  -  .254  mA/cm2) 


of  electrons  and  positive  Ions  and  the  relatively  weak  electric  field  in 

the  discharge  (Ref.  11:4996).  However,  at  current  densities  above  -25  mA/ 
2 

cm  ,  the  increased  electric  field  In  the  cathode  sheath  causes  signifi¬ 
cant  Townsend  ionization  to  occur,  producing  a  positive-ion  number  den¬ 
sity  maximum  which  is  greater  than  the  number  density  in  the  plasma 

region.  The  maximum  in  n+  exists  because  Townsend  ionization,  which  is 

l  /*  '*■ 

represented  by  the  term  a[  W  |n  ,  is  small  near  the  cathode,  where  n  *  0, 
and  also  near  the  plasma  region,  where  the  electric  field  strength  (and 
thus  the  electron  drift  velocity)  is  small  (Ref.  11:4997).  At  these  high 
current  densities,  the  electric  field  is  also  much  stronger  in  the  posi¬ 
tive  column.  This  causes  a  larger  portion  of  the  total  voltage  drop 
across  the  discharge  to  occur  in  the  positive  column  (6V/6x  =  -E) .  In 
addition,  these  high  fields  are  able  to  draw  positive  ions  away  from  the 
anode  sheath,  leaving  an  excess  of  electrons  there.  Thus,  the  anode 
sheath  has  a  negative  space  charge,  which  drives  the  field  further  nega¬ 
tive  in  the  sheath  and  reinforces  these  tendencies. 

For  both  the  low  and  high  current  Case  III  solutions,  integration 
Is  begun  In  the  positive  column,  where  the  values  of  all  of  the  independ¬ 
ent  variables  can  be  determined.  In  the  positive  column,  which  is  quasi¬ 
neutral,  the  number  densities  of  electrons  and  ions  are  assumed  to  be 
equal.  Also,  at  the  relatively  low  electric  fields  in  the  plasma,  Town¬ 
send  ionization  is  minimal  and  the  external  source's  production  of  elec¬ 
trons  and  positive  ions  is  balanced  solely  by  recombination  losses.  Thus, 
in  the  plasma: 

ne  -  n+  -  (S/Ry)1'2  (48) 

Also,  because  of  their  higher  mobility,  the  electrons  in  the  plasma  are 
assumed  to  carry  the  total  discharge  current  density,  J, 
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Figure  9.  Case  III  Discharge  (Argon,  240  Torr 


(49) 


at  a  constant  velocity,  W^, 


».  -  v«,  on 

since  the  electron  number  density  is  uniform  throughout  the  plasma  region. 
Using  tabulated  values  for  the  drift  velocity  versus  E/N  (see  Appendix  A), 
the  positive  column  field  is  also  determined.  Then,  since 


J+  -  n+k+E  (51) 

the  positive  ion  current  density  can  be  obtai ted  and  the  electron  current 
density  recalculated: 

ie  -  J/e  +  j+  (52) 

Thus,  using  equations  50  and  51  in  52,  an  expression  is  obtained  for 

-  J/e  +  nek+FLD(je/ne)  (53) 

where  FLD  is  a  tabulated  function  which  gives  the  electric  field  for  a 
particular  electron  drift  velocity.  This  expression  is  iteratively 
evaluated  until  the  values  converge.  Thus,  values  for  all  of  the  inde¬ 
pendent  variables  (ng,  n+,  jg,  and  E)  are  determined  at  the  cathode  end 
of  the  plasma  region  (see  Table  II).  With  these  starting  values,  the 
cathode  sheath  region  is  integrated  and  the  cathode  is  defined  to  be  at 

that  point  where  n  ■  0,  or  j  ■*  yj .  t  if  secondary  emission  is  considered. 
6  6  * 

The  starting  values  used  for  the  anode  sheath  integration  are  deter¬ 
mined  by  the  discharge  current  density.  At  low  current  densities 
2  16 

(J  <.  10  taA/ca  for  S  *  3.6  x  10  ),  the  field  in  the  positive  column  is 
also  low.  Thus,  the  positive  ions  produced  by  the  external  source  are 
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sot  removed  from  the  anode  sheath,  and  a  positive  space  charge  is  formed. 
This  positive  space  charge  causes  the  electric  field  to  become  positive 
in  the  anode  sheath.  Thus,  the  anode  acts  as  a  cathode  and  the  anode 
boundary  conditions  are  similar  to  those  at  the  cathode.  To  start  this 
integration,  the  field  is  assigned  a  very  small  positive  value  of  .1  V/ 
cm.  Then,  with  the  positive  column  values  for  n and  n+,  j+  is  calculated 
as  described  previously.  Since  the  total  current  density  is  constant, 
j  is  also  determined.  With  these  starting  values,  the  anode  sheath  is 
integrated  and  the  anode  is  defined  to  be  where  n  ■  0,  or  j  *  yj, ,  if 
secondary  emission  at  the  anode  is  considered .  (In  this  analysis,  second¬ 
ary  emission  at  the  anode  was  not  considered,  so  the  boundary  condition 
used  in  this  case  was  ng  **  0.) 

2 

In  contrast,  at  high  current  densities  (J  >  500  mA/cm  for  S  • 

3.6  x  10*^),  the  anode  sheath  has  a  negative  space  charge  and  the  electric 
field  goes  to  larger  negative  values.  In  this  case,  the  starting  value 
for  the  electric  field  at  the  positive  column/anode  sheath  boundary  is 
Just  the  field  value  In  the  positive  column.  Using  the  positive  column 
values  for  ng  and  n+,  j+  and  j  are  also  calculated.  With  these  starting 
values,  the  anode  sheath  is  integrated  and  the  anode  is  defined  to  be 
where  n+  ■  0. 

At  Intermediate  currents,  two  anode  sheath  solutions  are  possible 
for  each  value  of  the  discharge  current.  For  one  solution,  corresponding 
to  the  low  current  density  case,  the  generally  negative  electric  field 
reverses  and  goes  positive  in  the  anode  sheath;  for  the  other  solution, 
corresponding  to  the  high  current  density  case,  the  field  goes  further 
negative  in  the  sheath.  Using  the  low  current  density  method,  as  the 
current  density  is  increased,  the  value  of  the  field  at  the  anode 
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decreases.  At  high  current  densities,  this  method  can  no  longer  produce 
a  solution.  Using  the  high  current  density  method,  as  the  current  den¬ 
sity  is  decreased,  a  solution  is  obtained  only  if  the  starting  value  for 

\ 

the  field  is  perturbed  positively  from  the  positive  column  value.  When 
the  required  starting  value  approaches  zero,  the  set  of  equations  becomes 
unsolvable  (See  Figure  10) . 

IP  v 

To  calculate  the  total  voltage  drop  across  a  Case  III  discharge, 
these  different  methods  of  solution  had  to  be  considered.  Using  the 
high  current  method  for  high  currents,  the  field  is  assumed  to  be  con¬ 
stant  throughout  the  positive  column.  Thus,  since  <$V/6x  *  -E,  the  volt¬ 
age  drop  across  the  positive  column,  V  ,  is 

VP  -  -Ecs  ’  dP  <54> 

where  dp  is  the  length  of  the  positive  column  and  is  the  value  of  the 
field  at  the  cathode  sheath/positive  column  interface.  However,  the 
fields  at  opposite  ends  of  the  positive  column  are  different  when  using 
the  low  current  method  and  for  low  currents  when  using  the  high  current 
method.  In  these  instances,  it  is  assumed  that  the  field  varies  linearly 
between  these  values.  Thus,  the  voltage  drop  across  the  positive  column 
is  reduced,  being  given  by 

Tp  ■  -ECS  •  dP  -  <'as  -  ECS>  •  V2  (55) 

where  E  is  the  value  of  the  field  at  the  anode  sheath/positive  column 

as 

Interface  and  Eafl  >  E^g. 

At  Intermediate  currents,  for  which  both  the  high  and  low  current 
density  methods  produce  a  solution,  the  total  voltage  drop  across  the 
discharge  can  have  two  different  values.  As  a  result,  there  must  be  some 
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criteria  to  judge  which  is  physically  correct.  In  this  investigation, 
the  voltage  drop  given  by  the  high  current  density  method  was  used  when¬ 
ever  both  the  high  and  low  current  density  methods  produced  a  solution. 
As  a  result,  the  voltages  chosen  were  higher  than  their  alternatives, 
although  the  difference  was  usually  less  than  ten  percent. 
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IV.  Results 


The  validity  of  the  computer  codes  written  to  model  externally- 
ionized  gas  discharges  was  ascertained  by  comparing  their  results  to 
the  results  of  Lowke  and  Davies’  calculations.  For  each  of  the  codes, 
the  spatial  variation  of  the  number  and  current  densities,  the  electric 
field,  and  the  voltage  (the  discharge's  "profile")  was  calculated  for 
the  same  parameters  specified  by  Lowke  and  Davies.  The  profile  cal¬ 
culated  for  a  Case  I  discharge  is  shown  in  Fig.  5.  Although  the  profile 
shown  nearly  duplicates  the  profile  calculated  by  Lowke  and  Davies  (Ref. 
11:4995),  there  is  one  major  discrepancy.  In  order  to  obtain  this  pro¬ 
file,  Lowke  and  Davies'  cathode  boundary  condition,  ng  =  0,  was  relaxed 

A 

and  assumed  to  be  satisfied  if  n  <  1  x  10  .  Without  this  modification, 

e 

the  total  voltage  drop  across  the  discharge  was  about  an  order  of  magni¬ 
tude  too  large.  The  proper  boundary  condition  at  the  cathode  was  given 
previously  in  equation  37: 


e 


Although  this  more  exact  cathode  boundary  condition  was  not  incorporated 
in  any  of  the  computer  codes  written,  the  effect  of  its  neglect  was 
decreased  in  both  Case  II  and  Case  III  discharges.  Their  larger  elec¬ 
tron  number  densities  and  gradients  allowed  the  approximate  boundary 
condition  of  ng  s  0  to  be  quite  adequate,  especially  for  Case  III  dis¬ 
charges.  However,  this  particular  problem  was  completely  avoided  in 
most  of  the  calculations  because  the  secondary  emission  boundary  condition 
(je  *  YJ+)  was  usually  used  at  the  cathode. 

The  profile  calculated  for  a  Case  II  discharge  also  nearly  dupli¬ 
cated  Lowke  and  Davies'  results  (Ref.  11:4995)  (See  Fig.  7.)  In 
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addition  to  the  error  In  the  cathode  boundary  condition,  the  peaks  In 
n^  and  n+  also  exhibited  a  small,  sharp,  non-physical  peak  centered  on 
the  starting  point  for  the  integration.  Although  this  "cosmetic"  defect 
could  be  overcome  simply  by  taking  a  larger  step  at  the  start  of  each 
integration,  this  defect  may  be  the  result  of  starting  conditions  at  the 
peak  which  are  not  completely  consistent. 

The  validity  of  the  Case  III  code  was  checked  more  thoroughly  than 

the  Case  I  and  Case  II  codes.  Both  individual  discharge  profiles  and  I-V 

characteristics  generated  from  a  large  number  of  profiles  were  checked 

against  Lowke  and  Davies'  calculated  and  theoretical  results.  The  pro- 

2  2 

files  obtained  for  current  densities  of  .254  riA/cm  and  48.  mA/cm  ,  shown 
in  Figs.  8  and  9,  were  nearly  identical  to  Lowke  and  Davies'  calculated 
profiles  (Ref.  11:4996-4997).  Although  the  cathode  discrepancies  seen 
In  the  Case  I  calculation  were  also  present  in  these  profiles,  their 
effects  were  negligible  in  this  case.  A  profile  was  also  obtained  in 
which  electron  diffusion  was  effectively  neglected  by  setting  its  coef¬ 
ficient  to  one  thousandth  of  the  electron  mobility  (see  Fig.  11).  The 
calculated  cathode  sheath  thickness  in  this  case  was  0.0441  cm,  which 
compared  very  favorably  with  the  thickness  of  0.044  cm  predicted  by 
Lowke  and  Davies'  theoretical  development  (Ref.  11:4996). 

The  effect  of  diffusion  can  be  seen  by  comparing  Figs.  8  and  11. 

Ions  produced  in  the  plasma  region  by  the  source  function  are  lost  to 
recombination.  Since  these  ions  never  get  to  the  electrodes,  they  can¬ 
not  contribute  to  the  total  current  in  the  discharge.  Diffusion  shortens 
the  plasma  region  and  reduces  the  number  of  charge  carriers  which  must  be 
produced  by  Townsend  ionization  in  the  sheath.  Thus,  when  diffusion  is 
significant,  a  smaller  sheath  field  (and  voltage  drop)  are  sufficient 
to  supply  the  charge  carriers  needed  for  a  given  current  to  flow. 
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Figure  12.  Case  III  Discharge  (Argon,  240  Torr,  S  -  3.6  x  10  ,  J  -  .254  mA/ 


Townsend  ionization  also  decreases  the  total  voltage,  particularly 
at  higher  sheath  field  strengths  (which  occur  for  higher  current  densi¬ 
ties).  This  effect  can  be  seen  by  comparing  Figs.  8  and  12.  Figure  8 
shows  the  discharge  profile  for  a  Case  III  discharge  with  both  diffusion 
and  Townsend  ionization,  while  Figure  12  shows  the  profile  obtained  when 
Townsend  ionization  is  neglected.  The  Increase  in  the  total  voltage  when 

rp  +* 

Townsend  ionization  is  neglected  (Fig.  12)  is  most  easily  observed  at  the 
anode,  where  it  is  apparent  that  the  cathode  fall  in  this  case  is  larger 
than  in  Fig.  8. 

These  results  are  summarized  for  a  number  of  discharge  currents  in 
the  I-V  characteristics  shown  in  Fig.  13.  As  shown  there,  with  the  addi¬ 
tion  of  Townsend  ionization,  diffusion,  or  both,  the  total  voltage  drop 
Is  decreased  at  each  current.  Note  that  the  decrease  due  to  Townsend 
ionization  is  most  pronounced  at  higher  current  densities,  while  the 
decrease  due  to  diffusion  is  more  pronounced  at  lower  current  densities 
(and  therefore  lower  cathode  fields).  Although  these  results  agree  very 
closely  with  those  of  Lowke  and  Davies  (Ref.  11:4996),  both  sets  of 
calculations  underestimate  the  voltage  drop  observed  in  Leffert's  exper¬ 
iments  (Ref.  13).  Lowke  and  Davies  offered  two  explanations  for  this 
discrepancy.  First,  the  set  of  equations  does  not  properly  account  for 
the  diffusion  of  electrons  into  a  retarding  electric  field,  as  in  the 
cathode  sheath.  Although  electrons  should  lose  energy  when  moving 
contrary  to  such  a  field,  the  increase  in  the  diffusion  coefficient  with 
E/N  produces  the  opposite  effect  in  electrons  diffusing  from  the  low  field 
positive  column  through  the  high  field  cathode  sheath.  Thus,  the  cal¬ 
culated  effect  of  diffusion  is  overestimated  and  the  cathode  fall  is 
underestimated  (Ref.  11:4998).  Lowke  and  Davies  also  ventured  a  second 


Fig.  13.  Calculated  and  Analytic  I-V  Characteristics  (Argon, 

240  Torr,  S  ■  3.6  x  10^  l/cm^/sec,  with  and  without 

a  and  D  ) 
e 


explanation  based  upon  nonuniformities  in  the  external  ionization  in 

Leffert's  gas  discharge.  (They  noted  that  Leffert's  source  function  at 

the  electrodes  was  about  half  of  its  value  at  the  center  of  the  inter- 

electrode  gap.)  (Ref.  11:4996)  To  see  if  this  second  idea  could  be 

used  to  explain  the  discrepancy  between  the  calculated  and  experimental 

I-V  curves,  the  effect  of  varying  the  source  function  was  investigated. 

Discharge  profiles  were  calculated  for  an  argon  discharge  at  240 
2 

Torr  and  J  «  15  mA/cm  .  Townsend  ionization  and  electron  diffusion  were 
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included.  The  results  obtained  for  three  different  source  strengths 
are  shown  in  Figs.  14,  15,  and  16.  At  higher  source  strengths,  the 
cathode  sheath  is  compressed,  the  cathode  fall  (and  total  voltage  drop) 
is  decreased,  and  the  electric  field  at  the  cathode  is  decreased.  Since 
the  number  of  electrons  in  the  cathode  sheath  is  larger  at  higher  S,  the 
cathode  sheath  field  need  not  be  as  large  for  the  discharge  to  carry  the 
same  current.  With  a  decrease  in  the  field,  the  voltage  drop  also 
decreases  since  =  -E.  This  effect  is  predicted  by  Lowke  and  Davies’ 
analytic  solution,  although  it  did  not  include  the  effects  of  ionization 
or  diffusion.  The  expression  they  obtained  for  the  total  current  density 
in  the  discharge  is 

J  -  <4eok+e3S3)1/4  VT1/2  (56) 

Thus,  the  total  current  density  is  proportional  to  S3^4  and  and, 

-3/2 

at  a  given  current,  VT  is  proportional  to  S  .  However,  the  voltage 

drop  predicted  by  this  relationship  is  much  larger  than  that  calculated 

2 

for  a  current  density  of  15  mA/cra  .  At  this  current,  ionization  in  the 
cathode  sheath  is  substantial,  as  shown  by  the  peak  in  n+.  At  lower  cur¬ 
rents,  ionization  is  less  important,  and  Lowke  and  Davies’  relationship 

is  more  closely  obeyed.  I-V  characteristics  for  three  different  source 

2 

strengths  are  shown  in  Fig.  17.  At  a  current  density  of  .5  mA/cm  , 
increasing  S  by  a  factor  of  four  decreased  V^,  to  about  .13  times  its 
former  value,  which  is  close  to  the  value  of  .125  predicted  by  Lowke  and 
Davies’  relationship.  At  lower  current  densities,  the  decrease  in  V^.  is 
not  as  large,  so  that  the  relationship  between  and  S  is  given  approxi¬ 
mately  by  V^,  -  1/S.  Using  this  new  relationship,  the  discrepancy  between 
Leffert’s  experiment  and  the  present  calculations  can  be  partially 
explained.  With  a  50%  decrease  in  S  at  the  electrodes,  or  an  average 
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Figure  IS.  Case  III  Discharge  (Argon,  240  Torr,  S  ■  3.6  x  101  ,  J  ■  15.  mA/cm^ 


Figure  16.  Case  III  Discharge  (Argon,  240  Torr,  S  «  3.6  x  10  ,  J  -  15.  mA/cm' 


of  about  S8Z  within  the  cathode  sheath  (using  a  sheath  thickness  of  .044 
cm),  and  assuming  V^,  ~  1/S,  the  predicted  increase  in  voltage  is  approx¬ 
imately  1.72.  Extrapolating  the  experimental  I-V  curve  in  Fig.  13  to 
2  x  10  1/cm  / sec  (.32  mA/cm  ),  the  expected  experimental  voltage  is 
approximately  1.67  times  the  calculated  value,  which  agrees  favorably 
with  the  predicted  increase. 

However,  the  reason  for  this  departure  from  Lowke  and  Davies  *  theory 
is  not  completely  understood.  The  source  strength  does  affect  the  pro¬ 
cesses  of  ionization  and  diffusion.  These  processes  could  have  opposite 
effects  upon  the  voltage  drop  in  the  cathode  sheath,  and  also,  therefore, 
upon  the  total  discharge  voltage.  At  various  current  densities,  one  of 
these  effects  could  be  dominant,  which  could  explain  the  departure  from 
the  relationship  between  and  S  predicted  by  Lowke  and  Davies,  who 
neglected  these  effects. 

At  increased  source  strengths,  the  effect  of  collisional  ionization 
is  decreased,  as  can  be  seen  in  Figs.  14,  15,  and  16.  The  peak  in  the 
positive-ion  number  density,  which  is  a  direct  result  of  ionization  in 
the  sheath,  changes  little  in  magnitude  for  different  source  strengths. 

As  S  is  Increased,  greater  numbers  of  electrons  diffuse  into  the  cathode 
sheath  and  more  electron-ion  pairs  are  produced  in  the  sheath,  so  that 
ionization  need  not  supply  as  many  electrons  and  ions  as  before  in  order 
to  carry  the  same  current.  Thus,  the  cathode  field  is  reduced  and  the 
maximum  in  n+  is  reached  nearer  to  the  cathode.  Adjacent  to  the  positive 
colimn,  in  the  region  known  as  the  negative  glow  in  glow  discharges, 
recombination  occurs  as  the  electrons  lose  energy.  At  higher  source 
strengths,  the  cathode  field  is  smaller  and  larger  numbers  of  electrons 
are  present  in  this  region.  Thus,  recombination  increases  and  the  member 
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Curr 


Figure  17.  Calculated  1-V  Characteristics  for  Three  Source  Strengths 
(Argon,  240  Torr) 
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density  of  positive  ions  is  reduced.  This  effect  is  more  pronounced  at 
higher  source  strengths,  as  can  be  seen  by  comparing  the  cathode  sheath 
minima  in  the  n+  profiles  in  Figures  14,  15,  and  16.  At  these  higher 
eource  strengths,  the  cathode  sheath  field  is  reduced,  Townsend  ioniza¬ 
tion  is  thereby  decreased,  and  a  greater  portion  of  the  discharge  current 
is  due  to  the  external  source.  Thus,  the  conductivity  of  the  discharge 
is  more  completely  controlled  by  the  external  source.  However,  the  large 
energies  invested  in  producing  these  high  source  strengths  reduces  the 
gain  in  switching  applications.  High-gain  switching  seeks  to  make  use 
of  the  weakest  source  able  to  sustain  a  high-current  discharge,  forcing 
the  cathode  shes th  to  generate  the  necessary  charge  carriers. 

One  method  of  increasing  the  effectiveness  of  the  cathode  sheath 
in  achieving  this  goal  is  by  coating  the  cathode  with  a  material  which  has 
a  high  secondary  emission  coefficient  in  the  gas  being  used.  With  an 
Increase  in  y,  the  electron  current  density  (and  thus  the  electron  number 
density)  at  the  cathode  is  increased.  The  effect  of  these  additional 
electrons  is  magnified  by  the  cascade  of  ionization  which  each  produces. 
Thus,  the  cathode  field  required  to  produce  the  additional  required  elec¬ 
trons  and  ions  is  reduced.  Consequently,  a  given  current  can  be  carried 
at  a  much  lower  total  discharge  voltage.  This  effect  can  be  seen  by  com¬ 
paring  Fig.  18,  with  no  secondary  emission,  to  Fig.  19,  which  includes  a 
2Z  secondary  emission  coefficient  at  the  cathode.  The  peak  in  n+,  which 
is  a  result  of  Townsend  ionization  in  the  sheath,  is  shifted  nearer  the 
cathode  as  y  is  increased.  As  a  result,  an  increase  in  y  causes  a 
decrease  in  the  cathode  fall  voltage.  (Note  the  shift  in  the  voltage 
profile  from  Fig.  18  to  Fig.  19.)  For  discharges  in  which  electrons  and 
ions  are  produced  solely  by  Townsend  ionization,  von  Engel  predicts  a 
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linear  relationship  between  the  cathode  fall  and  the  natural  logarithm 
of  y  (Ref.  10:1077).  As  shown  in  Fig.  20,  this  is  nearly  the  case  for 
low  source  strengths.  However,  at  higher  source  strengths,  the  effect 
of  Townsend  ionization-  is  decreased  and  small  values  of  y  do  not  appreci¬ 
ably  decrease  the  cathode  fall.  Thus,  larger  values  of  y  are  required 
for  ionization  effects  to  come  into  play.  At  the  source  strengths 
required  in  electron-beam  controlled  switches,  the  cathode  fall  voltage 
is  relatively  insensitive  to  the  secondary  emission  coefficient.  However, 
some  gas/electrode  material  combinations  have  unusually  high  secondary 
emission  coefficients  (e.g.  for  cesium  electrodes  in  argon,  y  *  .4  elec¬ 
trons  per  incident  metastable  argon  atom)  (Ref.  17:84).  In  these 
Instances,  the  decrease  in  cathode  fall  voltage  could  be  considerable. 

This  decrease  in  voltage  with  increasing  y  is  shown  by  the  I-V 

characteristics  in  Fig.  21.  As  predicted  by  Lowke  and  Davies,  secondary 

emission  of  electrons  from  the  cathode  causes  the  I-V  characteristic  to 

rise  more  rapidly  with  an  increased  y  .  However,  at  current  densities 

2 

below  a  certain  value  (.3  mA/cm  In  Fig.  21),  ionization  is  minimal  and 
the  cathode  fall  is  nearly  the  same  as  with  y  =  0. 

All  of  the  I-V  curves  in  Fig.  21  were  calculated  using  the  constant 
recombination  rate  coefficient  given  by  Lowke  and  Davies  (Ref.  11:4994). 
Because  this  value  significantly  overestimates  the  magnitude  of  this 
coefficient  in  the  cathode  sheath,  the  influence  of  including  the  func¬ 
tional  dependence  of  the  recombination  rate  on  E/N  was  investigated. 

Two  functional  forms  were  considered  for  argon,  while  only  one  function 
was  used  for  methane.  One  of  the  argon  recombination  coefficients  was 
for  Ar+  ions,  while  the  other  was  for  A^"*"  ions.  The  total  discharge 
voltages  obtained  using  these  functions  differ  by  less  than  one  percent 
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Cathode  Fall 


Secondary  Emission  Coefficient 


Figure  20.  Cathode  Fall  Versus  Secondary  Emission  Coefficient 
(Argon,  760  Torr,  J  *  500.  mA/cnr) 
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10°  101  io2  io3  io* 


Figure  21.  Calculated  I-V  Characteristics  for  Three  Secondary  Emission 
Coefficients  (Argon,  240  Torr,  S  «  3.6  x  10*®) 
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2 

(only  current  densities  below  48  raA/cm  were  investigated),  and  both  are 
less  than  five  percent  below  the  values  obtained  using  a  constant 
characteristic  of  300°K  ions  (See  Appendix  A) .  Although  this  result 

seemingly  vindicates  Lowke  and  Davies'  assumption  of  a  constant  R^,  the 
functional  form  for  Ar+  ions  was  used  for  many  of  the  calculations  because 
it  can  be  calculated  more  easily  than  the  coefficient  for  AT^  ions. 

ip  w. 

(See  Appendix  A.)  However,  the  recombination  coefficient's  E/N  dependence 
requires  an  iterative  approach  to  obtain  self-consistent  starting  condi¬ 
tions  at  the  positive  column  interfaces  with  the  sheaths.  At  high  dis¬ 
charge  current  densities  (field  strengths),  the  rapid  variation  of  the 
variables  involved  makes  this  calculation  very  unstable.  In  those  in¬ 
stances  in  which  a  self-consistent  set  of  parameters  could  not  be  found 
in  this  way,  a  Large,  constant  Rf  was  used  to  start  the  integration, 
while  subsequent  calculations  in  the  integration  used  the  functional 
form  of  Rr . 

Using  this  technique,  or  using  a  constant  Rr,  the  argon  I-V  char¬ 
acteristics  of  Lowke  and  Davies  were  extended  to  higher  current  densities 
and  expanded  to  include  additional  effects.  These  I-V  characteristics 
were  compared  qualitatively  to  similar  calculations  by  Zakharov,  et.al. 
(Ref.  10:1077),  and  directly  to  experimental  data  by  Bletzinger  (Ref. 
5:1-2).  All  of  these  calculations  were  done  at  760  Torr  with  an  elec¬ 
trode  separation  of  .3  cm,  no  metastable  ionization,  and  with  y  -  .02. 

The  results  obtained  are  shown  in  Figs.  22  and  23.  At  low  current  den¬ 
sities,  almost  the  entire  voltage  drop  in  the  discharge  occurs  in  the 
cathode  sheath  and  the  discharge  current  density  increases  relatively 
slowly  with  increasing  voltage  (region  AB  in  Figure  23) .  At  a  parti¬ 
cular  ''ignition"  voltage,  however,  there  is  a  sharp  increase  in  the 
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slope  of  the  I-V  curve.  Beyond  this  ignition  voltage  the  discharge  cur¬ 
rent  density  increases  rapidly  (region  CD  in  Figure  23)  for  small  increases 
in  the  total  discharge  voltage.  This  sharp  increase  in  the  current  den¬ 
sity  is  due  to  the  rapid  increase  of  Townsend  ionization  in  the  cathode 
sheath  at  these  currents  and  voltages.  In  this  current  density /voltage 
regime,  the  discharge  current  density  is  determined  primarily  by  the  mo¬ 
bility  and  density  of  electrons  in  the  plasma  region,  and  is  approximately 
proportional  to  the  total  discharge  voltage.  As  a  result,  the  electric 
field  in  the  positive  column  must  increase  to  keep  pace  with  the  required 
increase  in  current  density.  (The  number  density  of  electrons,  the  pri¬ 
mary  current  carriers  in  the  positive  column,  is  not  allowed  to  increase 

1/2 

since  it  is  assumed  that  ne  =  ( S/R )  .)  With  a  larger  electric  field, 

the  voltage  drop  across  the  positive  column  also  increases,  becoming  com¬ 
parable  to  the  cathode  fall  voltage.  At  very  high  current  densities,  the 
cathode  fall  becomes  negligible  compared  to  the  positive  column  voltage 
drop.  The  cathode  fall  may  also  be  negligible  at  intermediate  currents 
if  the  interelectrode  gap  is  substantially  increased. 

The  transition  from  the  low  to  high  conductivity  regions  occurs 
smoothly  for  large  external  ionization  source  strengths.  However,  at 
lower  source  strengths,  regions  of  negative  differential  conductivity 
occur,  as  in  curve  BC  in  Figure  23.  In  this  I-V  regime,  the  effective 
resistance  is  negative  as  the  discharge  reconfigures  itself  into  a  more 
efficient  operating  mode.  This  behavior  is  similar  to  that  of  a  sub¬ 
normal  discharge,  as  described  previously,  in  the  transition  to  a 
stable,  glow  discharge.  Also,  the  subsequent  rapid  rise  in  current 
density  is  similar  to  that  observed  in  normal  glow  discharges  (See 
Fig.  1).  These  regions  of  negative  differential  conductivity  at  low 
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current  densities  and  source  strengths  were  also  observed  in  similar 
calculations  by  Velikhov  (Ref.  28:593)  and  Zakharov  (Ref.  10:1077)  for 
Nj  discharges.  Although  Zakharov  maintained  that  this  region  was  not 
physical  and  could  not'  be  measured  experimentally,  Averin  et.al.  reported 
such  a  measurement  in  1980  (Ref.  7).  Zakharov  described  this  region  of 
negative  differential  conductivity  as  an  unstable  "transition  region 
between  a  low-current  discharge  (in  which  the  charges  are  produced  only 
by  the  external  ionization  source)  and  a  discharge  in  which  the  cathode 
sheath  which  is  formed  is  an  electron  emitter"  (Ref.  10:1077).  As  men¬ 
tioned  previously,  this  negative-differential  conductivity  region  vanishes 
at  higher  ionization  source  strengths.  For  large  S,  the  effect  of  Town¬ 
send  ionization  at  low  current  densities  is  reduced  because,  as  seen  in 
Figs.  14,  15,  and  16,  the  cathode  sheath  fieli  and  the  sheath  thickness 
are  also  reduced.  Thus,  large  source  terms  "smooth"  the  transition  from 
external  ionization  to  Townsend  ionization. 

The  Case  III  method  was  used  to  obtain  the  argon  I-V  characteristics 

in  Figures  22  and  23.  The  low  current  density  variation  of  Case  III  was 

2 

used  for  current  densities  up  to  48  mA/cm  ,  while  the  high  current  den¬ 
sity  variation  was  used  at  higher  currents.  Thus,  the  observed  regions 
of  negative  differential  conductivity  are  not  simply  a  result  of  using 
different  variations  of  the  Case  III  method.  Also,  because  only  this 
method  was  used  to  obtain  the  I-V  characteristics,  the  curve  for  S  « 

3.6  x  1014  in  Figure  23  is  incomplete.  At  this  low  source  strength,  the 
Case  III  method  is  not  applicable  at  low  current  densities,  as  explained 
previously.  Thus,  this  characteristic  could  be  completed  at  low  current 
densities  by  using  the  Case  II  method. 
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In  order  to  experimentally  observe  this  negative  differential  con¬ 
ductivity  region,  the  external  circuit  load  line  must  be  shallow  enough 
to  avoid  simultaneous  operation  at  more  than  one  point  of  he  discharge's 
I-V  characteristic  (See  Fig.  3).  To  accomplish  this,  a  very  large  cir¬ 
cuit  resistance  must  be  used,  as  was  done  by  Averin  et.al.  (Ref.  7)  for 
discharges  in  nitrogen. 

Bletzinger  (Ref.  5:2)  was  also  able  to  observe  this  negative  *•.  _*f- 
ferential  conductivity  region  for  discharges  in  argon  at  760  Torr.  The 
I-V  characteristics  which  he  obtained  at  two  different  source  strengths 
(produced  by  an  external  electron  beam),  are  compared  with  the  calculated 
characteristics  in  Fig.  23.  Although  Bletzinger's  electrode  separation 
was  2.2  cm,  while  the  electrode  separation  in  the  calculational  model 
was  only  .3  cm,  it  is  still  reasonable  to  compare  these  characteristics. 

In  Bletzinger's  experiment,  the  positive  column  field  was  measured  to  be 
less  than  .01  Townsend.  As  a  result,  the  cathode  fall  constituted  nearly 
the  entire  discharge  voltage.  In  the  computational  model,  with  an  elec¬ 
trode  separation  less  than  15%  of  the  experimental  model,  most  of  the 
discharge  voltage  drop  will  occur  in  the  cathode  sheath  only  if  the  cur- 
rent  density  is  below  25  mA/cra  .  Thus,  the  inconsistency  between  calcula¬ 
tion  and  experiment  which  is  observed  in  the  positive  column  (i.e.  the 
calculated  positive  column  field  is  much  larger  than  the  observed  field) 
has  a  minimal  effect  when  the  comparison  is  restricted  to  low  current 
densities  and  small  electrode  separations.  As  the  current  density  is 
Increased,  the  calculated  positive  column  field  must  increase  and  the 
unphysical  behavior  of  the  model's  positive  column  in  accentuated. 

Comparing  the  calculated  characteristics  (Fig.  23),  at  source 
strengths  of  3.6  •  10^  and  3.6  •  10*^,  to  Bletzinger's  characteristics 
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Voltage,  kV 

Figure  22.  I-V  Characteristics-Calculated  for  an  Argon  Discharge 

(760  Torr,  y  -  .02,  d  -  .3  cm,  no  metastable  Ionization) 
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Figure  23.  I-V  Characteristics-Calculated  for  an  Argon  Discharge  (expansion  of  low  current  and 
voltage  region  of  Fig.  22.) 


it  is  seen  that  their 


for  source  strengths  of  6.2  •  10^  and  2.82  •  10*^, 
general  shapes  are  quite  similar,  especially  at  lower  currents.  As  the 

current  is  Increased,  the  calculated  positive  column  field  must  increase 
•  \ 
more  rapidly  for  lower  source  strengths.  Thus,  the  calculated  character¬ 
istic  for  S  ■  3.6  •  10^^  rapidly  diverges  from  the  experimental  curves. 
Comparing  the  more  well-behaved  calculated  characteristic  for  S  *  3.6  * 

if  if  (j*.  ** 

10  to  Bletzinger's  2.82  •  10  curve,,  it  is  seen  that  the  slope  of  the 
calculated  high  conductivity  region  is  not  as  great  as  that  observed. 

This  is  partially  due  to  the  previously  discussed  problem  in  the  positive 
column.  However,  it  could  also  be  the  result  of  using  too  small  a  value 
for  the  secondary  emission  coefficient  in  the  calculations.  In  the  cal¬ 
culation  done,  the  secondary  emission  coefficient  was  assumed  to  be  2Z, 
but  this  may  be  considerably  smaller  than  the  actual  value,  since  the 
secondary  emission  coefficient  for  metastable  argon  atoms  (which  may  be 
present  in  large  numbers)  is  40%.  With  a  larger  secondary  emission  coef¬ 
ficient,  the  slope  of  the  I-V  characteristic  is  increased,  as  shown  in 
Fig.  21.  Thus,  a  secondary  emission  coefficient  of  two  percent  probably 
underestimates  the  magnitude  of  this  effect  in  the  discharges  investi¬ 
gated  experimentally.  The  slope  of  an  I-V  characteristic  would  probably 
also  increase  if  the  effect  of  secondary  emission  at  the  anode  were  con¬ 
sidered.  However,  since  such  emission  would  only  occur  for  low  current 
densities  (for  which  the  anode  sheath  field  is  positive),  the  slope 
Increase  would  probably  be  minimal. 

It  may  also  be  observed  in  Fig.  23  that  the  voltage  separation 
between  the  experimental  characteristics  is  much  less  than  the  separation 
of  the  calculated  characteristics.  In  the  calculatlonal  model,  an 
increase  in  S  results  in  an  Increase  in  ng  in  the  positive  column.  For 
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a  given  current  density,  then,  the  electric  field  in  the  positive  column 
must  decrease.  This  results  in  a  weaker  electric  field  in  the  cathode 
sheath  and  thus  the  observed  decrease  in  discharge  voltage  (See  Fig.  17). 
If  there  were  another  process  (in  addition  to  the  external  ionization 
source)  by  which  free  electrons  were  generated  in  the  positive  column, 
and  if  the  number  of  electrons  produced  by  this  process  were  proportional 
to  the  local  field  strength,  the  decrease  in  discharge  voltage  resulting 
from  an  increase  in  S  would  be  less  than  when  S  is  the  only  source  of 
electrons  in  the  positive  column.  The  decrease  in  the  field  which  would 
normally  occur  with  an  increase  in  S  would  also  cause  the  other  source  of 
electrons  to  reduce  its  production,  negating  some  portion  of  the  increase 
in  n£  due  to  S.  Thus,  the  closer  voltage  spacing  of  the  experimental 
characteristics  seems  to  indicate  that  there  is  an  additional  source  of 
electron-ion  pairs  in  the  positive  column. 

Since  this  process  is  probably  operative  in  the  cathode  sheath  as 
well,  the  electric  field  needed  to  support  the  required  current  density 
would  be  reduced.  Thus,  both  the  cathode  fall  and  the  total  discharge 
voltage  would  be  reduced,  resulting  in  a  shift,  of  the  entire  character¬ 
istic  to  lower  voltages,  as  demonstrated  by  the  experimental  character¬ 
istics  in  Figure  23. 

A  very  likely  candidate  for  this  additional  process  is  metastable 
ionization.  Because  argon  has  a  very  high  energy  metastable  state,  even 
the  weak  electric  fields  in  the  positive  column  may  produce  significant 
ionization  from  these  excited  states.  With  the  addition  of  the  effect 
of  metastable  ionization,  an  I-V  characteristic  was  calculated  for  a 
source  strength  of  3.6  •  10*^  in  a  760  Torr  argon  discharge.  This 
characteristic  and  one  of  Bletzinger's  experimental  curves  (for  a  similar 
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source  strength  and  electrode  gap)  are  shown  in  Figure  24.  (The  basic 

I-V  characteristic  in  Figure  24  differs  slightly  from  that  in  Figs.  22 

and  23  because  a  constant  value  of  R  was  used  for  the  calculations 

r 

represented  in  Figure  24,  while  an  Rr  with  a  functional  E/N  dependence 

was  used  for  the  calculations  in  Figures  22  and  23.)  Although  the  meta- 

2 

stable  calculation  was  not  carried  to  current  densities  below  7  mA/cm  , 
the  experimental  and  calculated  curves  are  in  reasonable  agreement  where 
they  overlap.  Thus,  the  addition  of  metastable  ionization  shifts  the 
I-V  characteristic  to  lower  discharge  voltages,  while  leaving  the  slope 
of  the  high-conductivity  portion  of  the  characteristic  relatively 
unchanged.  Consequently,  to  reproduce  the  very  rapid  current  rise  which 
is  observed  experimentally,  the  secondary  emission  coefficient  must  also 
be  Increased. 

Metastable  ionization  may  also  be  used  to  correct  the  inflated  value 
for  the  electric  field  in  the  positive  column.  Assuming  that  the  rate  of 
ionization  of  metastables  is  of  the  same  magnitude  as  the  rate  of  recom¬ 
bination  and  the  rate  of  external  ionization,  the  number  density  of  elec¬ 
trons  in  the  positive  column  is  obtained  using  the  electron  continuity 
equation  (3).  In  the  positive  column,  9je/9x  =  0,  and  the  weak  electric 
field  allows  the  neglect  of  Townsend  ionization.  Thus,  a  quadratic 
equation  in  n  is  obtained: 


S  -  Yn 


+  R.  n  n 
in  e  n 


with  solution 


R.  n  + 
1m  m 


1 


4Sy 


2Y 


(57) 


(A  negative  value  for  the  square  root  would  not  be  physical.)  When 
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Current  Density,  mA/cm 


Figure  24.  I-V  Characteristics-Experimental  and  Calculated  for  an 
Argon  Discharge  (same  parameters  as  Fig.  22,  but  with 
metastable  ionization) 
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me tas table  Ionization  is  significant,  this  equation  will  yield  an  increased 


n#.  For  a  given  current,  then,  the  electric  field  is  reduced.  For  a  meta- 
stable  ionization  rate  of  R^m  **  1  *  10  **exp(-4.15/Te),  a  metastable  number 
density  of  n^  *  2.5  •  10  •  exp(-11.65/Te)  (See  Appendix  A),  a  source 
strength  of  S  *  3.6  •  10^,  a  recombination  coefficient  of  Rf  *>  8.81  • 

10'7,  and  at  an  electron  temperature  of  1  eV,  the  number  density  of  elec- 
trons  in  the  positive  column  is: 


2 

Then,  for  a  current  density  of  25  mA/cm  ,  slightly  above  Bletzinger’s 
largest  current  density,  E/N  in  the  positive  column  is  about  .004  Town¬ 
send.  Therefora,  although  the  rates  and  number  densities  used  are  quite 
crude,  and  the  decrease  in  Te  with  decreasing  E/N  was  not  accounted  for, 
it  is  apparent  from  this  analysis  that  the  minute  positive  column  fields 
observed  experimentally  in  argon  discharges  could  possibly  be  accounted 
for  by  considering  metastable  ionization  occurring  in  the  positive  column. 

Although  argon  is  not  very  attractive  for  switching  applications 
because  of  its  low  dielectric  strength,  the  effects  observed  for  secondary 
emission  and  metastable  ionization  may  be  applied  to  the  modeling  of 
other,  more  useful  gases: (Ref.  5:2). 

One  such  gas  which  is  currently  being  investigated  is  methane. 

Because  of  its  very  high  drift  velocity  and  high  dielectric  strength. 

It  seems  to  be  well  suited  for  switching  applications.  I-V  character¬ 
istics  were  calculated  for  methane  discharges  at  760  Torr.  In  Figure  25, 
these  calculated  characteristics  are  compared  against  experimental  mea¬ 
surements  made  by  Bletzinger  (Ref.  11:2).  Although  the  calculated  curves 
are  shifted  to  higher  discharge  voltages,  the  spacing  between  different 


source  strength  curves  in  nearly  the  same  for  both  pair  of  experimental 
and  theoretical  curves.  In  addition,  both  sets  of  curves  exhibit  transi¬ 
tion  regions  of  negative  differential  conductivity,  although  these  regions 
are  accentuated  both  in  voltage  and  current  density  for  the  calculated 
curves.  However,  it  should  be  noted  that  these  discrepancies  are  rela¬ 
tively  minor,  considering  the  inaccuracies  inherent  in  the  experimental 
measurements.  For  example,  the  values  assigned  to  the  experimental  source 
strengths  are  only  approximate.  These  values  were  assigned  based  on 
average  stopping  powers,  while  a  more  exact  calculation  requires  a  statis¬ 
tical  approach,  using  Monte  Carlo  techniques.  Thus,  the  values  for  the 
experimental  source  strengths  may  be  in  error  by  up  to  an  order  of  magni¬ 
tude.  Thus,  if  the  actual  experimental  source  strengths  are  a  factor  of 

ten  larger  than  those  listed  in  Fig.  25,  the  calculated  and  experimental 

18 

curves  for  S  -  1  x  10  would  nearly  coincide.  This  explanation  would 

also  reconcile  the  discrepancy  in  the  voltage  spans  of  the  negative  dif- 

* 

ferential  conductivity  regions .  :.:4s  seen  previously  for  argon  discharges 
(Fig.  22),  with  an  increase  in  S,  the  voltage  span  of  this  region 
decreases  and  eventually  disappears.  Thus,  if  the  experimental  source 
strengths  are  in  error  as  indicated,  the  voltage  spans  of  the  negative 
conductivity  regions  would  be  comparable  for  the  experimental  and  cal¬ 
culated  characteristics.  The  exaggerated  extent  in  current  density  of 
the  calculated  negative  conductivity  regions  could  possibly  be  the  result 
of  an  Inflated  secondary  emission  coefficient.  The  value  of  two  percent, 
which  was  used  in  these  calculations,  may  be  much  too  large.  As  shown 
in  Figure  21,  a  decrease  in  y  decreases  the  slope  of  the  high  current 
density  region  of  an  I-V  characteristic.  Such  a  decrease  in  y  could 
cause  a  calculated  I-V  cui*e  Lu  return  more  quickly  to  a  positive  slope 


Current  Density,  mA/cm 


figure  25.  I-V  Characters  tics-Experimental  and  Calculated  for  a 
Methane  Discharge  (760  Torr, y  -  .02,  d  »  2.2  cm) 
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as  Che  current  density  is  increased.  Thus,  it  may  be  that  the  actual 
value  for  the  secondary  emission  coefficient  in  a  methane  discharge  ia 
much  less  than  the  two  percent  figure  used  in  the  calculations.  It 
should  be  noted  that  these  discrepancies  between  calculation  and  experi¬ 
ment  cannot  be  explained  by  the  addition  of  metastable  ionization,  as 
was  done  for  argon.  For  one  thing,  only  the  rare  gases  and  some  mole¬ 
cular  gases  have  such  high  energy  metastable  states.  In  addition,  meta¬ 
stable  ionization  would  also  cause  the  experimental  characteristics  to 
be  more  closely  spaced  than  the  calculated  characteristics,  but  this  is 
not  observed  to  occur. 
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IV.  Conclusions  and  Recommem  tions 

The  behavior  of  argon  and  methane  gas  discharges  was  Investigated. 
Particular  attention  was  paid  to  processes  occurring  in  the  cathode  sheath. 
It  was  found  that  these  sheaths  became  strong  electron  producers  at  volt¬ 
ages  above  an  ',ignition,,  voltage  characteristic  of  the  discharge  gas  and 
the  external  circuit  parameters.  This  strong  electron  production  results 
in  very  high  discharge  conductivities.  The  I-V  characteristics  calculated 
for  argon  agree  qualitatively  with  those  observed  experimentally,  parti¬ 
cularly  if  the  comparison  is  restricted  to  the  cathode  sheath  region.  In 
order  to  decrease  the  large  calculated  positive  column  fields  to  experi¬ 
mentally-observed  levels,  it  is  hypothesized  that  metastable  ionization 
makes  a  significant  contribution  to  the  ionization  in  the  positive  column. 
Metastable  ionization  also  shifts  the  calculated  I-V  characteristics  to 
lower  voltages,  which  are  observed  experimentally.  The  experimental  curves 
also  have  a  larger  slope  in  the  high  conductivity  region  of  the  character¬ 
istics.  This  larger  slope  can  be  obtained  in  the  calculations  by  the 
inclusion  of  secondary  emission  at  the  cathode.  Metastable  ionization  does 
not  seem  to  play  a  role  in  methane  discharges.  The  calculated  I-V  char¬ 
acteristics  obtained  for  this  gas  agree  very  well  with  those  measured 
experimentally.  The  shift  in  the  calculated  characteristics  to  somewhat 
higher  voltages  may  be  due  to  an  underestimation  of  the  external  ionization 
source  strength  produced  in  the  comparison  experiments. 

These  observed  effects  in  argon  and  methane  discharges  may  be  applied 
to  the  design  of  more  efficient  high-energy  switches.  Some  of  the  desirable 
features  of  such  a  switches  include:  high  conductivity  when  "on,"  very  low 
conductivity  when  "of h±g<>  :£sii  (ratio  of  energy  switched  to  switching 
energy,  i.e.  the  energy  expended  by  a  sustaining  electron  beam),  and  rapid 
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1 

i 

switching.  To  achieve  high  conductivity,  gases  with  high  drift  velocities  j 

should  be  used.  However,  a  gas  with  a  lower  drift  velocity  can  also  pro¬ 
duce  a  highly  conductive  discharge  if  high  y  coatings  are  used  on  the 

\ 

cathode.  These  coatings  would  also  allow  the  use  of  lower  external  ioni¬ 
zation  source  strengths  to  achieve  the  same  high  conductivity.  (At  lower 
S,  the  I-V  characteristic  above  the  ignition  voltage  has  a  shallower  slope 

ft  .*• 

than  that  at  higher  S.  This  larger  slope,  i.e.  conductivity,  can  also  be 
achieved  by  increasing  y  for  lower  source  strengths.)  Thus,  the  switch's 
gain  would  be  increased  without  affecting  the  conductivity.  In  order  to 
decrease  the  time  required  for  a  discharge  to  turn  "off,"  additional  loss 
processes,  such  as  attachment,  could  be  produced  by  the  addition  of  small 
quantities  of  attaching  gases.  These  materials  would  quickly  eliminate 
any  free  electrons  in  the  discharge  once  the  external  ionization  source 
is  turned  off.  To  decrease  the  time  required  for  a  discharge  to  turn  on, 
external  excitation  could  be  used  to  quickly  populate  metastable  states 
of  a  discharge  gas.  This  would  cause  the  discharge  to  immediately  transi¬ 
tion  to  an  I-V  characteristic  at  much  higher  currents.  (The  amount  of 
this  shift  is  determined  by  the  external  circuit  because  the  operating 
points  of  a  discharge  are  at  the  intersections  of  the  circuit's  bad  load 
and  the  I-V  characteristic (s)  of  the  discharge. 

Although  this  investigation  answered  some  questions  about  electron 
beam  sustained  gas  discharges,  it  also  suggested  a  number  of  questions 
and  proposals  for  further  Investigation. 

1)  I-V  characteristics  for  a  760  Torr  argon  discharge  should  be 
recalculated  in  an  attempt  to  reduce  the  discrepancies  between  the  cal¬ 
culated  and  the  experimental  characteristics.  These  calculations  should 
include:  metastable  ionization,  in  both  the  cathode  sheath  and  the 
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positive  column;  more  accurate  values  for  the  rate  of  this  Ionization  and 
for  the  number  density  of  metastables;  more  accurate  values  for  the 

secondary  emission  coefficient  (taking  into  account  the  large  coefficients 

\ 

associated  with  metastables) ;  and  source  strengths  which  are  the  same  as 
those  used  by  Bletzinger  (Ref.  5). 

2)  Experimental  investigations  of  argon  discharges  should  be  extended 
to  lower  source  strengths  to  determine  if  the  I-V  characteristics' ootalned 
in  the  first  proposal  are  physically  realistic.  In  particular,  the  region 
of  negative  differential  conductivity  should  be  explored. 

3)  The  I-V  characteristics  calculated  for  methane  should  be  extended 
to  higher  current  densities  and  a  wider  range  of  source  strengths.  For 
these  calculations,  a  more  accurate  value  for  the  actual  secondary  emission 
coefficient  shot.ld  be  obtained.  Also,  the  validity  of  the  method  for 
determining  the  source  strength  produced  by  an  external  electron  beam 
should  be  investigated. 

4)  High  y  cathode  coatings,  which  could  reduce  the  source  strength 
required  to  access  the  rapid  current  rise  portion  of  discharge  I-V  char¬ 
acteristics,  should  be  experimentally  explored.  These  coatings  would 
increase  the  gain  of  an  externally-sustained  gas  discharge  switch. 

5)  I-V  characteristics  should  be  calculated  for  a  range  of  gas 
pressures.  The  effects  produced  may  be  a  function  of  the  discharge  cur¬ 
rent  density.  There  may  be  an  optimal  operating  pressure  for  a  given  gas 
discharge  switch. 

6)  The  diffusion  coefficient  of  electrons  used  in  the  computer  cal¬ 
culations  should  be  modified  to  correct  for  the  inflated  value  of  this 
coefficient  when  electrons  diffuse  into  a  retarding  electric  field. 

7)  Because  of  the  importance  of  rapid  "off"  switching,  the  effect 


of  small  quantities  of  attaching  gases  in  externally-sustained  discharges 
should  be  determined  by  calculating  the  I-V  characteristics  of  these 
mixtures..  Initially,  a  mixture  of  argon  and  oxygen  should  be  investigated 
because  Bletzinger  (Ref.  5)  has  already  measured  I-V  characteristics  for 
such  a  discharge. 

(/  ^  ' 
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Appendix  A 

Transport  Coefficients  and  Material  Functions 

Townsend  Ionization  Coefficient  (a) 

Argon:  a  =  N  •  2.9  •  10~^7  exp (-1.48  •  10  ^  •  N/E)(cm2)  This  ex¬ 
pression  for  a,  which  was  obtained  from  experimental  measurements  (Ref. 
35:515),  was  used  by  Lowke  and  Davies  (Ref.  11:4994).  This  expression 
was  used  for  Case  II  and  Case  III.  However,  recent  experimental  data 
indicates  that  this  expression  should  be  modified  at  low  E/N  (Ref.  32:167). 
This  modified  expression,  as  used  in  Case  I,  Js: 

o  -  N  1.76275  •  10-17  exp(-l. 45027  •  10-15  •  N/E) 

Methane :  Tabulated  experimental  values  if  Cookson  et.al.  were  used 
for  the  net  ionization  coefficient,  (a  -A),  which  implicitly  accounted 
for  attachment  occurring  in  methane  (Ref.  14:. .8).  Values  for  (a  -  A)  at 
E/N  <  72Td  (1  Townsend  =  1  Td  =  1  •  10  17  Volt-cm2)  were  obtained  by 
logarithmic  extrapolation. 

Recombination  Rate  Coefficient  (R_) 

. .  . r — 

Argon:  Rr  *  7.6  •  10  •  e  *  (cm  /sec)  This  expression  for  the 

recombination  rate  of  Ar2+  ions,  given  by  Lowke  and  Davies,  was  determined 
experimentally.  It  requires  a  knowledge  of  the  average  electron  energy 
e"  (in  eV)  as  a  function  of  E/N.  This  function  was  used  only  for  compari¬ 
son  with  another  function  (Ref.  11:4994). 

-7  3 

Rf  ■  8.81  •  10  (cm  /sec)  This  constant  value,  which  was 
used  by  Lowke  and  Davies,  was  obtained  from  the  expression  above  for  an 
average  electron  energy  of .  This  greatly  exaggerated  the  recom¬ 

bination  rate  at  higher  electron  energies,  and  thus  higher  fields,  i.e. 
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In  Che  cathode  sheath.  However,  since  this  rate  is  multiplied  by  n#, 
which  decreases  rapidly  away  from  the  plasma  region,  Lowke  and  Davies 
maintain  that  it  is  not  too  inaccurate  (Ref.  11:4994).  This  value  was 
used  in  the  present  analysis  for  code  verification  and  for  I-V  comparisons 
at  different  source  strengths  and  with  the  addition  of  secondary  emission. 

*  7.5  •  10  ®  •  (E/N)  *73  (cm3/sec)  This  analytic  expres¬ 
sion  for  the  recombination  rate  of  Ar+  ions,  given  by  Dzimlanski  and  Kline, 
was  used  for  most  of  the  argon  calculations  (Ref.  14:13).  The  value  for 
E/N  in  this  expression  must  be  given  in  Townsends.  Although  the  region  of 
validity  of  this  function  is  from  .1  to  10  Td,  cathode  sheath  fields  in 
this  investigation  sometimes  reached  100  Td.  Although  Lowke  and  Davies 
maintain  that  the  use  of  the  kx^  Rr  is  more  appropriate,  the  total  dis¬ 
charge  voltage  obtained  using  Lowke  and  Davies '  analytic  Rr  was  less  than 
5Z  different  from  the  voltage  obtained  using  I zimianski  and  Kline's  Rr . 

Methane :  Rr  «  2.4  •  lO"”7  •  (E/N)-"34  (cm3/sec)  This  expression  for 
the  recombination  rate  in  methane  was  obtained  from  Dzimlanski  and  Kline. 
The  value  for  E/N  must  be  given  in  Townsends.  This  function,  which  was 
only  valid  from  .2  to  10  Td,  exaggerated  %r  at  the  low  fields  in  the  posi¬ 
tive  column,  although  it  probably  did  not  appreciably  affect  the  calcula¬ 
tions  in  the  cathode  sheath,  as  indicated  in  the  preceeding  discussion. 

This  insensitivity  seemed  to  vindicate  Lowke  and  Davies'  assumption  of  a 
constant  Rr  (Ref.  11:4994). 

Drift  Velocity  of  Electrons  (VQ 

Argon:  Tabulated  values  for  Wg  (as  a  function  of  E/N),  calculated  by 
Engelhard t  and  Phelps  (Ref.  31:A377),  were  used  in  Lowke  and  Davies'  study 
(Ref.  11:4995)  and  in  the  present  analysis. 


Methane:  Tabulated  values  for  W  (as  a  function  of  E/N)  were 
■ "  —  ■  1  e 

obtained  from  calculations  and  experiments  by  a  number  of  Individuals 
(Refs.  22;  23;  24).  The  values  used  In  this  analysis  were  compiled  and 
presented  by  Dzlmlanski  and  Kline  (Ref.  14:17).  Values  at  low  E/N  were 
obtained  by  logarithmic  extrapolation. 

Longitudinal  Diffusion  Coefficient  of  Electrons  (D^) 

Argon:  Tabulated  values  for  D  /k  (volts)  were  obtained  from  the  cal- 
culations  of  Lowke  and  Parker  (Ref.  21:307).  Values  at  higher  field 
strengths  were  obtained  from  Lowke  and  Davies'  extension  (Ref.  11:4995) 
of  the  previous  calculations.  In  order  to  determine  D^,  values  for  the 
mobility,  ke,  were  obtained  from  the  similarly  tabulated  values  of 
versus  E/N. 

Methane:  Tabulated  values  for  D_/k  (volts)  were  obtained  from  the 
-  e  e 

calculations  of  several  authors  (Refs.  23;  24;  26).  Values  at  low  E/N 
were  obtained  by  logarithmic  extrapolation  of  the  values  compiled  by 
Dzlmlanski  and  Kline  (Ref.  14:17).  In  order  to  determine  Dg,  values  for 
k@  were  obtained  from  the  similarly  tabulated  values  of  versus  E/N  for 
methane. 


Metastable  Ionization  Rate  (Rim) 

-8  1/2 

Argon:  R^m  ■  5  •  10  •  Te  •  exp (-4. 15 /Te)  This  rate  was  obtained 

by  shifting  the  cross-section  for  ground  state  ionization  by  lower  electron 
energies.  Thus,  for  argon,  with  an  ionization  energy  of  15.8  eV  and  a 
metastable  energy  of  11.65  eV,  the  metastable  ionization  energy  is  4.15  eV. 
Thus,  the  rate  of  metastable  ionization  is: 


Rim  "  Weoe 


-4.15/Te 
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Assuming  a  Maxwellian  electron  velocity  distribution  and  a  cross-section 
linear  in  Te  (Ref.  34) 


R.  s  (6.7  •  107  •  Te1/2)  •  (2  •  l(f17  •  Te) 

lm 

R.  r  5  •  10'8  •  T.1/2  •  .-4'15/Te 


(2  + 


4.15.  -4.15/Te 

iT")e 


(Te  is  the  electron  energy  in  eV.) 

However,  this  rate  is  very  conservative.  Using  the  metastable  ionization 

-15  2 

cross-section  calculated  by  Duke  (o  =  1  •  10  cm  ) ,  the  rate  obtained  is 
about  20  times  that  shown  above  (Ref.  36).  This  higher  rate  was  used  in 
the  estimation  of  the  contribution  of  metastatle  ionization  in  the  positive 
column. 

Metastable  Number  Density  (n^) 

Argon:  n^  •  (Nq  -  n+)  •  exp (-11. 65 /Te}  The  number  density  of  meta¬ 
stable  argon  atoms  was  approximately  by  using  a  Boltzmann  distribution  of 
energy  level  population  for  a  total  number  density  of  (Nq  -  n+)  unionized 
argon  atoms.  (The  metastable  excitation  energy  for  argon  is  11.65  eV  and 
Te  is  the  electron  energy,  also  in  eV.) 

Electron  Temperature  (T&) 

Argon  and  Argon/Oxygen:  The  average  electron  temperature  was  deter¬ 
mined  as  a  function  of  E/N  for  pure  argon.  These  values  were  obtained 
from  Boltzmann  calculations  done  with  experimentally  determined  cross- 
sections  (Ref.  36). 
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Computer  Codes 
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ooo  ono  n  n  n  o  |o  n  n  o  non  ,o  o  o  loon  no 


PROGRAM  CAStK  INPUT. OUTPUT.  TAPE3  -  OUTPUT.  TAPE9V) 

DIMENSION  YQ(5).wl(2Q).J..HJt20).XI(20).T(150).JE(lSO).N6(150> . 

C  JP(1*0|»£(150).V(150).NI(150) 

COMMON/ CONST/ S. NO .EPS  1. A_PHAtGAMNA.ECHG.MU • MUE.K. EFLD.DE. T.O 
COMMON/ Gfc AK 9/HUSfcO, NUUSE J.NSTEP.NFE.NJE 
COMMON/ PLT/T. JE.Ne.Nl. JP.E . V.VMIN.M 
EX TERNAL  DIFFUN.PEDERV.W 

R£AL  ME, NO. NP. MU. MJc. MG.  MEM.JE.JP. NE.NI  _ v _ 

DISSPLA  PLOT  OPTION  (l-YES  0-N3) 

DPLT  *  1.0 


SET  CONSTANTS  FOR  GEAR  PROGRAM 
NUMBER  OF  OIFF  EONS  IN  SYSTEM  /*> 

N  ■  5 

LOCAL  ERROR  TOLERANCE  PARAMETER 

EPS  -  .0001 

STEP  SIZE  IN  GEAR 

HO  -  .0000033 

GEAR  METHOO 

MF  -  12  ~ . .  ~  '  ~  “  "  * . . 

SET  PHYSICAL  CONSTANTS  FOR  PROBLEM 

ELECTRONIC  CHARGE  (COULOMB) 

ECHG  »  1.6E-19  . " . .  ~ 

P6RMITIVITY  OF  GAS  ( COUL 3MB/ VOLT/ CM) 

EPS!  •  8.B5E-14  ” 

ELECTRON  NASS 

ME  -  .911E-27  ” 

BOLTZMANN  CONSTANT 

BR  -  1.3BE-1B  '  ” 

TRIAL  ANOOE  E  FIELD  IV/ CM) 

EFLO  •  -125.  “ 

C  -  4.2 E  — 8 


L  COUNTS  THE  NUMBER  OF  TRIAL  ANODE  E  FIELDS  USED 

L  -  0  '  ‘  .  .  '■■■ 

<  USED  TO  NEGATE  OIFF  EONS  IN  OIFFUN  ANO  PEOERV 
K  -  l 

SGN  USED  LATER  IN  SELECTION  OF  ANODE  E  FIELD  TO  PR00UC6  NE  £3  26' 
AT  CATHODE 
SGN  -  1.0 


,  SET  PARAMETERS  FOR  DISCHARGE 

30  CONTINUE 
C  INOEX  FOR  FIRST  CALL  T0~5XAR 
INDEX  •  l 

C  DISTANCE  BETWEEN  ELECTRODES  (CM)  1 

0  -  .33 

C  OISTANCE  STEP  SIZE  " 

DELT  -  .01*0 

C  POINT  AT  WHICH  lNIflAL~CaN3"ITrarNS~~6lV€N 

TO  -  0.0 

C  FIRST  OUTPUT  POINT  OESIREO  ' 


non  no 


C 1/CH2/SEC1 


TOJT  -  .01*0 

ELECTRON  CURRENT  UtNSITY  AT  ANOOE  (  1/CM2/SEC1 
Y0(ll  -  C/fc'CHG 

ELECTRON  NUMBER  Oe  NS  ITT  AT  ANODE  (  1/CM2/SEC 1 

ro(2)  -  o.o 

ION  CURRENT  DENSITY  AT  ANODE  ( 1/CM2/SEC1 
Y3(31  -  0.0 

ELECTRIC  FIELO  AT  ANODE  (V/CM1 
Y0(  41  -  E  FLO 

POTENTIAL  AT  ANODE  ARBITRARILY  SET  TO  ZERO 
Y0(51  -  0.0 

CAS  PRESSURE  (0YNE/CM21 
PRES  -  240. *( 1.0L325E6/760. ) 

CAS  TEMPERATURE  (K1  _  _ 

TC  -  273.  . 

ELECTRON  TEMPERATURE  ( K ) 

TE  -  II 60S. 

NUMBER  DENSITY  OF  CAS 
NO  -  PRES/BK/TG 

IONIZATION  SOURCE  STRENGTH  (ION  PAIRS/CM3/SEC1 
S  -  8.3E11 

POSITIVE  ION  MOBILITY  (CM2/SEC/V)  _ _ 

MU  -  A. 9E I9/N0 
RECOMBINATION  COEFFICIENT 
GAMMA  -  8.31E-7 

CALCULATE  ION  NUMBER  OENSITY  AT  ANOOE  _ 

NP  -  YD ( 3  I  /MU/ YQ 1 41 

CALCULATE  CURRENT  DENSITY  (AMPS/CM21  __ 

C  -  ABS(YOIl)-  Y  0( 3  1 1 *EC MG 
WRITE  (3,10001 
WRITE  (3,10011 

WRITE  (3,10021  T0,Y0(11,Y0(3I,C,Y0(2I,NP,Y0(A1, 


Y0(  51 


N  COUNTS  THE  NUMBER  OF  I NTGRATION  PTS  BETWEEN  ANOOE  AND  CATHOOE 
M  ■  I 

LOAD  BOUNDARY  CUNDITIQNS  INTO  PRINT  ARRAYS 
TCM1  -  0  -  TO 
OE ( Ml  -  Y0111 
" NECMi  •  Y0(21  " 

NICMl  -  NP 

'  0PCM1  «  Y0(  31  •  ■“  . 

6IN1  -  Y0CA1 

V( N 1  -  Y0(5i  ” 

SET  INITIAL  COMPARISON  VALUE  FOR  VOLTAGE 
VMIN  •  1000.  _  _ _ _ 

ii 

BEGIN  LOUP  TO  SOLVE  SYSTEM  OF  EQNS  AT  SEVERAL  POINTS  BETWEEN 
ANOOE  AND  CATHOOE 
AO  CONTINUE 

M  «  M  ♦  1  "  "  r 

CALL  ORIVE  (N,T0*H0,Y0,T3UT «£PS,NF« INDEX 1 

c  -  abs ( yo( ii  -  yo(3h«e:hg  __  _ 

NP  -  ABS( Y0( 31/MU/ Y0( 4)0 

_ DRIFT  •  W(Y0(A1/N01  _ _ _ _ 

MUE  -  DRIFT/Y0C41 


2 


I 


QIFF  •  OLC YOC 4 l/NO I 

DMJ  •  0  IFF/MUE  _ _ 

WR I  TEC  3  *  lOOd  I  M. ALPHA  .HOE, DRIFT  ♦  31 FF *DMU .GAMMA 

_  1008  FORMAT  C 12X* I J * 9X « TC1PE9 . 1* 6X 1 1  _  _ 

C 

c  fill  arrays  4ITH  OUTPUT  CANO  write  oyer  them  until  find  the  valu 

C  FOR  THE  ANODE  E  FIELD  WHICH  PkJUJCES  NE  .£Q.  ZERO  AT  THE  CAThODc 

_  c  _  _ _ _ _ _ _ 

TIN)  -  0  -  TOUT  . 

JECMI  -  YOtl) 

NECrt)  «  YOU) 

NIC  Ml  »  Np 
JPCMI  •  YOC  3 ) 

_ EC  M)  -  YOC  A)  _  _  _ 

. .  VC  HI  -  YOC  5)  .  ~  **  . . 

C  OVERALL  POTENTIAL  MINIMUM  WILL  8E  SUBTRACTED  FROM  POTENTIAL  AT 

C  EACH  POINT  TO  MAKE  POTENTIAL  AT  CATHQOE  ZERO 

_ _  _  VMIN  -  AMINIC YOC 51 iVMIN) 

WRIT6C 3.10021  TQUT.YOCll»YO(3)»C*YOC2l*NP«YOC4l*YOC5)*NQJSEO 

_  _ C  NORMAL  CONTINUATION  OF  INTEGRATION  WITH  INDEX  .Ea.  ZERO _ 

IF  CIN3EX.EJ.0)  GO  TO  100 
C  ELSE#  ERROR  TERMINATION  ANO  EXIT 

WRITE  C  3 • 1003 i  INOEX 
_ _  GO  TO  200 

C  LOOP  BACK  TO  LINE  AO  UNLESS  HAVE  REACHED  CATHODE  "  ~ 

C  OECREASE  STEP  SIZE  AS  GET  NEAR  CATHODE 

100  IF  CA3SC0ELT).LT.A8SC5.E-**O)>  GO  TO  120 
IF  CCD-TOUTI.lt. C2.*0£LT||  OELT  -  DELT/2. 

IFCCO-rOUTl.LT. COELT I I  OsL T  -  OElT/10. 

TOUT  -  TOUT  ♦  DELT 
GO  TO  AO  * 

_ _  120  WRIT  EC  3  *  100 A)  NSTEPtNF E .NJE _ 

C  DO  OISSPLA  PLOTS  OF  FIRST  ANO  EVERr  TENTH  ITERATION 

«  2F{ COPLT.EQ.l. | .ANO. CCL.eO.Ol.OR. CMOOCLvlOl .EJ.O I) I  CALL  PLOTClI 

C 

g.  - - - 

C  NOW*  CHECK  RESULTS  OF  CALCULATION  TO  SEE  IF  NE  AT  THE  CATHODE  IS 

C  ZERO  -  IF  NOT*  MODIFY  E  FIELD  AT  ANODE  AND  REPEAT  ENTIRE 

__ _ C _ I NTEGR A T ION  B  Y  LOUP  IMG  8  A CX  TO  LINE  30 _ 

•  L  •  L  ♦  l 

1  IF  CL. EO. 11  GO  TO  130  “  . 

_ _ I F  C  (_NEX  1  *  Y  0  C  2  ).LT. 0.01. OR. CNEGFL6.60.il  I  GO  TO  1A0 _ 

J  C  _  BEGINNING  AT  AN  E  FIELD  VALUE  WHICH  ASSURES  THAT  NE  AT  THE  CATrtU 

-  ~  C  WILL  8E  POSITIVE*  THE  E  FIELD  IS  INCREMENTED  CnEGATIVcLYI  UNTIL 

'  C  GOES  NEGATIVE  -  THE  E  FIELD  IS  THEN  INCREMENTED  POSITIVELY  USING 

_  c  yHE  SCN  variable  -each  succeeding  Increment  is  halved  and  is  ma 

C  POSITIVE  OR  NEGATIVE  BASSO  UPON  WHETHER  NE  CHANGES  SIGN  OR  NDT  - 

•M  C  WHEN  THIS  INCREMENT  BECOMES  LESS  THAN  ONE  PERCENT  OF  THE  E  FIELD 

C  VALUE*  THE  OPTIMIZATION  .03?  IS  EXITED  TO  LINE  300_ 

c  "•  *"  . . " 

130  ME  XI  -  Y0C2I 
EFlOI  -EFLO 

_  _ IF  CABSCY0C2)  I.LE.L.E4I  30  T0JI03 _  _ _ 

EFLO  -  EFLO  -  10." 
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DELE  -  EFLO  -  cFLOl 

GO  TO  30  _ 

"140  DclE  -  AdS!0ELE/2.) 

NEGFLG  -  l 

IF  (N£<l*YO!2).LT.J.O>  SGN  •  —SGN 
NE* 1  -  YO! 2) 

EFLOl  •  E FLO 

_  EF..D  -  EFLO  -  DELE *S GN  _  > 

IF  !A8S!Y0!2) ) .LE.I.E4)  GO  TO  *00 
IF  (A8S(0ELE).LT.AdS(.01*EFL0n  GO  TO  300 
GO  TO  30 

200  WRITE  !  3  , 1004 )  NSTEP • NFE  » N  JE 
GO  TO  500 

300  CONTINUE  _  _ _ >; _ 

350  CONTINUE  ' 

500  CONTINUE 

1000  FORMAT (  lHl,3X"0IST4NCEHl3X"jeM13X"J^"?X"CURRENT  DENSITY" iX» 
C"NE"13XHN>"I1X"E  FI£LQ"8X,'VQLTAGE"4X"QRDER,,) 

1001  FORMAT!  bX"CMM10X,,l/C«2/SEC"bX"l/CM2/SEC"6X,,AMPS/CH2"8X, 

C*l/CM3,*10XMI/Cn3"dX,*V0LrS/Ctt"9X"Y0LTS"///)  _ _ 

1002  FORMAT! 8! 1PE 12. 4, 3 X), 1 4) 

1003  FORMAT! //26H  ERROR  RETURN  WITH  INDEX  -,I3//) 

1004  FORMAT ! // 21H  PROBLEM  COMPLETED  IN,I5.bH  STEPS/ 

C  2 IX, I  5  «  L4H  F  EY ALUATIONS/ 

C  21X, 15* 14N  J  EVALUATIONS///) 

CALL  exit  _ _ _ _ _ 

"END  “  . . .  '  . .  ‘  . . . ” 


SUBROUTINE  DIFFUN! N,T,Y,YOQT) 

COMMON/CONST/S*  NO,  EPS  I ,  ALPHA, GANNA, ECHG. MU,  HUE ,J(, EFLO. J3E-T,0 _ 

"  DIMENSION  YIN), YOOTIN) 

_ EXTERNAL  W,OL  _  _ _ : _ _ 

REAL  MU, NO 

_ IF  !1.48E-15*ND/A8SIYt4)).GT.300.)  GO  TO  10 

ALPHA  -  N0*l.7b2?5E-17*E(Pt-l.45027E-15*N0/ABSiY!4) ) ) 

_  GO  TO  20  _ _ * _ _ 

T6  ALPHA  V  'o.'O 

20  CONTINUE  _ _ 

C 

C  SIGN  CONVENTION  IS  THAT  THE  DRIFT  VELOCITY*  W,  IS  POSITIVE 

C  WHEN  THE  ELECTRIC  FIELD, E,  IS  NEGATIVE 

___  _  SGN  •  -1.0  _  _  _  _  _ 

IF  IYI4I.LT. 0.0)  SGN  -  l.O' . 

C  * 

C  INITIAL~"VALUES  FOR  DIFFERENTIAL  EQUATIONS  AT  ANODE  I  ROUNDOFF 

C_  ERRORS  MOULD  PREVENT  YOUTH)  AND  YD0TI3)  FROM  BEING  EXACTLY  S) 
IF  (T.EQ.O.O)  GO  TO  50 

GO  TO  100  _  _  ___  .. 

50  YOOTtl)  -  S  '  * 

Y00TI2)  -  — Y! 1)/0LIYI4)/-N0) 

Y00TI3)  -  S  . . 

YOOT 14)  -  0.0 
Y00T(5)  -  -EFLO 
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GO  TO  200  _ 

100  CONTINJE  . .  " 

TOOT  11)  -  S*Y(2)*AlPHA*A3S(M(Y(4) /NO)  )-GAMM A*Yf 2 ) *Y ( 3 ) / Yl  4)/MU 
Y03TC2)  -  (Y(2)*rt(Y(4)/N3)/SGN  -  Y ( 1 ) ) /OL ( Y i 4) /NO) 

Y00TI3)  -  Y00TC1) 

Y03T ( 4 )  -  (6CHG/EPSl)*(Y(3)/Y(4)/HU  -  YI2)) 

Y00TI5)  -  -Y(4) 

OIFF  t  ON  S  ARE  NEGATED  TO  ALLOW  F3R  INTEGRATION  FROM . 

SHEATH  dACK  TO  CATHODE 
200  CONTINUE  '  '  ^  ‘  ~ 

YD3T( 1)  -  YQ3T ( i)*( — 1 ) **< 

YOOT ( 2 )  -  Y0UT(2)*<-1)**< 

YD3TC3)  -  Y00TC3)*t-l)**< 

YOOT (4)  -  Y00T<4)*(-1)**<  . . 

YOOT  <  5 )  -  Y0UT(5)*(— !)•*< 

RETURN  . . . 

END 


FUNCTION  VilXlN) 

THIS  FUNCTION  CALCULATES  ELECTRON  DRIFT  VELOCITY,! CH/SEC).  FOR  __ 
PARTICULAR  VALUES  OF  E/N  IGRIO  VALUES  FOR  ORIFT  VELOCITY  ARE 

FROM  ENGELHART  AND  PHELPS)  _ _ _  _ _ _ 

DIMENSION  WI( 33 ) tXI (33) 

COflrtON/ CONST/  S  •  NO*  EPS  I  *  ALPHA  t  GA1MA*£CHG  tNU  f  HUE  «  <  _ 

REAL  NO  , MU* HUE 

DATA  FROM  ENGELHAROT  A NO  FHELPStPHYS  REV, ( 133) , JAN  64.PA377 
DATA  *1  / 3169. 2. 4*94.3 *76*7 .9* 13002. ,22293.,  30209. , 5 369 1 . , 5369 1 .  , 
C 72 756., 90267., 104566., 11 309?., Lid  391., 133171., 148334.,  1/3  524.  , 
C195821., 222439., 250207., 291743., 343518., 390211., 450547., 54456 5.  , 
C6930 72., 1029840. *1350700. ,1731450 .*2323430. , 3037390. , 39 70b40.  , 

C5 61 42 30., 6524640./ 

DATA  XI /5 .37656-2 1,9. 1264 E-2L, 1. 4096-20*2.1 5436-20  »  3.02S2E-20, 
C3. 3814c -20, 4. 2638E-20 ,4. 2638E— 20, 4.8117E- 20 ,6.0674E— 20, 8.4372E-2 
C1.2126E-I9,i.7595t-19*2. 51226—19*4.305 36-1 9,8.66 16 E— 19, 1.46056-1 
C2. 479LE-1 8 ,4. 405*66-1 8*8.  24716-18,1.64826-1 7,2.98686— 17,4.84236-1 
C6. 16746-17,7. 75136-17,1. 1396-16* l. 55926- 16, 2.02496-16, 2. d9076- la 
C3. 95556-16*5.5 1826-16  *  8. 4372E-16* 9.93576-16/ 

XIN  -  ASS  ( XI N )  _  _  _  

CALL  INTRP(XI*NI,33*2*XI N*YQUT) 

*  •  YOU! 

RETURN 


FUNCTION  OL  (  X  I  r«  I 

THIS  FUNCTION  CALCULATES  TrtE  LONGITUDINAL  U  IFFU3IQN  COEFFICIENT 
FOR  PARTICULAR  VALJES  OF  E/N  (GRID  VALUES  FUR  OLrtJ  ARE 
FROM  LOaXE  AND  PARKER) 

DIMENSION  UL 1U ( Jd  )  »  X I ( 3d  I 

CONNON/  CUN ST/ 3, NO, EPS  l , A_PHA ,G AMR  A , ECHO, MU , NUc , K , 6  FLO 
REAL  NOfrtU»nUt 

EXTERNAL  a,INTKP  _ 

OATA  FRO*  LORRE  AND  PARKER*PHYS  R£V,181,P  307  ANO  LOMKE  ANO  OaVI- 
,J  APPL  PHYS.48C 12),0EC  77*  P4995 

DATA  0LMU/.00o9242,  .0069242, .00 78  65 «•  008751  *.015234 «. 03 ID 7 » 

C. 061229  *  •  U901  , .207L39,.29o053,.l9269,.  11  i 9  56,  . 082  89 ,  .  08 00  8:>  , 

C. 084912  ,.  111456,.  171 903,.  257249,.  3618  3,.  484  85*.  571389*.  6*7*63, 

C. 725177,. 867117, .901086, 

Cl. 20595, 1.580 18, 1.95908, 2. 3 1774, 2. 75962 *3. 1288, 3. 4653, 3. 5777, 

C3. 6314, 3. 70162, 3. 93726,4. 24 172, 4. 0186/ 

DATA  XI  /3.254  3  7E-21,6.04  2  2  7£-2l,l.llQ96E-20,1.42495E-2i), 

C l. 87932E -20, 2, 20343E-20, 2. 5  39226-20, 2 .9 1608 E— 20, 3.2  4&38E-20, 

C 3. 7 1535 £-20, 4. 462 726-20, 5. 90609E- 20, 7. 8 4332 E-20, 9 . 587386-20, 
C1.61065E— 19*3.206276— 19»3.50556-19»2.19432E— 18*4.11 9G3E-18* 
C7.9195<»£-18, 1.023536-17,1. 50216-17, 1.653986 -17,2.26366-1 7, 

C2 .938336-17, 

C3. 2961c -17, 3. 85666-17, 4. 252 16-17, 4. 69796-17,5.48536-17, 6. 54 18 
C8. 80446-17,1. 15296-16,4. 09266-16,6.02146-16,9.54336-16,1.5816 
C2. 84716-15/ 

_XlN  -  ABSCXINI  _ _ . _ 

CALL  INTRPtXl  , OLRU , 38  ,2 , XI N , YOUT) 

□l.  is  dl/mue  Tines  hue  _ _ _ 

&L  •  A8SI YOUT *M(XlN J/ (XI 9*401 > 

RETURN 
END 


,UJ  Ui 
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SUBROUTINE  PEDERVCN,T,Y»  ?0«  NQO) 

OI MENS  I  UN  PD< NOO.NOO) ,YC M) 

EXTERNAL  OL,R 

COMMON/ CONST/S ,N0, EPS  I , A.PHA.G AMMA.ECHG, HU , HUE , K.EFLO 
REAL  MU «N0 .HUE 
C 

C  SION  CONVENTION  IS  THAT  THE  ORIFT  VELOCITY,  W,  IS  POSITIVE 

C_  WHEN  THE  ELECTRIC  FIELD. E.  IS  NEGATIVE 
SGN  »  -L.0 

IF  CYC4).LT.0.0)  SGN  -  1,0 

HUE  •  RCYC4l/NQ)/YC<il  .  " 

ALPHA  •  N0*1 . 76275E-17*E<PC“1.45027E-15*N0/A8SC YTNI I ) 

GAMMA  -  8.8IE-? 

POCL.U  -  0. 

P0C1.2)  -  ALPHA*ABSCmCYC4)/N0H  -  GAHMA*Y  (  3  I /V  (  4  I /MU 
PDC1.3)  -  — GAMMA*YC2 )/YC  4) /MU 

P0C1.4)  -  GAMMa*YC2)*YC3)/MU/YC4>**2  ♦  YI2I *ALPHA*MU6 

_  P0C1.5)  -  0.0  _  _ 

POC  2*1)  •  -l./OLCYC  4)/.N3) 

P0(  2.2)  «  SGN*rfCYC4)  /NO)  /  OH  ft  4)  /  NO  ) 

PD C  2, 3)  -  0.0 

P0C2.4)  -  YC2)*MUE/0LCYC4)/N0) 

POC  2.5)  -  0.0 
P0( 3. I)  -  0.0 

POC  3,2 )  •  ALPHA*ABSCwCYC 4) /NO) I  —  6AMHA*Y| 3 )/Y ( A ) /MU 

_ POC  3,3)  -  -  GAMHA*YC2)/YC  4)/NU 

POC  3. A)  -  GAMMA *YC2)*YC  3) /MU/ YC  4)  **2  ♦  YC2) *ALPHA4MUE 
POC  3.5)  -  0.0 

PD'  *,L )  -  0.0  .  ~  '  ""  “ . . .  ~ 

POx  4,2)  -  -ECMG/EPSI 
POC 4*3)  -  ECHG/EPSI/Y C4)/HU 
P0C4.4)  •  POC  4»2)*YC3I/NJ/YC4)**2 
POC  4,5)  -  0.0  ~  ' 

POC  5  •  I )  «  0.0 
'P0C5.2)  -  0.0 
POC  5,3)  -  0.0 
P0C5.4)*-  -L. 

POC  5,5 )  -  0.0 

C  PARTIAL  DERIVATIVES  ARE  NEGATED  TO  ALLOW  FOR  INTEGRATION 

C  FROM  SHEATH  BOUNDARY  BAC<  TO  CATHOOE 

POC 1,2)  ■  POC 1»2)*C— I)**<  v 

__  POC 1,3)  -  POC 1,3)*C-1)**< 

POC 1*4)  -  POC i.4)*C— !)••<  . 

POC  2, 1 )  •  P0C2*l)4C— l)*4R 
P0C2.2)  •  POC 2»2)*C-1 )*•<  * 

P0C3.2)  -  P0C3,2)«C-L)**< 

POC  3,3)  -  P0(3,JI*{-l)M(  . 

POC  3 ,4 )  •  PDC3,4)«C-L)**< 

P0C4.2)  -  P0C4,2)*C-1)**<  -  * - 

_  POC  4,3)  •  POC  4,3 )*C-l )**< 

POC  4,4)  -  P0C4,4)*C-1)**«  * 

POC  5,4)  -  P0C5,4)*C-1)**« 

RETURN 
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SU3RQUT  INE  INTRP 


PURPOSE 


c 

r 

INTERPULAT 

E  BETm 

EEN 

data 

POINTS 

TO 

c 

DESCRIPTION 

OF  PAR 

AMETERS 

c 

X 

array 

OF 

DATA 

POINTS 

FOR 

c 

Y 

ARRAY 

OF 

DATA 

POINTS 

FOR 

M T S  NUMBER  Of  P  4 1 R  S  UP 

N  TERM  S  NUMBER  Of  TERMS  IN 

XIN  INPUT  VALUE  OF  X 

TOUT  INTERPOLATED  VALUE 


EVALUATE  A  FUNCTION 


INDEPENDENT  VARIABLE 

dependent  variable 

DATA  POINTS 
FITTING  POLYNOMIAL 

OF  V 


SUBROUTINE  I NTRP I  X . Y . NP TS , N TER MS  *  XIN, TOUT ) 
DOUBLE  PRECISION  OELTAX. DELTA, A.PRQU.SUM 


of  xm 


DI MENS  I  UN 

X  (  1 )  •  Y  (  I 

) 

DIMENSION 

DELTA! 10 

) «  A( ID ) 

SEARCH 

FOR  APPR 

3PRIATE  VALUE 

11 

03 

19  I  - 

l.NPTS 

. . 

IF 

(XIN  - 

XU))  13 

•  17.  19 

13 

11 

-I  -  NTERMS/2 

IF 

(ID  15 

.15*21 

15 

11 

-  1 

. 

GO 

TO  21 

17 

TOUT  -  YC I 

) 

18 

GO 

TO  61 

19 

CONTINUE 

11  -  NP  TS 

-  NTERMS 

♦  1 

^1  12 
IF 
f2 
II 

~25~  IF 
26  II 


27 


-  II  ♦  NTERMS  -  l 
C  NP  T  S  -  12)  23.31.31 

-  NP  TS 

-  12  -  NTERMS  ♦  1 
CH)  26.26.31 

-  1 

NTERMS  -  12-  IV  ♦  1 


"EVALUATE  DEVIATIONSDELTA 


31 


- 35 


DENOM  -  X(  ID1I  -  XC II) 

DELTAX  -  (XIN  -  X(I1))/0=N0M 

03  35  I  -  1. NTERMS  . . . 

IX  -  II  ♦  I  -  I 

OElTA(I)  -  « X  C IX  t  -  XIIU)  /  QENOH 


4 


ACCUMULATE  COEFFICIENTS  A 


"All)"-  Yd  1) 

00  50  K-2,  NTERMS 
"PROO  -l,  ' 

SUM  -  0. 

IMAX  -  X  -  l  .  « . 

IXMAX  -  II  ♦  IMAX 
"DO  49  I  -  1,  IMAX " 

J  -  K  -  I 

PR30  •  PR  00* ( 0 tLT A ( X )  -  DELTA! jf) 


non1 


49  SUM  -  SUM  -  A  4  J )  /PROD 

50  44X1  •  SUM  ♦  Y  4  IXMAX )  /POO 

ACCUMULATE  SUM  3F  EXPANSION 

51  SUM  -  A  4 1 ) 

00  57  J  -  2*  NTERMS 
PROD  -  1. 

IMAX  -  J  -  l 
00  56  I  -  1*  IMAX 

56  PROD  -  PRU0MD6LTAX  -  0E-TA4I)) 

57  SUM  -  SUM  ♦  A4J)*PR3D 

60  TOUT  -  SUM 

61  RETURN 

ENO  . .  " 


SUBROUT INE  TRIMAT 


SUBROUTINE  TRIMAT 4  A f  N « NOyPRjSA VE 1 ) 

RETURN 

ENO 


SUBROUTINE  PEO 

< 


SUBROUTINE  PE04 N*T , Y.PW, NO, CON, IER I 

RETURN 

ENO 


no  loon 


subroutine  plokli 

01  MENS  I  UN  T(150)«J:(150),NE(150),NI(150),JP(16U),E 
CPT  (  160 )  ,  P  Jt  (  130)»PN£(150)*PJP(150),PE(  150) ,PNl t 150 
C1PAK( 150)  ,X(2)  ,2(2) 

COMMON/ PL T/T  ,  Jt  ,NE,NI,  JP,E,V,VMIN,M 
REAL  JE i JP»NE »NI 
CALL  COMPRS 

CALL  BGNPL(L)  _ 

CALL  3ASALF(6HL/CST0) 

CAuL  Ml XALF(3HSTAN0ARD) 

CALL  ri  TL€<"  S'*. -100, 

C"(3)ISTaNCE  -  CMS,,,100,,,(N)6»IN)*  -  1/CH3S"  ,  100, 8  . 
CALL  CKOS  S 

CALL  GRAF(0.0,.l,.33,-2.  =  8,l,68,2.E8)  _ fQ 

2(1)  -  2(2)  -  0.0 
X(l)  -  .30 
X(2)  -  .33 
CALL  CURVE  (X,2,2,0) 

CALL  MESSAG("(C)ArHOOESw,100,-.5,-.5) 

CALL  NESSAGC(A)N0  3Eiw»130»7.7,-.5) 

CALL  8LNKK2. 4,5.65,0.03,1. AO) 

FI.L  IPaK  ARRAY  KITh  LEGEND  INFORMATION 
I0JI1MY  -  L1NEST(  IPAK,  150,80) 

CALL  LlNES("(N)E  E-  NUM  DENSITY  (  X )  IOCS"  ,  I P  AK,  1 ) 
CALL  L I NE  S (" ( N ) ♦  IJN  NUMBER  OENS IT YS", IPAK, 2 ) 

CALL  LInES("(J)£  ELECTRON  FLUXS", IP AK , 3 ) 

CAlL  LINES("(J)»  ION  FLUXS",  IPAK,  4) 

CALL  LINES  (M(V)  VOLTAGES", IPAK ,5 ) 

call  lines  ("(E)  electric  fields", ipak, 6) 


(  150) 
)  ,PV( 


, V ( 153  )  , 
130)  • 


CALL 
2(1) 
X(l) 
X(  2  ) 
CALL 
CALL 
CALL 
CALL 
FI.L 


0,6.0] 


CALL  LINES  (M(V) 
CALL  LINES  ("(6) 


REVERSE  ANODE  ANO  CATHODE,  CHANGE  ZERO  IN  POTENTIAL  FROM  THE 

ANODE  TO  THE  CATHODE,  ANO  MAGNIFY  NE  BY  100 

00  350  MM  -  1,M 

MMN  -  1  -  MM 

PT(MM)  -  T(MMM) 

P  JE (MM)  -  JEfMHM)  _ 

PNE (MM)  -  NE ( MMM) *  100. 

PJP(MM)  -  JP(HMH) 

“PE (MM)  -  E(MMM) 

PV(MM)  -  V(MMH)  -  VMIN 
PNKMM)  -  NI(MMM) 

CONTINUE 

00  NUT  PRINT  ERRORS  IN  PLOT  VALUES  (OUT  OF  RANGE)  FOR 
PLOTS  OTHER  THAN  THE  FIRST 
IF  (L.GT.O)  CALL  N0CH6K 

CAlL  SPLINE  _ 

QALL  CURVE(PT,PNE,M,10) 

CALL  CURVE(PT,PNI,M,10) 

CALL  YSRAXS(-4.G£li,l.EU»*4,0Ell,6#0«”(  JIE.4J)*  -  1/CM2/SEC3 
C-100, 8. 0, 0.0) 

‘CALL  DASH  -  -  , . 

CALL  SPLINE 

CALL  CURVE(PT,PJE,M,10) 

CALL  CURVE(PT,PJP,N,10r 

CALL  YGRAXSI— 60,0*20, , 60 .0,6.0, "( V)OLT AGE  -  ( V ) S" ,-100, 8. 7,0, 
CALL  CHNOOT 
CA-L  SPLINE 


CALL  CJRVE1PT  »PV»H»LO) 

CAWL  YGRAXS(-20O.»i00.t20O.f6.t"(E)  IF)I6L0  -  I  V  ) /CM"  #  100^^.  7 , 0 . 

CALL  CrtND  SH 

CALL  SPLINE 

CALL  CURVEiPT ♦P6*M»10I 

CAlL  RE  SE  T ( "3LnK  L" ) 

CAlL  LEGENUt IPAK»6»2.i*0.20l 
CALL  6NOPL (L  ) 

Call  oqnepl  '  .  . ‘  "  '  '  ’ 

RETURN 

ENO  ""  .  ~  ■  ' .  "  . . . . 


m i>  - 


t 
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"PROGRAM  CaTeI  I  (  INPUT, OUTPUT, TAPE3  -  OUTPUT) 

DIMENSION  YO< 5> * T< 200 >• J5 ( 200) *NE( 200 >«NI! 200) «JP( 2001 • 

rr<  200  ),V  1200)  ' 

COMMON/ CONST/S.NO. EPS  I. ALPHA *GAMMA*ECHG«MU,MUE.IUEFLO*NP 

“Common/ plt/t.je»ne.nT  *jp,e.  v.Vm  in  ,diishth»oashth,  magmas  *  o' 

C0MN0N/GEAR9/  HU SEO* HOUSED  »NSTEP»NFE*NJE 
TxTer  na  l  '  0  IF  IHJn.PEDERV,  wTf  l  o 

REAL  ME»NO,NP,MU,MJE,MG,JE.NE,NI . JP ,MUCHEP , MGAM 


Y0(2I  -  SQRn8.0Ell/8.aiE-7)*.85 
STERNE  -  YO ( 2 ) * • 01 
20  CONTINUE 


EES  * I 
P33SQ3SI 

I  1  Vi  i 


N  MA55 
ME  -  .911E-27 


N 

BK  -  1.38E-16 


[;14ri44: 


INDEX  FOR  FIRST  CALL  TO  GEAR 


GAS  PRESSURE  ( 0YNE/CM2) 
>RE$  -  2V0.*(  1.013E6/760.  ! 
GAS  TEMPERATURE  <K) 


-  293.  • 

ELECTRON  TEMPERATURE  (X) 


112 


c 

mSlFCOlar  WEIGHT  OF  GAS  ! GM ) 

MG  -  40.*1.67E-27 

c 

POSITIVE  ION  MOBILITY  ! CN2/SEC/V) 

MO  -  4. 9E 19/NO 

c 

"RECOMBINATION  COEFFICIENT  ! CM3/SEC ) 

GAMMA  -  8.81E-7 

c 

IONIZATION  SOURCE  1  ION  P Al RS /CM3/SEC ) 

S  -  8. Ell 

c 

NET  CURRENT  IN  OISCHARGE  ! AMPS/CM2/SEC) 

C  -  .022E-6 

CC  -  C 

c 

ELECTRON  NUMBER  DENSITY  AT  START  ! 1/CM3) 

NUCHEP  -  MU*ECHG/EPS I 

MGAM  -  HUCHEP-GAMMA 

C - 

INITIAL  DISTANCE  STEP  SIZE 

TSTEP  •  .01*0 

NP  •  !MGAM*YO!2)M!HGAM«YO!2) l**2M.*MUCHEP*S)**.5)/!2.*MUCHEP) 
WRITE  ! 3* 1002)  YO ( 2 ) * NUCHEP * MGAM, NP 

c 

ION  CURRENT  OENSITY  AT  START  I1/CM2/SEC) 

EFLD  -  YO! 4)  -  -l.E-6 

Y0( 3)  -  0.0 

YOU)  -  C/ECHG 

C 

"POTENTIAL  ! VOLTS)  AT  START  SET  ARBITRARILY  TO  ZERO 

YO! 5)  -  0.0 

- c - 

"FIRST  OUTPUT  POINT  ufeslRED 

TOUT  -  TO  ♦  5.E-4 

wR I TEl TtTOOOl 

WRITE  13*1001) 

1009 

WRITE  ! 3*1009) 

FORMAT  !  13X••M*•i2X"ALPHA•llX"MUE•l^X»W•14X*,OL^NlOX,,OL/MUE,*/  ) 

WRITE  (3*1005) 

WRITE  ! 3* 1002)  T0*Y0!  1) *Y0! 31 *C*Y 0! 2>*NP*V0 !4) * YO ! 5 ) 

T - STDR1~™XMETER ’VAUJETKT  ftOUTOWf 


355  roio  -  YO!l) 

- Y020  -  Y0«zr 

Y030  »  YOI3I 


t  il  I  Cl  1 


Y050  -  Y0I5I 


FSS 

Ll'TYi 


M  IS  THE  NUMBER  OF  THE  INTEGRATION  POINT 


n  -  u 

SET  INITIAL  VALUES  OF  ME  AMO  J« 


Y02  - 


AR 

MIN  -  1000. 


SCN  -  1. 


1000  FORMAT!  1H1*3X*»0ISTANCE«10X,»JE*13X«4^*7X"CURRENT  DSNSITV"6X» 
t"NE"l'iX"M*"llX"fe  ElELD»ax"V0LTACE"4Xb0RbER,*l 

1001  FORMAT!  6X"CH»10X,tl/CM2/SEC',6X»l/CH2/SEC«6X,'AMPS/CN2"9X* 

C1»i/CM3Ml6xl,l  /EmT**  X" VOLt  S / C  M  Y  OlT  S"  /  /  /  » 


TOUT  -  TOUT  ♦  TSTEP 

rtnre 

390  CONTINUE 


I 

Mi 


DASHTH  -  TOUT 


400  CONTINUE 

- OKS  TOUT - 

C  NK  IS  THE  INDEX  NUMBER  OF  THE  CATHODE 


WRITE  (3*1006) 


Y0I4)  -  -EFLO 


r/tiL'ii 


ii  iVfii 


NP  -  Y060 


Y0( 1)  -  Y0(3)  ♦  C/ECHG 


EPS  -  .001 


TOUT  -  TO  ♦  5.E-4 


imitx  •  i 
WRITE  (3*1000) 


500  WRITE  (  3*1007) 

[TOTtor wat~ttx  "anw — 

00  »  OASHTH  ♦  OKSHTH 


IF  (L.CT.10)  GO  TO  300 
IF  t(bELO*( 00-0). LT. 0 . ( 


io  TO  70< 


i 
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H«ITE<3«100ai  H« ALPHA, HUE, QRIFTtOIFf.DHU 


.TOI<T»roill»r<>^»»c.yot?».w.>o«4i.yoisi.i.QU»n 


IF  C INDEX. EQ. 01  GO  TO  103 


-°°3  F0Rh.4T<//26H  eRR0R  RETURN  WITH  INDEX  ■•13//) 


VU  IU  c  \JsJ 

100  IFIK.EQ.l)  GO  TO  350 


350  CONTINUE 


IF  CC2.»Y0(2)~r02).LT.0.Q)  TSTEP  •  TSTEP/10. 


ro2  -  ro<2i 


/*>! 


I'M 


305  CONTINUE 


rrii 


m 

I 

■ 

RTifia 


ro4  -  roi*) 


TOUT  -  TOUT  ♦  TSTEP 


OtLO  •OO-I) - 

V0I2J  -  Y020  -  STEPNE 

TTO  TO  20  - — 

700  STEPNE  -  A8S ( STEPNE/2 . ) 

NEGFLG  -  1  - 

_ IF  <DELD*C00-0).LT.0.0»  SON  -  -SGN 

Feld  -  dd~=~o  - — - 

_  TO 12)  -  V020  -  STEPNE*SGM 

fP  ( A8^S(00-0r.LE. lT£-3*0)  GO  TO  900 
IF  <ABS(STEPNE).LT.ABS(.005*VO(2m  GO  TO  900 

gg— Yg  20 

900  CONTINUE 

If  TO PLT.EQ.1.0)  CALL  PLOT(MK) 

1004  FORMAT < //21H  PROBLEM  COHPLETEO  IN,I5,6H  STEPS/ 

~C  2ix* is* 14H  f  Evaluations/ 

_ C  2 IX , 1 5 * 14H  J  EVALUATIONS///! 

200  WRITE  I3tl00¥)  NSTEP* NFEtNJE 
300  CONTINUE 


END 


OUT INEOIFFUNI N,T,Y, YOOT ) 
DIMENSION  Y { N  )  , YOOT INI 


XTERNAL  W,OL 

COMMON/ CONST/ S, NO, EPS  1, ALPHA, GAMMA, ECHG, MU, HUE, K.EFLD.NP 
REAL  mu.no, mue 

CHECK  ON  MAGNITUDE  OF  EXPONENT  TO  PREVENT  UNDERFLOW  OF  RESULT 
rF  1 1. A8E-15*N0/ABS(Y(4)Y.GT.250.*K)  GO  TO  10 
ALPHA  -  N0*2.9£-17*EXP(-l . 48E-15*N0/ A8S (Y  ( 4 ) )  ) 

GO  TO  20 
10  ALPHA  -  0.0 
2HHC0 NT  IN U E 


CONVENTION  IS  THAT  THE  DRIFT  VELOCITY,  U,  IS  POSITIVE 
WHEN  THE  ELECTRIC  FIELD,  6,  IS  NEGATIVE 


IF  <Y(4).LT.0.0)  SGN  -1.3 

YOOT n I  -  S  ♦  (AL?HA*A¥S(  W<YH  )/N0)  )-GAMNA*YI3)/YI*)/HU>*Yi2) 
YOOT (  2  )  -  (Y<2)*W(Y(4)/N0)/SGN  -  Y( II ) /OL ( Y C 41/ND) 

YOOT  f  3  )  -  S  ♦  TaT7»HA*ABS  (WlY(4)/N0))  -GAMHAPY  C  3  ) /Y  C  4  f/HU  I  *Y  C  2  > 
YOOT <41  -  (ECHG/EPSI )*IYI3)/Y<4)/MU  -  YI2)> 


IF((T.EQ.O.O).AND.(ABS(Y( 2 ) — Y<  3 )  /  Y(  4) /MU) .GT.l. E~6*Y(2 ) )  )  GO  TO 


THIS  IF  CHECK  PREVENTS  ROUNOOFF  ERROR  FROM  MAKING  OJE/DX  AND 


Y00TC3)  -  YOOT Ill 


■il'HITB 

rmi 

EEB 


IF  (IT.EO.O.O) • ANO. (K.E0.2) )  60  TO  90 


50  YOOT  <11  -  -1. E-4*SQRT(S/ GAMMA )*Wf  EF.LO/MOI /E  FLO 


Y00TI3)  -  YOOT ( 1 ) 


90  CONTINUE 


9  ^ 

OIFF  EONS  ARE  NEGATED  TO  ALLOW  FOR  INTC6RAT ION  FROM 
ShTaTh  BTaCk  fCTCATMOoE 
TOOTH)  -  YOOT  1 1  )*I“l  •  )P*K 


-  YUOT 1 2 1-1*1 • )**X 
YOOT C  3)  -  Y00TC3)*C-1.)**K 
yoDTm)  -  yOoti  4)*  f — 1.  f*~*i 


FUNCTION  FLOtME) 


DIMENSION  WI ( 33 ) «Xl ( 33) 

COMMON/ CONST/ S, NO, EPS  I, ALPHA, GANNA, ECHC»HU,MUE,K,EFLO,NP 
REAL  NO  «  MU*  MUE  * NP 

DATA  *1/3169.2, 48 94. 3  , 76 5 f . 9 , 1 3002. • 222 93. » 3 02 09. , 5 369 1775 369 L. • 
C72756., 90267,, 104566., 113097,, 118391,, 13317 1., 14833 4., 173 524., 
C195821., 222439. ,250207. .291743., 343518., 390211., 45054777544565., 
C698072., 1029840., 1350700., 1731450., 2323430. ,3037390. , 3970640. • 
C56 14230. ,652484077 

OATA  XI /5.3765E-2 1,9. 1264E-21,!. 4096-20, 2. l543E-20,3.0262E-20, 
C3.3814E-20, 4. 26386-20,4.26386-20, 4. 81176-20, 6. 06746-20,8.43726-2 
2 126E-19, 1.75956-19, 2. 6122E-19, 4. 3053E-19 ,8.66166-19, 1.46056-1 
4T916— 18, 4. 40566-18,8. 2471E— 18, 1.64 82E— 17,2. 9868E— 17, 4 .84286— l 
1 6746-17, 7. 75186-17, 1. 1396-16, 1.5592E-16, 2. 02496-16, 2. 89076-16 
9555E-16,5.5182£-16,8.4372E-16,^.935'7E-16/' 

( ME. EQ .0.0 )  GO  TO  10 

-1.0  ~  " 

IF  IME.GT.O.O)  SGN  -  -1.0 

.  ABS(WE)  ■ 

CALL  INTRPI W1,XI,33,2,NE,X0UT) 

FLO  -  A8TrXOUt*NO)*SGN  ' 

GO  TO  20 

HTTLO  -  0.0 

20  CONTINUE 


~FUN<rr  I  ON  0 L  (  X I  Nl  ! 

_THI$  FUNCTION  CALCULATES  THE  LONGITUDINAL  DIFFUSION  COEFFICIENT 

for  particular  Values  of  e/n  cgrio  values  for  dlnu  are 

FROM  LOWKE  ANO  PARKER  I 
01 NENS 1 ON  DLMUC 38 1 «  X I  ( 38  I 

CON NON/ CO NST/ S, NO, EPS I, ALPHA* GANN AtECHG*MU*MUE*K»EFLD*NP 


EXTERNAL  W* INTRP 
REAL  NO  *HU*MUE*NP 


OATA  FROM  LOmKE  ANO  PARKER, PHYS  REV,161*P  307  ANO  LONKE  ANO  OAVIE 
*J  APPL  PHYS  *48(12)*0EC  77*  P4995 

OATA  OLHU/V 0069242 *• 0069242* .007865* . 008751  *.015234*. 03107* 

C. 061224  *.  12401*  .207139*. 2 *6053*  •  19269*.  1119  56*. 0«'>  99*.  08008$*  j 
C.084912 ». Ill 4  56, • 17 1403*. 267249*. 36183* .48485* .57x389*. 647963*  ? 

C. 725177 *.867117*. 901086* 

Cl. 20595  *1 .58018* 1.95908*  2. 31774*2.75962,3. 1288*3.4653*3.5777, 

C3. 6314, 3. 70162* 3. 93726* 4. 24 172* 4. 6186/ 

~1» A  fm7  3 . 2'6 4  3  7E^2 1 , 6 . 0 4 2  2  7 E -2 1 ,  l .  1 1 0 9 6E -2  0  *  1 .  42  49 5  £-2  0  * 

Cl. 87932E-20*  2 .20343E— 20*2 . 53922 E-20, 2 .91608  E— 20, 3.2  4638E-20* 
"C3U71535?— 20,  4. 46 27 2 £-20*  5 .90669E -20* 7.84332 E— 20*9. 58738E-20* 
C1.61065E- 19*  3.2062 7£— 19*  8. 5055E— 19*2. 19432E— 18*4. 11908c— 18* 
T779l954E-T8,1.0235TE-l7*1.502l6-i>*l.65348E“l7,2.2636E-l7, 
C2.93833E-17, 

~CXr296lE'-l7,3.d5666-iT,4.i521E-l7*4.69>9E-17,5.4853E“17,6.5418E-l7 
C8.8044E— 17* 1. 1529E-16*4.3926E-16«6.0214E-16,9.5433E-16* 1. 58166-15 


I  / .  N«  r  \  F*  I  s  # 


IF  (XIN.EQ.O.O)  GO  TO  10 
■~~AB'S  f  X  l  N 
CALL  INTRP<XI,0LNU*38,2,XtN,Y0UT> 


GO  TO  20 


ED 


20  CONTINUE 
1 
ENO 


FUNCTION  W(XIN) 

C  THIS  FUNCTION  CALCULATES  ELECTRON  DRIFT  VELOCITY  (CM/SEC)  FOR 

Z  Particular  Values  of  e/n  (grid  values  for  drift  velocity  are 

C  FROM  ENCELHART  AND  PHELPS)  _ . 

DIMENSION  Mll  (  33  )  •  X I  ( 3 3  ) 

_ COMMON/CO N S  TV S ♦NO, EPS  I , ALPHA , GAMMA, ECHG, MU, MUE,K,EFLO,NP _ 

REAL  n6,MU, HUE ♦ NP 

C  OATA  FROM  ENGELHARDT  AND  PHELPStPHYS  REV,( 133) ♦ JAN  64.PA377 

VST  A  W I  / 3 169 . 2  ♦  48  94 . 3  ♦  765  7 . 9  ♦  1 3  0  0  2.  ,  2 2293  .  ,  30209.  •  5  369 1 .  ♦  5 369 1  •  ♦ 

C 72 756., 90267. ,104566. ♦ 113097. ♦ 1 18 391. ,13317 1. , 148  33 4. , 173524. • 

C 195821. ,2  2 2  43 9. , 2 50  20  7. , 2 9 1 743. ,343518. ,390  211. ,450547. ,544565. , 
C698072. . 1029840.,! 350700. , 1731450., 2323430. ,3037390., 3970640. , 
C561 42 30. ,6524840./ 

DATA  XI/5.3765E-21,9.1264E-21«1.409E-20,2.1543E-20,3.0262E-20, 

£3. 3814E— 20, 4. 2638E-20, 4. 2638E— 20, 4.8117E-20 ,6.0674E— 20, 8. 4372 E”20 
C 1. 2126E- 19 »1.7595E-1 9,2.61226—19,  4.30536— 19, 8.66 16E— 19,  1.4605E— 1'> 
C27T7V'lE-18,4.4056E-1878.2  47lE-t8,  1.  64026-17 ,2.98666-17,4.84286-1 
C6.1674E— 17,7.75186-17,1. 1 39 E-l 6,1.55926—16, 2.02496—16,2. 8 907E— 16 , 
£37^5556^16, 5. 5182E-16, 8. 4372E-16 , 9.9357E-16 / 

IF  (XIN.EQ.O.O)  GO  TO  10 
501*  -  ABS  ( XI N) 

CALL  INTRP(XI*MI ,33, 2, XI N, TOUT) 

- ins-YOUT - - 

GO  TO  20 

u  „  n  n  -  ■  ' 


20  CONTINUE 

- RETURmT 

END 


“SUBROUTINE  PEDE RVlUiT*  Y,  P  6  , NOO I - - 

JJINENSION  PO( NOO* NOO ) ,  Y(  N ) 

EXTERNAL  OL,W  — - 

_A9**9N^CONST^  S*N°*EPSl,ALPHA*GAHRA*ECHG,HU,nUE*K,EFLD,NP 
REAL  N0,nU,NUE,NP~“  "  “ - — 

SIGN  CONVENTION  IS  THAT  THE  DRIFT  VELOCITY*  W«  IS  POSITIVE 
WHEN  THE  ELECTRIC  FIELD,  E,  IS  NEGATIVE 

sGn  — i.o  - - 

_IF _ CYMI.LT. O.O)  SGN  -1.3 

ff't  41.  EQ«  0. 0 )  Y  (4  )  -  -I.E-10 - — ~ - 

"7u>6  -  W(Y(4>/N0)/Y(4) - 

CHECK  ON  MAGNITUDE  OF  EXPONENT  TO  PREVENT  UNDERFLOW  OF  RESULT 
IF  ( 1. 48E—15*N0/ ABS ( Y ( 4) ).GT.250.*K)  GO  TO- 10 
ALPHA  -  NQ*2.9E-17*EXP(—1.48E— 15*NQ/ABS CYC4D) 


ALPHA  - 


ONTINU 
PD< 1,1)  -  0. 

P0TL,2I  -  ALPHA*A8S(W(  YC41/N0H  -  GANHA*Y<  3  i /V  (  4  l/HU 
POC 1,31  -  — GANN A*Y 12) /Y  <41/ NU 


l? 

PD( 1* 5 )  -  0.0 


P04  3,1)  - 


■— 


PDf 3,5)  -  0.0 


U> 


POC 4,4 )  -  POC  4,2) *Y<3 l/NJ/Y(4l**2 


rDf  4*4 J  -  rOl 4 1 4 )  * (— l.)**K 
P0<5*4)  -  P0(5*4)*(— 1«>**K 
RETURN 
6  NO 


SUBROUTINE  >E0«t«,r,y7FM,N6.tdN*ttftJ 

RETURN 


I 


SUBROUTINE  PlotTmIO 

_  01  MENS  I  ON  T(200r.J£(200),NE(200).Nl(2qO).JP(200).£(200)  *Y(200JL* 

CPT( 200) »PJE( 200 ) ,PNE ( 200 ) , P JP( 200 )  ,PE ( 200) , PNI ( 200 ) ,PV( 2001 , 
CIPAK( 150) *X(2),Y(2) 

COMMON/ CONST /S.NO, EPS  I, ALPHA, GAMM A, ECHG, MU, HUE, K.EFLO.NP 
COMMON/ PLT/T.JE, NE.NI.JP ,E , V, VM I N .OKSHTH, OA SHTH.MA , HAS , 0 


REAL  JE , JP.NE.NI ,NP 
H  ■  MA 
CALL  COM  PR'S 
CALL  BGNPL(MK) 

CALL  8T$ALF(6HL7CSTD) 

CALL  Ml  XALF(BHSTANDARD)  /*>• 

C'AinTTLEC 
C.-100, 

t"(  D  ) ISTANCE  -  Cn$”,100,*(N)E,(N)*  -  1/CM3S" . 100,8.0,6.0) 
CALL  CROSS 

C ALL  GFaFTO .0,. 1,. 33,-8. E 8,2 .£8, 8 . 68 > 

X( 1)  -  .3 


X(2)  -  .33 
YCl)  -  0.0 
Y (2)-  0.0 
CALL  CURVE(X,Y,2,0) 

- CalL“mESSag(-(C  1  ATMOOEi*  ,  100',-. '5, -.5  * - 

CALL  ME  SS AG  (**(  A  ) MOOES*,  100, 7.7,— •  5) 

CAUnOoTli  2.4, 5. 65,0. 05,1.26) 

CALL  LINES("(N)E  E-  NUMBER  DENS XTYS*, IP AK ,1 ) 


2 

CALL  LI N£S("( J) £  ELECTRON  FLUXS" , IPAK ,3 > 

- call~ltnes<  - 

CALL  LINES  <"<V)  VOLTAGES", I PAX, 5) 

- ULL~urNEs(wm  ~a'EcTRTrTaEro'fw,  ifak  ,6> - 

00  350  MM  -  1 ,MK 

- REVfft SE  POSI T IONS_07“‘8oU'M0 ART^AMo'  CATHODE,'  SHlET  ZEftO  "IN - - - 

POTENTIAL  FROM  BOUNOARY  TO  CATHODE,  AND  CHANGE  E  FIELD 

- DIMENSION  FROM  V /CM  TCTXvVCM 

MMM  -  MK^l  -  MM 
FTTHMV  •  Dk'SH T  H  -  Y  (  M K Ml 

P JE (MM)  -  JE(MMM) 

PNETmMI  ■  N£(  MMMl 
"  PNI(MM)  -  NI(MMM) 

~  FJPfNMT~-  JPYMMM)  ‘ 

PE(MM)  -  E(MMH) 

PVTMH)  -  ST  (MMM  i~t«Tn 
350  CONTINUE 

r~MKi  -  MK~V  1  T 

'  DO  450  MM  -  MK1.HA 

- FT  (  MM )  TTM  M  >  ♦  D  "-"DAS  NTH - - ! - r 

-  PJE(MM)  -  JE(MM) 

‘  FNETMM)  -  NE  ( M'Ml  i 

‘  PNI(MM)  -  NI(MM) 

PJP(MM)  •  JP'(MM)  I 

PE(MM)  -  E(MM) 

POTEnTTAL  AT  ANOoe  ShEaIM  BOUNDARY  IS  SHIFTED  BY' CATHODE 
FALL  (VMIN)  AND  BY  THE  VOLTAGE  DROP  ACROSS  THE  POSITIVE 
COL U MN  TwHIC H  IS  LINEAR  IN  X) 

PV(MM)  -  V( MM )  -  VMIN 

Continue  ' 


t 


1 3  « 101  HK 

WRITE!  3*  10)  HAS _ _ _ 

WRITE (3* 10)  HA 

10  FORHATI 1X*I5*7IIPE13.5*2X) )  _ _ _ 

00  20  HH  -  1 « HA 

WRITF(3.10)HM.PT<HM).PJEIHH)«PN£I  HH) »PMI(HH  )»PJPIHH)»P£IHH)»PVIH,- 

20  ConTTnUE 

CALL  SPLINE  _ 

Call  curve  (Tpt  ,  pnETh7T67  ~ 

CALL  CURVE«PT*PNI,H»10)  _ 

CALL  YGRA~XSC-8.0Ell*2.Ell*»8.0Eli.6.0*"<J)E  •!.!)♦  -  1/CH2/SEC*** 
C-100,8. 0.0-0)  _  ' _ 

Call  dash 

CALL  SPLINE  _  _  _ _ _ 

CaXlT CURVETpT.PJE*H,10) 

CALL  CURVE!PTtPJP*H»lO)  __  _ 

CAUL  YCRAXS (-8.0,2. ,8. 0,8.0. ■(VIOLTAGE  -  IV ) $".-100, 8. 7,0.01 

CALL  CHNOOT _ _ _ _ 

CalITTPlTnE 


CALL  CURVElPT.PV.H.lO) _ _ _ 

Call  yg r ax S( -T6o77*o. ,i60.*6.f"(E)  IF)  ielo  -  (VF7cTn1SToro7^777oT ) 

CALL  CHNOSH _ _ _ _ 

CATTCSPLInE - “  ~ 

CALL  CURVE(PT,PE,H,10I  _ _ _ 

rALL"REsn,{irffLNi?r*a  “ 

call  LEGENOI  IPAK.6,2.5,0.20) _ _ _ 

CALL- ENDPLiMiO 

CALL  OONEPL  • _ _ _ 

TfETWN - -  ~ 


ENO 


DIMENSION  YO(‘j|  «T(AbO)tJE  t  ASQ  >  *  NE  I  A  50  )  «  N  If  *  *0  I  *  J?  C  450 1  * 

UJf  30)  *V  L*J>QlI _ 

COMMON /CONST/S. NO* EPS  I,  AL PHA,G A MMA* ECHG *  MU* HUE ♦ < ♦  EFLD 
COMMON/ PL T/T  ,  Jt  tNE  .  Nli  JP  ,F*V*VHlNtQ<<SHTH.;)ftSHTH*nA.HAS* 
C0MM0N/GEAR9/  HUSED* MOUSE  D*NS  TE  P »NF E  »  N  J  E 


REAL  ME  *  NO«NP  «  MU  *M J  £  «MG  *  J  E  *Nfc  »N  l »  JP 


3ISSPLA  PLOT  OPTION  (l-YES  O-NO) 


SET  CONSTANTS  FOR  GEAR  PROGRAM 


NUMBER  OF  0  IFF  EOUNS  IN  SYSTEM 

N  «  S 


LOCAL  ERROR  TOLERANCE  PARAMETER 


L*I*J 


STEP  SIZE  IN  GEAR 


GEAR  METHOD 


EKICK  «  0.0 


ECHG  »  1.6E-19 


911E-Z7 


Pi 


1.306-16 


GJPJEl 


INDEX  FOR  FIRST  CALL  TO  GEAR 
N 


GAS  PRESSURE  IOYNE/CN it 


I A  :l  AW4IHU 

mmmmm 


GAS  TEMPERATURE  (K! 


L25 


C  -  <*8.E-3 


ELECTRON  NUMBER  06NSITY  AT  START  (1/CM3) 

AMMA 


POSITIVE  ION  NUMBER  OENSITY  AT  START  (1/CN3) 


mam 


ESTIMATE  Of  ELECTRON  CURRENT  DENSITY 
YOf 1 )  - 


00  12  JJ  - 


IF  (ASS(YU(1>  -  CE).LT.1.E-2Q*Y0(1I)  CO  TO  15 


k««l 


WRITE  (3,999)  Y0(1),Y0(2) ,EFL0, GAMMA 

—9.99  P QAM  AT(1QX,5( 1PE2Q.12)) _ 

12  CONTINUE 
_ 15  CONTIN 


Y0( 1)  -  CE 

: _ FROM  JE  ANO  NE , 


WRITE  (  3,1002)  Y0(  1 )  ,Y0(  2  ) 

WJE— »  .YQ-C 1  LLtJlU.  J _ 

THUS,  ESTIMATE  E  FIELD 


MUE  *  WE/tFLO 


V0( 3  I  *  MU*EFLO*NP 

- RE-RTUR6_^FJEL0  TJ  BEGIN  INTEGRATION  TOWARD  CATHOOE _ 

Yom  -  6FL0  ♦  .OL*fcFLO 

- ?. 0.1XNTIAL.  (VOLTS)  AT  START  SET  ARBITRARILY  TO  ZERO _ 

Y 0 ( 5 )  -  0.0 

- wRJ  Ii-iJLtilil  y^IUa.YD.I  2JjiiajJLUI0.(  AlxlSUll _ 

WRITE  (3,9991  NP,MUfc»MU,WE,EFLO 

_ f XR S T_.gu.TPUl  POINT  DESIRED _ _ 

TOUT  •  TO  ♦  1.  E-4* 

_ LL  L_UjlA915J _ ; _ 

WRITE  (3,1001) 

_ wfiJJL.13LdLQim _ 

1009  FORMAT  (  13X*»M,,12X*,ALPHAH1  lK“«UE**l3X"W“iA)<-DL"iOX''DL/nOE"/» 

_ w  R  I  TE_  J  3  ,_1.QJ> 5J _ _ _ __ _ 

WRITE  (3,1002)  TO,YO(l)*YO(3),C«YO(2 ) . HP , YO ( A) , f 0 ( 5) 


355  Y010  -  YO(l) 


Y030  •  Y0( 3) 


26 


*  IS  THE  NUMBER  OF  THE  INTEGRATION  POINT 
1  » 


SET  COMPARISON  VALUE  FOR  VOLTAGE 
VM LIN  » 


SGN  -  1.0 


FORNATI  1H1»3X**DIS  T  ANCE"  10  X*ME«l  3X*Ja"7X«CURREnT  OEMSI  TV»6X. 
•,NE"13X"N*"11X"£ 


1001  F  QRrtATC  bX^CM^lQ  X**l/  CH<i/SE  C**6X"l  /CN2/SEC**0X'*AHPS/CN2*,9X* 


1005  FORMAT  Iix-CATHOOE  SHEATH  BOUNDARY" 1 


AO  CONTINUE 


BEGIN  LOOP  TO  SOLVE  DE  SYSTEM  AT  VARIOUS  POINTS 


C  -  ABSIYOd)  -  Y0( 3 ) )*ECHG 
NP  »  Y0< 31/MU/YOIa 1 


ORIFT  -  tf ( YO ( 4 ) /NO ) 


i  XRIT6  (3*1002)  TOUT*  YOU).  Y0<  3 ).  C.Y3U  ),  NP.  YG( 4 1 .  Y0<  51  .N0US6  0 

-L.Q-Q.2_F  QAM  AT  ( 8.U.PE12  »4.3X)«i4)  '  , _ 

IF  (INDEX. EQ.O)  GU  TO  100 

- *ALTE  L3j  1Q03)  INDEX  ' _ _  _ 

1003  F ORMAT ( /  /£6H  ERROR  RETURN  MITH  INOEX  -*I3//> 


9B  EKICKU  »  EKICK 


37*  CONTINUE 
roi  *  vodi 


Y07  •  YG(2» 
Y03  «  YOI3) 


YO*  •  YO  <  *  I 


TOUT  -  TOUT  ♦  TSTEP 


375  CONTINUE 


INTEGRATE  ONLY  PART  WAY  THROUGH  POSITIVE  COLUMN,  THEN 


IF  C  TOUT »G£. I •8*0-DKSHTNI I ?GO  TO  *10 
UT  -  " 


GO  TO  *0  , 

>5  CON T I N 


TIM)  -  TOUT 


PI4UJ 


NEIMI  ■  YOU* 
_N  I  (  Ml  «  NP 


128. 


IF  ITtM>.LT.S.*TSTEP>  CD  TO  3B9 

J F _ .( XQ12  I .  Cl .JCA2 )  GO  m  3a 8 _ 

IF  <Y0<3)/MU/Y0<4).GT.Y03/MU/V04)  CO  TO  98 

_IF  (YO(2).GT.l.E3Q)  CO  TO  388 _ 

CO  TO  389 

EK  irxft  *  E  K 


GO  TO  99 


Y04  -  Y014) 

YQS  ■  Y  Q ( b ) _ 

TOUT  -  TOUT  ♦  TSTE? 


I1IJKTI] 


390  CONTINUE 


OASHTH  -  TOUT 
CO  TO  5 


400  CONTINUE 
KSHTH  - 


MK  IS  THE  INDEX  NUMBER  OF  THE  CATHODE 


mini 


WRITE  13*100 81 
MAT  (  1X“CATH 


SET  K  •  2  FOR  INTEGRATION  TOWARD  ANOOE  SHEATH  BOUNDARY 


RECALL  PARAMETERS  AT  PLASMA/SHE  ATH  BOUNDARY 


. 

Y0I2I 

_ _ 10121 

«  Y  020 
»  Y030 

YO  (  4) 

»  EFLO  ♦  l.E-8 

! 

• 

■■Hirni 

f 

HP  » 

YO 1 3 ) / MU/E  FLO 

; 

r 

_ £  -  CC _ 

K  • 

►  » 

HO  • 

.00001 

WRITE  13*1001) 


ffnjrjnEifTjHl 

■CThI  viTwitiTi  I 1  nT  zanTl  O  fvTT  U 1 


GO  TO  30 

ULcomny. 

M  •  M  ♦  1 


JEIM)  «  YO(l)  *  C/cCHC  *  YOt  3) 


WE  «  W ( YO ( 4 )  /NO  ) 

a  ■  YO(2i» wl _ _ 

WRITE  (3*10121  YOt  L)*Y0(2  )  *Y0< 3 ) «Y0( 4) * Y0< 5 ) • WE* B 

JLGKJIAI _ LI  Xj  7_(_Ui.EJLi»_Z.Ll _ 

TIM)  ASSIGNED  A  VALUE  OF  ZERO  AT  MAS  SINCE  ACTUAL  POSITION  OF 


T(HI  -  TO  - 


T STEP  -  i.E-3 


EPS  -  .001 
runt*  -  t 


WRITE  (3.1000) 


IJIMMIJlOi 

Ktnaw 


WRITE  (3*1011) 
nRMAT  tlx" 


WRITE  (3.1002)  TO*YO(1)»YO(3I*C*YO(2)*NP*YQ(4)»YQ<5) 

cn  rn  30 


500  WRITE  (3*1007) 


IF  (OPLT.EO. 1.0)  CALL  PLOT(MK) 


2 IX  *  1 5  *  14M  F  EVALUATIONS/ 


M 


200  W 

l  _ 300  C 

RITE 

ONTI 

(3*1005)  NSTEP.NFE.NJE 

NUE  '  ...  ,  'v 

i\ , 

CA 

_ EM. 

LL 

Q _ 

EXIT 

0  I  MENS  ION  Y<N» * YOOT <NI 

■lAT6jHAL.Wj.0i _ 

COMMON/CONST/S.NO.EPSI .ALPHA, GAMMA, ECHG, MU, MUE.K, EFLO 


CHECK  ON  MAGNITU06  OF  EXPONENT  TO  PREVENT  UNDERFLOW  OF  RESULT 


nrvi'i 


ALPHA  »  NO*2.9E-17*EXP(-1-.48E-15*NO/A9SIY(A))) 


11* 


SIGN  CONVENTION  IS  THAT  THE  OR  I  FT  VELOCITY.  W.  IS  POSITIVE 


SGN  — 


YOOT  ( 1 )  -  S  ♦  <  ALP  H  A*A  8S  (W(YI4I/N0I)  — G  AMMA*Y(3)/Y(4)/MUI*YUI 


klilU 


wvarfvv 


YOOT ( 3 )  -  S  ♦  ( ALP H A*ABS IWIYIAI/NOII — GAMMA  *Y (3)/Y(4)/MUI*YU) 


PM 


LM.ir 


YOOT 15)  «  -Y(4l 


THIS  IF  CHECK  PREVENTS  ROUNDOFF  ERROR  FROM  MAKING  DJE/DX  AND 


[I  vt  t  •  fTta  «  rw  •  «>li 


IF  ( (T.EO.O.O). AN0.IK.E0.1>>  GO  TO  30 
GO  TO  40 


30  YOOTUl  - 


liffi 

nlSnfnBSn!] 


YOOT 13)  -  0.0 


■mvmmH 
r  f  I'fflr  ,  .Y^1  J  * 

H*mi*m3*i 


YOOT <  5  I  -  -EFLO 


IF  (K.E0.21  GO  TO  50 


SPATIAL  DERIVATIVES  OF  JE»  HE.  J+t  AND  E  SET  TO  ZERO  IN 


50  YOOT (II  ■  YOOT 12)  «  YOOT <  3 )  -  YOOT I A )  - 


l*l*l 


60  CONTINUE 


GO  TO  80 


riiM*C*r»ilaTI: 


M3* 


YOOTUl  «  ♦1.E8 
Y0QT(.2>  -  -1.E8  _ 

*00T<3>  -  YOOT (II 
TI4)  -  ♦ 


YOOT (51  -  -EFLO 


SHEATH  BACK  TO  CATHOOE 


U*I»il 


11*1*1 


YOOTUI  -  Y00TUI«(-1.I**K 


YOOTUl  -  YOOTf  4IM-i.l**K 


kl*l*ll 


FROM  ENG6 LHAR  T  ANO  PHELPS  I 

JLLMEiiSlO  N-_aLU2  JUXLU^J _ _ 

COMMON/ CONST/  S  , NO,  EPS  I « ALPHA *2 AMMA ,ECH G, MU ,MUE ,K ,E FLO 


DATA  FROM  6NGELHAROT  AND  PHELPS, PHYS  RE  V, I 1 33 » *  JAN  64,PA377 

•  V  >  I  :  •  »  •%  •  A  -%  /  ,i  .\  J  -»  t  I  r  r»  /•«  1  rt  /%•»  "l  A  «  n  A  r  A  /  A  «  C  ^  i.  m 


C 72756. ,90267. ,104566. ,113  097.  , 1 18  39  1 .  •  1  331  71 .  «  1*8  3  3  4.  ,  1  7352  4., 

mean  1  i  /  1  >1  1  ^  A  1  A  1  -»  r%  1  *¥  ■>  T/.  1  C  1  U  lOrt  1  1  1  A  C  nc  A  1  C  A  A  £  aC 


C6980  72.  . 102 93 AO. « 13  30  700.  ,173 1* 50. *2323430. « 3037390.*  3970640.* 

5614730. .6624840./ 


OAT  A  XI  /6.3765E-21,  9. 126  A  E-2 1, l . 409E-2 0 , 2 . 1 54 3E-2 0, 3.0262E-20 

•t  -*  >  1  ,  r  t  A  i  ,  a  f  -»  A  1  *  <  A  r.  r*  n  A  s  a  t  i  V  r  ^  A  a  >\  /  1  i.  f  "*A  A*  i.  *»  "* 


C1.2126E-19, 1. 7:  L9.2.6  122E-1 9 ,4. 305 3 E-19 ,8.66l6E-19, 1.4605E-1M 


XIN  »  ABSIXINI 


W  -  Y0UT 


END 


DIMENSION  MI(33)«XI 1331 


REAL  N0*MU«MUE 


C72756. ,9026  7. ,104566. ,113 097..1 18391. ,133171. ,143  33  A. ,17352  4., 


C69807 

_ £56142 

2. *10 
30.  .6 

29840.* 13 50700. « 17314 50.* 2323430. *3037390. ,3970640. • 
524840./ 

OATA 

XI  /5. 

3 7656-21, 9. 1264 E-21,l .4096-20,2.15436-20,3.02626-20, 

Cl. 21266-19, 1. 75956-19,2.6 1226-1 9,4. 30536-19*8 .66166-19, 1. 4605E-18 , 


C6.167AE-17, 7. 75  L8E- 17,1 . 1  396-15 *1.53926 -16, 2. 0249 E-16, 2 . 8 90 76-16, 

%  AC  C  A  r  a  .  C  r  .  A  .  .  A  A% 


SON  -1.0 


ME  -  ABSIME) 


FLO  -  A8S(X0UT*N0I*SGN 

IF 


ENO 


32 


_ LUNCTlflN  QLtXlll _ 

C  THIS  FUNCTION  CALCULATES  THE  LONCITUOINAL  DIFFUSION  COEFFICIENT 

£ _ FOR_P  ARTICULAR  VALUES  QF  E/N  IGRIO  VALUES  FOR  QLHU  ARE _ 

C  FROM  LOMKE  AND  PARKER) 

_ Q-l-flE-NS-IQM  QLMU (  38> , Alt  38) _ 

COHMON/CONST/S,NO, EPS  I .ALPHA, GAMMA, ECHO, MU, HUE, K,E FLO 

_ E.XIERMAL  W,  INTKP _ 

REAL  NO , MU  *  HUE 

£ _ DATA  FROM  I  OH  K£_JlNQ_PAR_K  E  R.PRYS  RFY.lBl.P  307  AND  LOMKF  AND  DAV  I  l 

C  ,J  APPL  PH  YS ,4  d  1 12  ) , DEC  77*  P4995 

- Q  AT.A  QLflU /.  00-69.242^.110 fa 92 hi* l *JLUJUt±*+  Ofl-fl  75 1  *  .  QI 52 lit* . Q 31Q  7* _ 

C.06U^,  .  12*01  ,  .20  71  39,.  2  460 53,  •  19269,  .11 1456,  .082 39,. 08008 5, 

_ C^llft_49l2..1l  14 5h..  IT  1409,.  >672  49  .  .  36  1  ft  3 .  . 484 88 . . 4 7 1 3 89. . 647463 .  _ 

C. 72*177, .867117, .90108b* 

- Cl.  20595.  I.  SBQ  18,1.98908.2 .3177*  .2.759  62  .9  .  L  28  ft  .  3.  46  S3 . 3 .5  77  7, _ 

C3. 6314, 3. 70162, 3. 9 3 726, 4. 2  4172,4.6186/ 

- DAIA.  -X1/3.2643  7E-2  1.6.  0  42  2  7E-2  1.1.  11096C-23 . 1 . 42495E-20 . _ 

Cl.87932E-20.2.20343E-20,2.53922E-20»2.9l608E-20,3.24638E-20, 

- Cl .  71 5 15  £  -2  Q ,  4 .  b  2  7  2  E.r2  Q  O  rfQh03£rZQ^J, JLA  J12&- 2  Q  ,..9.-58  738  £  r2.Q.» _ 

C1.61065E-19, 3.2062 7E-19, 8 . 505 S£>19, 2. 1 94325-13,4. 1 19036-18, 

- C.7...919  55£  -la  ,1 .0  2  351.E.-17  «JL  *3flZl£  -17  «X.  65J4BL-1 1 , 2,  Zl>3bE.rX7ji _ 

C2.93833E-17, 

- C-3-.2Sb.lE  -1 7  « 1.  85  faf>£  r  1 7 . 4 .2  S2UE=JL  7,4.697.9-£-.17jl5 . 4953E-1 7  »  ha  54JLBE-17 

C8. 80446-17,1 « 1 5296-16,4 .09266- 16  ,6.021 4E-1 6,9. 54  3  3£— 16,  1.58 16E— 15i 
_ C2.S-57U.-15/ _ i _ : _ 

KIN  -  A  8  S ( X I N ) 

- CALL  1NTRJ*  IXI,DLMU,38,2«X  INaYQUT) _ 

DL  •  A8S ( YOUT* W  <  XIN)/(XIN*N01) 


IF  <Y<4)  .LT.0.0  I  SGN  -1.0 


HUE  »  WIY(4)/N0I/Y<4) 


.*ri:Lri  nnroinri  him  ci  ti:  u  i 


Pnl 


P0(1.2)  -  ALPHA  *A8S ( Wl V  <  4  ) /NO ) 1  -  GAMNA*Y( 3 )/Y( 4) /MU 
.it  ■  -r.  amma*/  <  ?  i/vi  <.  »  /mu 


POU^t  -  GAHMA*Y<2 )*Y< 3) /HU/Y( 4)**2  ♦  Yf  21  ♦  ALPHA  ♦HUE 
pnri.si  •  n.n 


P  0  (  2  »  1  >  -  /OL  IY(4)/N0) 


P0<2*3>  -  0.0 


PDC2.51  - 


P0( 3.2)  -  ALPHA *ABS t HIV  14  I /NO) )  *  GANNA *Y (3  I /Yf  A) /HU 

pnn.^i  «  -  r.iNNA*¥(?}/y(4)/ttu 


P0I3.4)  -  GAMrtA*Vt2)*YC3l /NU/VI4>**2  ♦  VI 2) *ALPHA*HUE 

pnn.si  a 


P0I4.1)  - 


P0I4.3)  -  ECHG/EPSI/V«4»/MU 


P0I5.2  )  - 


P0CS.4)  -  - 


PARTIAL  DERIVATIVES  ARE  NEGATED  TO  ALLOW  FOR  INTEGRATION 


PD (  1.2)  -  ?D(l«2>M-la»**« 


|J»ICTT4| 

EXXQE3 

imwr 

MRMNMBM 

MMM 

|J*XfrTfc» 


rm 

□E 

Pr; 


PD( 2*2 ) 


PDC3.3) 

pnn.4i 


P 0 <  4* 2 ) 


P0I4.4) 


r^'iprrAUDnuj. 


P0t3*3)*<-1.)**K 
*(-1.1*4 


P0C4.2)*C-1.I**K 


yatrviU! 


•  P0<4.4)*(-1.)**K 


I 


ENO 


51  MENS  ItJN  T I  4501,  JEU50I  ,N£|4531  ,NI<430»  .JPI450)  ,E<4501,  VI 4501 
CPTUSOl  .PJht  460  1  ,PNE  <4501  .PJP  <4501  ,  PE  1450),  PNH  fL.XQJ.iiiUiO.U_ 
ClPAKIlbOl 

_ COM  NO  NidJMlXii  NO, EPSli  ALPHA, GAMMA,  f  CHQ , MU, MUE , K , E F LQ _ 

COMMON/PLT/T  ,  JE  ,NE  ,  Nl,  JP  ,  E  , V,  VN  l  N,  l)*SH T  H,D  A  SHTH,  M  A  ,  MAS  *D 


M  •  MA 

-TALU_£QM?_R5 _ 

EKCK  -  EKICK 
_£jUJL_  jMLNPLJ  M  JCJ _ 

CALL  PAGE! 11.30, b. 5) 


CALL  8ASALFI 6HL/CSrO) 
CALL  TITLE!"  ‘SO, 


“•51 


■in  ■  m  numrimti  rw  wqsmi  ttiinwn-wn  i 

ryw^ffcnnnjnpncw[y4w|ti#f^ 

■QQvnni 


CALL  LINES! "!J1 1  ION  FLUX $?*  1PAK .4 1 
CALL  LINES  l"!V>  VOLTAGES  ",IPAK, 51 


CALL  LINES  !"!E>  ELECTRIC  FIELDS". IPAK.61 
NQCH 


35 


V  -  6.75 


CALL 

STRTPTIA.at 

CONNPTCAiY) 

NLi 

CALL 

CONNPTI Z,Y1 

_ CALL 

CQMNPT ( Z  ,d  1 

CALL 

CONNPT ( A,8 ) 

: 

►  i 

-  ■  i  ■—  ■  ■ 

C 100, < 
_ SIQ.JU 

►  1.0, -1.01 
>0  MM  -  1,MK 

REVERSE  POSITIONS  OF  BOUNDARY  ANO  CATHOQE,  SHIFT  ZERO  IN 


DIMENSION  FROM  V/CH  TO  KV/CN 


UUiACTiKtJlLJiH 

PiHr* 


I A  I  111  InVJBJH  k  I 


350  CONTINUE 
MK1  ■  MK  ♦ 


NASi  •  MAS  -  l 
MM.  -  MK 


C 

FILL  REMAINDER  OF  PRINT  ARRAYS  R  ITHQUT  INVERTING  ANY  POSITIONS, 

c 

BUT  SHIFTING  THE  POTENTIAL  2ER3  AND  CHANGING  V/CN  TO  KV/CH 
PT(MM)  -  TIMM)  ♦  OKSHTH 

PJE(MM)  -  JEtMMI 

PNEIMM)  •  NE(MM) 

P  N I ( MM )  «  Nil  MM  ) 

PJPIMM)  »  J  P  (  M  M  ) 

PEIMM)  *  £ ( MM  )  /  1000 . 

PV(MN)  -  ( V  <  MM )  -  V  M IN  1 

PEtMM)  -  E ( MM ) / 100D 


muvlem 


IflQJ 

FALL  (VHIN)  ANO  BY  THE  VOLTAGE  DROP  ACROSS  THE  POSITIVE 

III 

00 


E  FLO* ( 0  -  DASHTH  -  OKSHTHI) 


450  CONTINUE 


RR  ITE  <3,101  NAS 


DO  20  MM 


10  format  <  IX,  I5,7UPEL3.5,2X)> 


CALL  CURVE  <PT,PNE,M,0» 


CALL  CURVE  <PT,PNI,M,0* 

— CALL  .YgKAX-SI-J.0El7.1.ei7t*3.0E17«6.0.*t.»E«(J>»  -  1/CM2/SEO", 

C-100, 8.0,0,01 

_ CALL.  .QAJSH _ __ _ 

CALL  SPLINE 

_ CALL  CUJLY&  .  <  P  T  ,  P  JP  .  m  ,o  1 _ _ _ 

CALL  REStT<"OASHMI 

— CAJJL. -CURVE...  1 EU2. JLt » .01 _ 

CALL  YGRAX$<-3.00, L .00,  3. 00,6.  ,  "(VIOLTAGE  -  K ( V I %«*-100 ,8 . 7 ,0.0  I 

— CALL-C-UNPlil _ _  ■  _ 

CALL  SPLINE 

_ CALLJLURVX  <PT,PV.M.Q)  •  _ _ 

JALL  YGRAXS(-100.  ,  2  5.  ,1 00  • ,6. «H  (E I  (FI  I  ELD  -  KC  V|  /CM",  IDO  ,-.7,  0 

_ CALL— C..HN  Dib _ ■  _ _ 

CALL  SPLINE 

.LALL— CJJJj  VI.  IPT,PE>M,C>) _ ' _ 

CALL  ENOPL ( MK I 

_ CALL  .0J3NEPJ _ _ _ ; _ 

RETURN 


FROGRA  M  CASEIII<INPUT,CUTPUT,TAP£3  -  Cl  .'PUT) 

01  H£  NSIQN  Y  17)  .T<EJ  ■  ) , JE(b  )  t NE (6 - C )  .NI (6U" >  . J P ( 6J S )  . _ 

l  £  (  fc  -  . )  *  V  (  €  "  ’ }  iMIlEi") 

IQMHCN/CONST/StM'.  .  £°SI  .A  LPrA  .uaMHA  ,  ICHS  «  MU.  *UE  » K » EFL 3  .  NOX  »N£X  *p; 


CCMMCN/PL-  / / ♦ JE .NE  ,NI,  JP  .E  .V.VMIN.O  KSr.  rH ♦ lAShTH , *A  .MAS  ,C  .Nil  .  JI . 
CSKICK 


COMMON/GSARJ/  hUSED » NGUSED , NSTEP . NFE ,  N JE 
EXTERNAL  OTFFUN.PE CERtf .W.FLC 


REAL  M£,N; . NP «  MJ »M UE «MG . JE ♦ N£ , Nl « JP  » JI .Nil. NUN . NN ,  NOX . NEX 


-ISS°LA  PLOT  OPTION  (1-YE5  ^-NOl 
CPLT  *  l.r 


SET  CONSTANTS  FOP  GE*R  PROGRAM 


NUMBER  OF  0 IFF  £CUNS  IN  SYSTEM 
N  = 


LCCAL  EP9CR  TOLERANCE  PARAMETER 
EFS  *  ,jli 


GEAP  METHOD 

Mc  = 


PERMIT  I VI TY  OF  GAS  <C0ULCM8/VCLT/CM> 
EPS  I  =  8.95E-1A 

ELECTRON  mass  “  'V 

ME  =  .StlE-2? 

EOLTZMANN  CONSTANT  ‘  •  ' 

EK  =  1 . 39 E  - 1 6 


INITIAL  DISTANCE  STEP  SIZE 
T STEP  a  l.S-4  ' 

TSTEP  =  ?  •  £-6 


/STEP  s  5. £-5 . 

index  fqp  first  call  tc  gear 


DEX  3  1 
CAf  PRESSURE  (0YNE/CM2) 


c 


TG  *  273. 

ELECTRON  T£»pg3ATURI  f =V) _ 

TE  =  '.''4 

G _ NUH3E*  OSNSITy  CF  GAS _ 

N*J  =  PRES'HK/TG 

C _ NUHE^R  DENSITY  OF  AT~acHlNG  GAS _ 

NOX  =  i.*(i.  12E6/76.  .)  /3K/7G 

_ NCY  =  l.£-?;«NC;< _ 

G  MOLEv.ui.AR  U  £ I^HT  OF  GAS  (GM) 

_ =  4<~ .* i . q7£-27 _ 

C  FCSITIV£~ICN  MOBILITY  (0M2/SEC/V) 

_ HU  =  U.Vill/N1: _ _ _ _ 

G  NEGATIVE  ICN  MOBILITY 

_ »UN  a  -  *10' _ 

C  REC0K3INA-t'Icn  COEcFICIENT  CCM3/SEC) 

_ GAMMA  =  UUE-7 _ 

C  IONIZATION  SOURCE  (ICN  PAIRS/CM3/SEC) 


_ S  *  3.  6£lfe _ _ _ _ 

C  NET  CURR-:NT  IN  DISCHARGE  ( AMPS/CM2/SEC) 


C _ ELECTION  N  UM  5  £  E  C£  N SITY  AT  STA  FT  < 1/C  M3) _ 

WRITE  (3.E99)  PG  AS  »TG  t  NOX  f  NO  *S  » GAMMA 
_ E  ~  1. " 

YG(2)  =  SGRT(S/GAMM*/U.*B)  > 

_ <0(1)  =  C/ZCHG _ _ _ _ _ 

EF  =  FLC(Y„<2> /Y3( 1) ) 

_ RA  a  OATCHUP/NO) _ 

CN3  *  l.E-**NCX 

_ RO  a  ?.E«1C  _ 

AN  =  5  *E-2*NCX 

_ CO  12  JJ  =  1,20 _  •  _ 

EF  =  eTE'(Y«UJ /Y  X2i) 

_ GH  =  r»AHA(£F/NC) 

KA  =  0 Af CH  <  Ee/Nw)  ■  •  ' '  ' 

_ E  »  NQX»RA/(FQ»OND  +  Vn(2)»RA> 

YG  (2)  =  i GRT  (S /GM/  U  •♦2) >  '  '•  ' 

_ AN  =  5.e-2*NCX  -  B«YC<2)  _ 

IF  (  AN »L7  •  '■•71  AN  «  T  T~*  '• 

EF  *  FLO  (Y  (iXY  3<  2)  )  - 

YOU)  =  O/ECHG  *  MU*?F  *Y&  ITT  ♦  2  4  T  (2) 

_ WRITE  (3,SS>9)  YOU)  tYO  (2)»GM>£F»BtRAtAK _ 

12  continue  " 

15  CONTINUE 

FORMAT  UX  »7(l#£iy  .9>)  ~  r  ’  ’ 

C  POSITIVE  ICN  NUMGE R  DENSITY  AT  START  U/CM3) 

C  NEGATIVE  ICN  NUMBER  C&NSITY  U./CM3) 

NN  a  3*V  i  (  2 )  ‘  "  *  • 

WE  *  YG(1) /Y: (2) 

C - 1  HUS' I- AT r-£  FIET3 - - - - - 1 - * - 

EFLO  *  FLO (Wf) 

C - c rRU^TTTr£LC  'TQ'  8g6lN  TfFrgCR'XTIS'N  TCWABP ’JATHCPi: - 

YC  (4)  =  EFLO  ♦  .Cl’EFUJ 

*u£---=  "  1  1 

G  POSITIVE  ICN  CURRENT  DENSITY  IN  POSITIVE  COLUMN  <  1/C  M2 /SEC) 
foT3)  a  hu*Tf Ll^Nf5  .  r 


o|o  o  rt  o  o 


C 

kEChTIVE  ICN'  CURRENT  OEN3I7Y  IN  P03  IT  I « £  COLUHN  tl/CM2/3£0) 
r;(8)  =  NN*MUN*£FLC 

/■% 

'J 

O’AL  FcUX 

C  =  iYH!)  ♦  Y r.  (b)  -  Y2I3))*ECHG 

c 

D  ETERil  IN  £  ELECTRON  TEHFERhVURE  FCR  THIS  FIELD 

TE  =  T  £MC  (  E^LC/N  )) 

.v A i  E  Cf  ICNI7A*,ICN  CF  f-ETAST  A8LES 

RIP  =  5.  ;E-e#SO?T (TE) *E<P(-h.15/TE) 

c 

NUMBER  D £ N E  I  Y  OF  HE IAS  74 9L EC 

N£x  =  (l.-.*NO  -  YD(2)/HU/Y.I4))*i.*EXF<-il.fc5/TE) 

r 

FOTEN;  IAL  WCLTS)  AT  START  SET  AR3I*RARILY  TC  ZERO 

Y  0  1 5 )  =  J  .  • 

R  IT  £  <3,CyJ)  r  Hi  )  »Ta  <2)  »  vf.  <  2)  »Yu<  4i  *Y  J  1 6>  »  NN#NP 

KRIT£  43,999)  HUE, WE ,EFlC, TE,R£X ,RIH.N£X 

c 

FI-riT  CU T  PUT  pcint  DESIRED 

TOUT  =  r:  ♦  i.E-5 

•’RITE  (3,ij)*) 

WPITE  CJ.i'Cl) 

rtRiTc  {  3  »  i r*  3  9) 

ICO  9 

FORMAT  ( 1 2 X“V**l2X“  ALFH  A**iiX **HU£“13X **W**14X**DL**1,?X**CL/NUE*V) 

WRITE  (3,1':.JS> 

HRITE  to,lJ02>  TQ*Yj(i)  ,YJ(3)»YC(6),C,Y0(2)  »NP»NN*Y3(4),YI(5) 

C 

STORE  PARAMETER  VALUES  AT  8CUNCARY 

’55 

r*ii  »  r.d) 

Y-)2j  *  Y„<2> 

Y0  33  a  YT  (  3) 

YO  43  3  YH4J 

YC5C  =  Y3<5) 

Y 063  *  YM6> 

Y070  a  NN 

C 

►  IS  THE  NUH8ER  OF  THE  INTEGRATION  PCIKT 

H  =  0 

C 

set  In  it  ilTTIiutt "  6?~h S.  "An 6  'J*  :  * 

Y93  =  +1.J 

• 

C 

Y02  *  ♦i.O 

SET  COMPARISON  VALUE  FCR  VOLTAGE 

VMIN  s  i  J  «j  1 « 


Vow  “  wit  h  K  E(i  1,"  (fiThdCI"-  YTOTTWan?  So  2* 

INTEGRATE  TC  ANODE  SHEATH  BOUNDARY  -  LATER*  WITH  K  EQ  4, 

IN7 EG »AT  £  FfcCH  mNODS  SHEATH  d<$UNOA*Y  10  ANdOI 
<K  EC  OHO  NUN5cR  NEGATES  DIFFERENTIAL  EQUATIONS  AND  PARTIAL 
CEfilVAl  I7tS  1C  ALLOW  IKTEURAv1CN  &CV  Tfi  CAtrfOOgi 
K  =  1 

3C*  CONTINUE  ■ 

C  - 

c  sSn  us-o  Lat*R  to  *In6  o$sItI<Jn  AT  nmICR  nE  65E5_"Tfl’ ZPftO  ~ 

_ SGN  =  i.’l  _ _ 


V 

15fl«  eORHAT  (iul  *3X**CI5T  ANC£**SX**  JE**liA***l»**llX**J-**5X*"CURRENT  OENSITY‘%*  , 

- g^N^iix «~i  1  x*'n - ,T^\wr  ristC^VOCASn - : - 

lCfll  ^CRHAT  (  6X**OH**3X“,l/CM2/S  EC’44X“1/CH2/$£C“4Y**1/"N2/SEC**4X**A*'PS/CK«J” 

- fTIVC>j,t  6  <*  I7C  V  !  M  BPTTCT  3"oVMV  a:TS7CV"7r,v  CUTS"  7/1 1 - 

1C05  FQRHAT  (lX^CATHOCE  SHEATH  8GUN0ARY**) 

z - - - - - 


4 o  CONTINUE 


o  o  O  O 


SET-IN  LOOP  TC  iOL Vt  uE  sy:- 
BETWEEN  *;NCDE  ANT  CAThOOE 


AT  VARIOUS  POINTS 


0  «Y  .»T  OUT  «  EPS  tWF*  INCcX) 


,  =  <Yi.ll)  ♦  Y.(F)  -  Yi,(3)i*SLrtS 

N  =  -Y  j  (  Ei  /YU  /Y  ;  <  O 


a  Y'  <  J)  /VU/V3 
IFT  =  W I Y  j ( 4) /N 3  ) 


=  OIFF/YUt 

=  T  EM° (Yj(4) /n2 ) 


..-?*SQPT  <TE)*tXP(-4.l5/TE) 

:’*Ni  -  Yr,(2)/'1U/Y-(4<>*l.*FXP<-11.65/TE) 


kritc  <3«i:ce>  m »a lpha t hue ♦ dr ift» oiff »te »nex*ri* 
8  format  <12X,I!,9X»7<iOt9.1»6X)) 


RITE  <3  » 1  •  1 2)  TGUT,Y  J  (1)  ,Y3<3>  ,Y01S> *C • YO < 2) ,N? tNN,YC <4) ♦V9<5> 
1302  PCcMa  <9<1PE12.4«1X>  «lFEl2.4) 


110  3  FORMAT  <//2bH  EFRCR  RETURN  WITH  INOEX  «*I3//> 
40  TO  20", 


13C  IF<K.£a.lJ  GC  TO  35-3 
:f<K.EQ.2>  go  TO  375 


"  IF  CK.EQ.4I  CO  TO  385 
253  CONTINUE 


TOUT 

=  r:<i) 


JP (M)  «  y:<3) 
E(N)  =  Y3(4) 


V  <  M)  *  Y  j  <  5  ) 
JI  <  M)  »  Y  n  (  6) 


NII(M)  *  NN 

FINQ  ilNIMUM  IN  POTENTIAL 


*NIN  «  A'lINKtfKlN,  V<NI 
OJPJE  »  l  AM2*YC<3)  ♦  V3ti) 


NE( M)  *  Y  3  <  2) 
Nl  (  M)  a  NJ 


IF<  ( 2«*CJc JE-GJPJE1)  .LT,  a.C 
GJPJEi  a  GJFJE 


C.*«L  i  •  1»E- 
IF  <<?.**'  <2)-Y02)  .LT.C.O)  ,TSt*P 

iTu^N  riNUi - : - : — 

YCl  a  Yu< 1) 


Y03  *  Y  0  (  3) 


•  TSTSP/5 


YJh  * 

Y  C  (  4i 

f  >05  « 

YC <  5) 

|  T0UT  * 

-  TO  IT  "i  ' 

T5TTP  ” 

140 

I/MIN  =  A  M I N 1  ( V  H I N  *  V  <  M 
IF  (YU  <^)  .L-.6  .Cl  GO  TO  293 


IF  ( (2  .*Yc ( 2) “Y02I .LT.C. 

IF  (TST£O.uT.l.£-a*0)  GO  TO  393 


Yd 1  *  YJfl) 
Y02  a  YD  <  2) 


s  Yfl  <  5) 

TOUT  *  TOU“  ♦  TSTEP 


GO  TO  4C 
390  CONTINUE 


►  A  »  M 

OASHTH  =  TOUT 


CKSH1M  = 

HK  IS  THE  INDEX  NUMBER  OF  THE  CMHOOE 


MK  *  M 

KRITE  (3.1JU6) 


<11  *  YO 
Y0C2)  *  YC21 
TTO>  a  Vu'TT 


FERTUR8  ?  FIELC  SLIGHTLY  TO  BEGIN  INTEGRATION  THRU  FOS  CCLUHN 


nn  s  y  :*r* 

N p  *  V  •;  <  3 )  /MU/EFLO 


»  uc 

a  .  CUT  i 


CESIREO 


ICUT  *  0.0  ♦  1 « • 
INC  EX  =  1 


WRITE  13, 1  .  • 
WRITE  <3,1  *1) 


WRITE  (3,  lJiE) 
WRITE  <3,1  *5} 


WRITE  <3,13  2)  T. 
SC  TC  3'j 


‘•ii  CONTINUE 
rt  a  M  4  1 


INDEX  NUMBER  CF  ANCO 
MS  =  M 


,YC(1>  ,ra  13)  ,Y0<6>  »U  *Y  J  (2)  »N?»NN,Y(j  14)  ,  Y0  (?) 


SHEA TH 


=  CC 


£  I  M)  a  YP  C  4)  *  ♦  .1 
R  A  =  n.UCH(Y'  <  A>/N0) 


G^  a  GAMA  (Y  j(4)/NU) 
YC  12)  *  Y  0  2 3 


AN  a  5.E-2*N0X 
e  =  NOX*RA/(PD*ONU  ♦  Y0(2)*RA> 


Y012)  a  S CRT (S/GM/C1.48)) 

9  =  NOX*R  A  /  1  RQ*QNC  ♦  YM2)*RA> 


Y3  12)  a  TCs'MS/GN/ U.*3>) 

a  =  nox»ra/(rq«ono  ♦  v;i2>*rai 


Nt(H)  =  YO 1 2)  a  SORT (S/GH/li.* 

JE1H)  a  YO  (1)  a  C/ECHG  ♦  Y ?  1 21 *«L*Y  „  U> * ( l.»Z.*B) 


NI(H)  a  NP  a  YC  <  2)  *  1 1.  ♦  9) 
JP(M)  a  va<3>  *  MU*YC(4)*NP 


C  *  !VuH 

WRITE  <3,11.1£)  Y3Cl>  ,YC  (2)  ,Y£  <3)  ,V315>  ,V)16)  ,NN 


FORMAT  (1)T,T(1F£i5  .7) 

TIM)  ASSIGNEE  A  VALUE  CF  ZERO  AT  HAS  S*NCS  ACTUAL  POSITION  OF 


TART  CF  AN 
MM)  a  TO.  *  9.<i 


mm 

nM^Si  :f  ■?:!  im  i-  j  -*  j ;; ;  *  xa  m  j 


TSTEP  a  5.E-6 


TSTEP  a  ?.£-4 


EPS  • 


N 

WRITE  (3,1000) 


WRITE  (3, IP  11) 
~ill  FORM  A 


sr 

.v» 


.  RITE  {3.  1C. 2)  T".Y3 (1) .1 J 13) *Y0  C6>  *C.Y3(2)  .W.NN.Yj  (4)  .Y3C5) 
GO  TC  32 


1 


o 

C 


21  £  HR  I1 
1J*-  CCN1 

CALL  EXIT 

_ iNO _ 


21X.I5.14H  P  EVALUATION!/ 
21X.IS.14H  J  EVALUATICNS///I 


SUBROUTINE  OIFFUNfN.T.Y .YOOT) 

DIMENSION  Y  { U )  » V  00  T ( Nl _  - 


EXTERNAL  U»CL 

COHHCN/CONS  T/S  .  N  .•  «  EPSI  *  A  L°H  A  .G  AH  EA.ECHG  »MU  .  KUc  . K  .EFLS »  NOX  »N£X 


IF  (i.^6i-i5*N-/A9S(Yl4n.tT.4-t.)  iO  TO  15 
£(. Ph A  =  Nr»2.FF-l7*£XP(-l.  -♦flE-l5*Nrf/AESIVt4))) 


AN  =  5,£-2*N0X  -  A  5S  JY  (  6> /MU/X  C  4)  ) 

if  (t.N.r.1.;)  an  »  'i.c  _ _ 

FT  A  s  CA*  :hH  C4)/NC)  *AN/m9S  CWIY(  4>/NC)) 

S A NNA  =  GAMAtY  f  4)/  NO ) _ _____ 

INj  =  l.£-J*NCX 

FD  a  3  •£•  1  ?  _ _____ 


NN  =  -Y  16)  /Mli/Y(4> 

T£  =  T  EMP  ( Y  <** )  /NJ)  _  ~-  -- 


RIM  =  5.v.0E-i*SCRMTE)«EXP<-4.15/Ti) 

“  *  ft.7F7»SCoriTF)-*7.E.ie»|2«»T£»ll«fc5)  *£XP( -11 «65/TE> 


SIGN  CONVENTION  IS  THAT  THE  DRIFT  VELOCITY,  W,  IS  POSITIVE _ 


WHEN  THE  ELECTRIC  FIELC*  £>  IS  NEGATIVE 

SGN  a-i.O  ■ _ 


IF  CY  <  4)  «LT  •  D  SGN  »1.3 

YOOT <1 )  =  S+Tf2l«f  tALPHA-ETA)*A8SfWCY<4)/N0>)-SAHfrA»V<3)/HU/yf4) 


C<RC*NN*CN:  ♦RIM*Y<  2)*NEX 
Y00T<2)  =  <Y(2)«W<Y  <4)  /NO  /SQM  Y  (1)  )  /OL  CY  t4)  /NO  ) 


N-) l-GAMHA^V 


Y00T<5> 


i/TA  j-'M.T'TE  i  »M  i 


CJ«VOX  POSITIVE  WHEN  STARTING  IN1E6RATIQA 


IF  ( (T  «EQ*  J  • 
GC  TO  4 


3w  Y30TU)  -  . 

YOOT (2)  *  0  ,c.  : 


jYCOTU)  *  J.C  •  -  »  •  ■  ... 

Y00T44)  «  5.E2 


YC0V<5>  =  -EFL 
YOOT  <6)  = 


C  S 


t  FpSniVE  COLUn 

57  YOOT(t)  *  YOOT <21  «  Y0CT<3)  *  VOOHM  *  0*9 


73  CONTINUE 
VCCT(l)  ; 


YOOT (2)  * 
•'COT  (  ?)  * 


YDOK4)  a 
YOOT (5)  * 


♦  l.£8 


-l.ES 

YCOT(l) 


♦  1.E3 
-EFLC 


CIFF  EONS  A^E  NEGATca  to  ftLL0W  FOR  INTEGRATION 


YOOT (1) 


YOOT (2) 
YOOT  <  3) 


YDCT(a) 
YOOT (5 ) 


YOOT <6  > 
RETURN 


Y COT  <!)*  (-!•) 


Y0CT(2)*C-1.)**K 

YGOT(3>M-i.)**K 


YCCT 

YOCT (51* 


YOCT(6)*(-l.) **< 


FUNCiICN  P  AMA (X INI 
GAMA  a  8.81E-7 


FUNCTION  CLNfXlN) 
OLN  =  .15FS5 


FUNCTION  3ATCMYIN) 
ATTACHMENT  PATE  IN  CM3/SEC 


CIMENSICN  YII23) ,A  T 1 23 > 


T5E-!2,2.1C91E-l2.I.39t-l2» 

C4.93gfc£-1?  .9.27 9E- 12,1. F^gE-ll,  2.1 -,83 f-1 1,2. 475 lE-ll, _ _ 

C2.77  52E-li.2.i261E-ll,2.*,<lE-ll,  2. 672E-11 , 3. i26 7 *-11 , 3 .*7*2£-ll. 


YIN  s  A 3S  (YIN/1  .£-17) 

CALL  INTPPfYI,AT,23,2,YIN,ATCUTI 


CATCH  =  At  OU 
RETURN 


£N 


FUNCTION  7£MP<XNIN> 
CIMENSICN  XN(ie) ,T<l€> 


DATA  XN/.j  ji,. Cl,. 1,1. ,2. 
CeC.,lTC.  ,113c./ 


>NI  N  =  A9S<XMN/l.£-l7) 

CALL  INTR°(XN,T » 16  ,2»XNIN»T< 


IF  (XNIN.L7  ..Cl)  T CUT  *  • 
TEMP  *  .C4 


f EMP  a  TOUT 

RETURN  _ _ 


£NO 


FUNCTION  W  ( X I N ) 

TH 13  FUNCTION  CALCULATES  ELECTRON  ORIFT  VELOCITY  tCM/S 


PARTICUL^o  V^LLcS  CF  c/N  <<»RIC  VALUES  FOR  CRIFT  VLLOC 
pc  ON  €N3ELHA;t  ANO  PHEL°S> 


Cl  ME  NS  TON  W I ( 32) tX 1(23) 

COMMON  /CON'S  7  /S  » N  . ,  E^SI , ALPHA  ,GANH A*  EOHG  ,NU»NU£  , K , EFLu , 


SEAL  N', MU, HUE 

EA.'A  F;-CM  ENGELHAR^T  ANU  PHELPS»*HYS  REV*  <133)  , JAN  64* 
CAT  A  Wl/3169.2  ,4894.  2,7657. 9,l3:  22.  *2229  3  **  3  32757  *$369 
C 72  756.  ,9.26  7. , 1145 66**11 3j 97., 11 8 341.  ,  2 2 3171 . ,  14 0 224 . , 


C 195621. ,2  22429. ,25  .227. ,291743,, 247516. ,39)211.  , 45154 
C  6 3 3 C 72 .  ,  1  ‘  246 4 ;  .  ,1  35»  7.',  i, .“,  17 2145  2.  ,23234 39.  ,30  3  7  3 EL  •  ,3 


Cp6 1423C . ,6^2-84 J./ 

_ 0  AT  A  X  1/5. 3765  F.»  21  ,9  .1  264£ -21 ,1  .409E-2., , 2.1543E-24,  3  .  9 

c3.3til*,£-2'',4.2626E-2C,4. 263  SE-  24,4.811  7E-2C  ,6.0  6  74E-2C. 
C  1.2126E-19.l.7595E-19,2.6122E-19,4.3ii53£-19,3.66ieE-19 


EC)  FOR 


IT Y  ARE 


0 , NQX ,N£X 


4, PA377 
F9X775 36917  . 
.,172524., 


44565. , 
,3971640.  , 


.9262E-2L, 
20.3.4T72E-2 
19,1. 462  5E“i  1 


CALL  I NTRP ( XI , HI ,3 3* 2,XI N, YCUT) 
W  =  YOUT 


RE*  URN 
NO 


FUNCflON  FLO(ME) 
CIMENSION  NIf 22) ,X 1(33) 


N/CONST/S ,NC , E 
REAL  NO.MU.MUS 


C  72756.  .P**  267. ,10  4566. ,1133  97., lie  39 


C 19 58 21 . , 2 2 2439 . , 25C2C7 . , 2917437, 
CE98072., 1C  2QS4C.,13597C0.,  1731 458. 


i  o  o 


u 


0 

V 


!♦ 

| 


* 

r 

t 


FUNCTION  9L<XlM 

this  FUNOVICN  CALC  LuAT£3  the  LONGITUDINAL  0IFFU5I0N  COEFFICIENT 
FOR  PAPTICCLAR  VALUES  CF  £/H  fiiPID  l  ALUES  FOR  DL  ►O  ARE 

_ F_ROH_L OMKc  AnL  °AP.K£P) _ 

CIM£NS ICN  CLKU  (33)  ,XI< 2d) 

CC-HCN/lQN-^/G  .Nj.  gPSI  .M-PKA^CAHWAtcCHi.  »wU«HU£»K«£FL3tNOXtN£XtP3 
EXTERNAL  rt.INTFP 

REAL  N  "  .MU  .MU  I _ _ _ _ 

CATA  FRO')  LCnKt  AN  0  FAR<EP,PHYS  REV.181.P  3J7  ANO  LOWKE  ANO  OAVI 

,J  A PPL  P i !  ?S, 43(12) «:£Q  T7 .  P4555 _ 

CATA  C  L  MU  /  .  J  . 1 9  2i2  » •  '0  €  92-»2  .  .  C  .  7  865  . .  8751  *  .S  15  234,  .C  31 T  7  . 

C  .  .16122  4. «  12  4  1..2C  71*9. .24R..53. .15269.. 111456.. 

0.^491?,.  lil456,.l7i4J2,.267246.  .361 &3 ,. 43485, . 571369. .647463 . 


G.725 17  7,  tc67ll7.  .9  :i',d  _ 

Cl.2b695il.58  16,1.95?.  6.2.31 77 4.2.75 96 2.3.12 88 .3.4653  *2.57/7  * 

03.6314,3. 7  162 ,1.9 3726. *.2 4172. 4. 6156/ _ 

CATA  XI/3.26427E-21.6. : 422 75- 21 ♦ 1. 11 - ?6£-23 , 1. 424 E5E-2 J. 
C1.P7932E-21.2.23342E-2  .  .2.  5Z522E-2r.  . 2 . 316.. 8E-2T  .  3.246?  g£ -2  T. 

C  ?  •  7153 5£- 2  _■  ♦-.4o2-»2c-2u  .5 . 9 -6C  9E-2C  , 7. 84J32E-2 0  .9 .58738E-2-  • 
CJ_.  61 . 655-19 .2. 2362  7E-19. 6»  5F  55E-15 , 2. 1S432E-18 . 4 . 119^8  £-18  » 
07. 919545-13.1.  2 35  3 E - 1 7 , 1 .  5a  2l£-~17 . 1 . 6 53 <♦  8E - 17 , 2 . 26i 6E-1  7 , 


C2. 3383 35-17. _ _ 

C2.2961E-1" .3. 6563E-17 »h»  2521£-17  ,V.c975£-i7,5.485  2E-i7  ,^.3418E-1 
Cc.cJ 44 £-17.1.1  629F  -16.4.  C926E- 16 ,6.  w2l4£-:i6 .9. 54  33E-16  ,tV5616£  - 1 


C2.b47l£-15/ 

XI N  =  ABS(XIN) _  . 

CALL~TvT9P  <  XI , CL MU  •38*2*XZN«Y60n 
CL  =  ABS(YCUT*h(XIN)  7IXIN*N0>1 
PET  URN 


•*'V 


SUnRC'JT  I\E  FFDERVI  NfT»Y*PO»NCr) 

:iM£f^TON  ?ctN ' ;»n:j>,y (n>  _ _ 


EXTERNAL  CL.K 

CCMMCN/C'IN:>'’/S  » N'j  t  EPS  I  tALPHfl«G  AM(*A»£CHv-  «KUf  MU£tK  t  £FLC«NOXtNeX,P. 


SON  =- 1 . 1 

IF  (Yt4>  .LT.  ..  ')  SGNai.J  . _ _ 


FUE  =  W(M<«>/Nt)  /Y  J4) _ _ 


CHECK  ON  MAGNITUDE  CF  EY  PONSNT  TO  PREVENT  UNDERFLOW  CF  RESULT 
Ie  (i.46E-l5*N.?/A8S(Y{i»)  >.r*^.2£0.*K)  f-C  TC  lr  _ _ 


LPHA  =  Nr*2.9E-l7*EXP(-l,U3E-l5*N./ABS(YC4)i) 

C  TC  2j  _ 


1J  ALPMA  =  J.j 

2  J  CONTINUE  _ _ _ _ 


po<i,d  = 

PC(1,2)  =  ALPHA* 40 S(*(>  (4>/NL)>  -  G  AMPA*  Y  (  3)  /V  C  4) /HU _ 

F0(l,3>  =  -£APNA*Y  (2)/Y  (•*)  /KU 

_ FOCI, 4)  =  GAMMA*Y( 2)*Y<3) /HU/Vl4l«*2  ♦  Y ( 2) »ALPHA»MUE _ 

P  C (1 *5 )  =  j.i 
FD(2,1>  =  -l./CL  (Y  <4*/K3> 


PD<2»2)  =  SGN*H(Y<4)/NO)/OL(Y<4)/N9 
PC  <  2 »  3  J  = 


P0<2,4)  =  Y<2I*MU£/0L<Y(4)/NI)> 
PD (2  *5 )  a 


FC<3  »1)  a 

PDC3,2)  *  ALPKA*ABS(H<Y<*)/NC>)  -  GkMHm*V (3) /Y (4) /KU 


FC (3  *3  )  =  -  GAPMrt*YC2) /YW> /MU 

PD  <  3  »  4 )  =  r,AHNA*Yf  2)  *Y<3)/HU/Yt4>**2  ♦  Y  (  2)  *AL  PHA*MUE _ 


FO<3,5)  *  3.*? 

F0<4,1)  =  3.0  _ : _ ; _ ' 


PD(4,2>  a  -ECHG/EPSI 
PO<4,3)  *  ECHC/E»SI/Y(4>/MU 


PC  t  *♦  »4)  a  F0(4,2)*Y  (3)/KU/Y  CM* 
F0<4,5>  =  3.? 


PDt  5»1 )  *  * 
PD  1 5  »2)  = 


FCC5,5>  =  w. 

°APTI AL  JEe IVAT IVES  ARE  NEGATED  TO  AU.OH  FOR  INTEGRATION 
FROM  SHEATH  PCUNOAfiV  9AC K  TO  CWHWc  ’ 

PC  <1»2)  =  PC<1.2)*«-i.) **K 


PC  <  1 »4)  =  °D l 1 «4>* (•!•) **< 


PD  <  2  »1 )  *  PC<2,1>* (-I.)**K 
PC ( 2  >2 )  =  POC2.2** <-!.>**< 


SU3RCUTItlr 
Cl  ME  NS  r_OM_ 
CFT  ( VI)  tPJ 
J-IPAK  <1  5.  I_» 
CC'mMCN/CSN 

cov  4Cn/pl  r 

"C i<ICK 

c;al  j'tJp 
'  K  =  MA 
EKCK  =  £FL 
CALL  CCKPP 
CALL  *C-N=L 
C A LL  P  NU 1 1 
C_ALL  PhYSC 
”  CALL  R ASA L 

_ C  ALL_  w  1 X  A  L 

C  A  LL  ~  IT  LE 
C"<  S) IS  TAMC 
CALL  C  *OSS 
CALL_GHAF( 

CALL  M~£SGA 
CALL  HESS a 
CALL  BLNKi 
CALL  L  INES 
CALL  LINES 
CALL  LINES 


CALL  L INC S 
CALL  LINES 


CALL 

CALL 


L  IN  rS 


FLCT(MK) 

T  ( € '.  v  )  .NEI6C  I  .Nile  I  )  *JP<6CT)  . £ C6C g  I  t'M  63  3)  « 

r  (b  ;  V)  ,PNt  (b.  ’  I  *  PJP  <  6K  :)  ,F£(cJ  J)  ,PNl(6J  :)  ,PV  <6jJ>  • 

J 1 1 e ^  •Nll»6:c>  ♦  CNII(5.’'> 

3T/S.NJ,  E  P  S I  ,  A  L  P  H  m  ,  C »,  M »» *“,  E  C  K  u  ♦  ,  NUE  »  K  ,  E  C  »  n6  X  ,  NE  X  *  P  i 

✓  Tt  J-  ,NE  ♦  N  I»  J  P»E»V, VKIN??i<SH"H»CASHTH«MA,M  AS, D, Nil  ,JI, 

t  JIfNE«Ml,MI _ _ 

—  — «  1 

5 

(mk> _ _ _ 

i  1  •  j  1 3«  5 ) 

-  (1,5.1.  5) _ _ 

F  reKL/CSTu) 

F(triSTMNCAFD) _ 

I**  **  »ii » 

E  -  CMS”  ,1CG ,“<N) E, (N) -,(N)  ♦  -  l/CM3f,  U  3  ,  S  .  0  ,6. 3  ) 

ih.Im  3,-6.3£ll,2.CEll,+6.:Ell) _ 

3  <**(0  ATHOE£S*‘,iCJ  .-.5  t-35) 

G  (“(A)  NO  CEVtluj  f7.7f.5) _ _ 

(Z«n5i65io«  .5  il •  7 v I 

(**{N>£  N  UNEER  uENailTV  CX)  1*?  g  f,  1PAK,  II _ 

C\*N>  +  NUMEER  DENSITY*”, I®AK, 2) 

NUMEER  CENSITY  tX>li0l“*IPAK*3> 


<  “  <  J  >  £  ELE-TSON  FLUXS-,IPAX,*»» 
(”tj)+  PCS  ICN  FLUyS”.rPAK,5) 


ME"  ION  FLUX  tX>L-JT*”,IPAK,6) 
C*(V)  V  CL  TAC  ES**  » IP  AK  *  *) 


CALL  ST9T?r<A,«) 

CALL  '  CNNP^  <>  _ 

CALL  L  CNN eT (  7  ,  Y ) 

CALL  C  CNN° " (  Z  »  5 ) 


I 


CALL  CCNNF7<  A,<?> 

C*LL  M  ES^  AN  (**  ( F  >  IG  3. 


1  •.« v  t  ♦  1  ♦  L  ♦*!#') 
r.C  35'  MM  =  l,P< 
t  REVERSE  =0*3fTCNs  C*  SC 
potential  fpcm  boundary 
i  p£  f  s  i  on  frp*  V7C1TT0T 
MM  P  s  VK*1  -  *M 


FT  (MM)  =  KSc,Th  ■ 
FJE(Mil)  =  JE(KMM) 
FnE(  MMi  =  NtTHFMl 
FNI(MM)  *  M<  M**M) 


FJP<M)  = 

FEtM.M)  *  E(NMM) /l?  CO. 


fUOM  (Cl  URGENT  (OIENSITY  (PILASMAS”, 


;hunm|  my  u  vwiruuwi  <- 

TO  CATHOCC*  ANC  CHANGE  E  FIELD 
- - — 


FJI(PM)  =  lOUg  . 

ONII(HM)  *  Nil ( MMM) *1J? • 
FV  (MPf*  TvTFTFm)  -  VMM 


5i 


K  R I T  £  (  3,  in  )  VM,PT  (MM)  ,PJ£(MM)  *PNE  (MM)  »  FM  O'M)  ,  *  JP  (  MM)  »P£(MM)  tpv  (M  " 

»5»  CONTINUE _ _ _ 

KK1  =  MK  ♦  1 

MAS1  =  MAS  -  1  _  _ _ 


CC  4'.?  MM  =  MKi  ,  *A$1 

FILL  REM»*INCcR  CF  PRINT  ,-RRhYS  WITHOUT  INVERTING  ANY  POSITIONS. 
3HT  Sri IFT  T N-i  ThE  POTENTIAL  ZERC  AND  CHaNG I Nti  V/CM  TO  KV/".M 
FT  (MM)  =  T(MM)  ♦  CKShTH 


PJEIMM)  =  JE(MM) 
r  N£  ( MM )  =  NEfMM) 


F 


F£(MK>  =  E(MM)/i!jJn. 

FJIIKMJ  =  JI(M»)  »1  A-iL. _ _ _ - 

FNII  (MM)  =  Nil  (MM) 

PV(KM)  -  (V(MH)  -  V  M IN)  _ 

•* .ofT£(  3,  lv  )  MM  t  FT  (MM)  .PJEl.MM)  »PN£(  MM)  .FNI  O'M)  tPJP(  MM)  ♦PE(MM)  *PV  (M 

AL  F  CONTINUE  _ _ _  — 


CC  4511  MM  =  MAS  *  MA 

KPItE(3»1j)  VH»PT(MM)  ,PJ£(MM)  »PN£  (  MM)  »  PM  (  MM)  ,PJP(  MM)  .PE(MH)  .PV(M  ' 


FT ( MM)  *  T(MM)  ♦  0  -  OASHTM 

FJE(WM>  =  JE(MM)  _ •  . . . . 


F  NEC  MM)  =  N-.(M*> 

FN!  (  MM)  -  M<M*)  _ ' _ 


*  JP(MM) 

=  c(MM)/l  ??8. 


OLUMN  (WHICH  IS  LINEAR  IN  X) 

PV(MM)  =fV(M*)  -  V  MIN  -  zFLO*  fu  -  CASHTH 
C  ♦  C  EKCK ) *  <  2-D AS»1TH  -  ORSrlfli)  /£* 

<♦50  CONTINUE _ _ 

WRITE  (3.1  )  MK 

WRITE  <3.1  3 >  mas  _ _ _ 


WRITE  <3*1  )  MA 


-  CKSHTH)) 


WRITE C  3, 1  )  MM .  FT  (MKI  «°  J  E (  MM)  ,PN{ 
FORMAT  ClX»I5»7(iPtl2«5*  EX)) 


CALL  CURVE  (PT,PNE»M,3) 

CALL  PASM _ _ _ •  _  .  _ 

/CALL  CURVE  (P7»PNI,M,)>  _ 

CALL  VC-RhXS  (-2.1 £15.  .5 El 5.  ♦  2.fl£H.».Q  J)£.CJ)  ♦ 


KH) .FNICMM) .PJPCKM) fBECMM) tPVCMM 


1  /C  M2 /SECS 


C-1C0  .3  .u  » L,  •  >) 

CALL  CASH  _ _ 


CALL  SPLINE 

CALL  CURVE(PTtF JI.M.10) _ 

CALL  ffU KVE  (PT.P-/fc.M.fl)  * 

CALL  RESET  (“CASH")  _ . 

Zall  Curve  (pt.pjz.m.^J  ~ ~~~~ 

CALL  TtoPAXSC-3  .  1. 13.0  »3).C- .6.  »M  (V)  CLTA65  -  (V)  SM«-130  .3  *7»_Q  »0> 
CALL  'fHNud T 


CALL  YCR’.XS  (-4.0,2.,4.t  ,6.  ,M(E>  <F)  ItLC  -  <  (V) /C  •*.  J*'»  10  •-  ,  -  .  7,  Ct  . ) 

CALL  C! <N Jj>M _ 

CALL  S=>LIK£ 

CALL  CUg^f  (PT,P£,  :i _ 

CALL  END=>L<*<K) 

CALL  DCNEPL _ j _ 1 _ 

AETUPN 

:NC  _ _ 


X&£ L 


« 


f  •  • 
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